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In systems possessing spatial or dynamical symmetry breaking, Brownian motion combined with 
unbiased external input signals, deterministic or random, alike, can assist directed motion of 
particles at the submicron scales. In such cases, one speaks of "Brownian motors". In this review 
the constructive role of Brownian motion is exemplified for various physical and technological 
setups, which are inspired by the cell molecular machinery: working principles and characteristics 
of stylized devices are discussed to show how fluctuations, either thermal or extrinsic, can be 
used to control diffusive particle transport. Recent experimental demonstrations of this concept 
are surveyed with particular attention to transport in artificial, i.e. non-biological nanopores and 
optical traps, where single particle currents have been first measured. Much emphasis is given to 
two- and three-dimensional devices containing many interacting particles of one or more species; for 
this class of artificial motors, noise rectification results also from the interplay of particle Brownian 
motion and geometric constraints. Recently, selective control and optimization of the transport of 
interacting colloidal particles and magnetic vortices have been successfully achieved, thus leading 
to the new generation of microfluidic and superconducting devices presented hereby. The field has 
recently been enriched with impressive experimental achievements in building artificial Brownian 
motor devices that even operate within the quantum domain by harvesting quantum Brownian 
motion. Sundry akin topics include activities aimed at noise-assisted shuttling other degrees of 
freedom such as charge, spin or even heat and the assembly of chemical synthetic molecular motors. 
Our survey ends with a perspective for future roadways and potential new applications. 
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I. INTRODUCTION 

Over the past two decades advances in microscopy and 
microscale control have allowed scientists and engineers 
to delve into the workings of the biological matter. One 
century after Kelvin's death, today's researchers aim to 
explain how the engines of life actually work by stretching 



thermodynamics beyond its 19th-century limits (Haw 
[2007] ). 

Carnot realized that all engines transform energy from 
one form into another with a maximum possible efficiency 
that does not depend on the technology developed or 
on the fuel utilized, but rather on fundamental quanti- 
ties such as heat and temperature. Kelvin and Clau- 
sius came up with two "rules of the engine" , later known 
as the first two laws of thermodynamics. The first law 
states that energy cannot be destroyed or created but 
only transformed; the second law sets fundamental lim- 
itations to what energy transformation can achieve in 
practical terms. Just as the first law was centered on the 
notion of energy (from Greek for "work capability" ) , the 
second law revolved around the new concept of entropy (a 
term coined by Clausius also from Greek for "change ca- 
pability"). 1 When expressed in such terms, the second 
law states that entropy cannot decrease during any spon- 
taneous or natural process. Notably, within the whole, 
virtual manufactory of all natural processes the first law 
takes on the role of an account clerk, keeping track of all 
energy changes, while the second law takes on the role of 
the director, determining the direction and action of all 
processes. 

The fathers of thermodynamics developed their laws 
having in mind macroscopic systems that they could de- 
scribe in terms of state (i.e., average) quantities such as 



1 Rudolf Julius Emanuel Clausius throughout his life allegedly pre- 
ferred the German word "Verwandlungswert" rather than "en- 
tropy" , although his colleagues suggested him to choose a name 
for his new thermodynamic quantity "S" (his chosen symbol for 
labeling "entropy", possibly in honor of "S"adi Carnot) that 
sounded as close as possible to the word "energy" . 



pressure and temperature, a reasonable assumption when 
dealing with the monster steam engines of the Victorian 
industry. A typical protein, however, is a few nanome- 
ters in size and consists of just a few tens of thousands 
of atoms. As a consequence the movements of a pro- 
tein engine are swamped by the fluctuations resulting 
from the Brownian motion of its surroundings, which 
causes the energy of any its part to fluctuate continu- 
ally in units of kT, with k denoting the Boltzmann con- 
stant and T the temperature. The effects of such energy 
fluctuations were brilliantly demonstrated by Yanagida's 
group (Nishiyama et al.\ 12002 2001), who observed ki- 



nesin molecules climbing the cytoskeletal track in a jud- 
dering motion made up of random hesitations, jumps and 
even backward steps. Similar results have been reported 
for a host of protein engines. The key question in mod- 
ern thermodynamics is therefore how far energy fluctua- 
tions drive micro- and nano-engines beyond the limits of 
macroscopic laws. 

A revealing experiment was performed recently by |Bus- 
tamante et al. (|2005[), who first stretched a single RNA 



molecule, by optically tugging at a tiny plastic bead at- 
tached to one end, and then released the bead to study 
the effect of random energy fluctuations on the molecule 
recovery process. By repeating many identical stretching 
cycles, these authors found that the molecule "relaxation 
path" was different every time. In fact, the bead was 
drawing useful energy from the thermal motion of the 
suspension fluid and transforming it into motion. How- 
ever, by averaging over increasingly longer bead trajecto- 
ries that is, approaching a macroscopic situation, |Busta- 



mante et al. (20051 were able to reconcile their findings 
with the second law. These results lead to the conclu- 
sion that the straightforward extension of the second law 
to microscopic systems was ungranted; individual small 
systems do evolve under inherently nonequilibrium con- 
ditions. 

However, a decade ago Jarzynski (1997) showed that 
the demarcation line between equilibrium and non- 
equilibrium processes is not always as clear cut as we were 
used to think. Imagining a microscopic single molecule 
process, Jarzynski evaluated not the simple average of the 
change of the (random) work of the underlying perturbed 
nanosystem, as it was pulled away from equilibrium ac- 
cording to an arbitrary protocol of forcing, but rather 
the average of the exponential of that tailored nonequi- 
librium work. Remarkably, such a quantity turned out 
to be the same as for an adiabatically slow version of 
the same process, and most remarkably, equals the ex- 
ponential of the system's equilibrium free energy change. 
This result, also ex perimentally demonstrated by |Busta" 



mante et al. (2005), came as much of a surprise because 
it meant that information about macroscopic equilibrium 
was somehow buried inside individual, randomly fluctu- 
ating microscopic systems far from equilibrium, see also 



( Crooks 


1999 Gallavotti and Cohen 


2007 Tal 


kner et al. 2008| 20071. 



1995 Jarzynski 



There is an additional limitation of 19th-century ther- 
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modynamics that is potentially even more significant in 
the design and operation of engines at the sub-micron 
scales. Kelvin's thermodynamics was based on the sim- 
plifying notion of an isolated system. The laws of macro- 
scopic thermodynamics therefore apply only to systems 
that are either separated from their environment or cou- 
pled to it under controlled operating conditions, that is, 
measured in terms of the state variables of the system it- 
self. However, at variance with a cylinder inside a steam 
engine, protein engines do not (cannot) work in isolation. 

Very much on the footsteps of the 19th-century scien- 
tists and engineers, modern experimenters have probed 
the proteins that play a crucial role in the cell individu- 
ally by feeding them with energy by injecting some chem- 
ical fuel "by hand" (e.g., ATP molecules) or exerting 
mechanical actions of some sort (Astumian 1997 Bus 



tamante et al. 



2005). In their natural setting, however, 



life engines are just parts of a closely interconnected func- 
tional web that keeps a cell alive. The great challenge of 
systems biology is therefore to put our understanding of 
isolated life engines back into the real world of the cell. 



be artificially manufactured by many orders of magni- 
tude. Nevertheless, inorganic nanodevices are attracting 
growing interest as a viable option due to their potential 
simplicity and robustness, without forgetting that inor- 
ganic nanodevices may provide additional experimental 
access to the molecular machinery itself. 

With this review the authors intend to pursue further 
the inorganic approach to nanodevices, based on three 
main assumptions: (1) That, in view of the most re- 
cent developments on nonequilibrium thermodynamics, 
the science of nanodevices, regardless of the fabrication 
technique, is inseparable from the thermodynamics of 
microscopic engines (Hanggi et al. 2005); (2) That the 
fabrication techniques on the nanoscales become more 
and more performing following the trend of the last two 
decades; (3) That a better understanding of the molecu- 
lar machinery can help devise and implement new trans- 
port and control mechanisms for biology inspired nan- 
odevices. In other words, we bet on a two-way cross- 
fertilization between the biomimetic and the inorganic 
approach. 



1. Artificial nanodevices 



2. Brownian motors 



Nanotechnology has been intricately linked with bio- 
logical systems since its inception. Fascinated by the 
complexity and smallness of the cell, Feynman ( 1960 ) 



challenged the scientific community to "make a thing 
very small which does what we want" . In his visionary 
response, Drexler ( 1992 1 proposed to focus on protein 



synthesis as a pathway for creating nanoscale devices. 
Both Feynman and Drexler 's propositions have been met 
with much skepticism as accurate manipulations at the 
nanoscale were deemed impossible. However, in view of 



the recent advances in systems biology ( Gross 1999 1 , cel- 



lular mechanisms are now being cited as the key proof 
of the nanotechnological viability of devices with atomic 
precision. In spite of their established complementarity, 
a fundamental difference between systems biology and 
nanotechnology is their ultimate goal. Systems biology 
aims to uncover the fundamental operation of the cell in 
an effort to predict the exact response to specific stimuli 
and genetic variations, whereas nanotechnology is chiefly 
concerned with useful design. 

Manufacturing nanodevices through positional assem- 
bly and self-assembly of biological components available 
at the cellular level is the goal of the so-called biomimetic 
approach - as opposed to the inorganic approach aimed 
at fabricating nanomechanical devices in hard, inorganic 
materials (e.g., using modern lithographic techniques, 
atomic force and scanning tunneling microscopy, etc). 
Nature has already proven that it is possible to engi- 
neer complex machines on the nanoscale; there is an ex- 
isting framework of working components manufactured 
by Nature than can be used as a guide to develop our 
own biology inspired nanodevices. It is also true that the 
molecular machinery still outperforms anything that can 



Nature provided microorganisms, with characteristic 
sizes of about 10~ 5 m, with a variety of self-propulsion 
mechanisms, all of which pertain to motion induced by 
cyclic shape changes. During one such cycle the config- 
uration changes follow an asymmetric sequence, where 
one half cycle does not simply retrace the other half, in 
order to circumvent the absence of hydrodynamic inertia 
at the microscales, i.e., for low Reynolds numbers (Pur 



cell 19771. A typical example are motile bacteria that 



"swim" in a suspension fluid by rotating or waving their 
flagella (Astumian 2007| Astumian and Hanggi 20021. 



As anticipated above, a further complication arises when 
the moving parts of a (sub)micron-engine have molec- 
ular dimensions of 10 _8 m or so. In that case, diffusion 
caused by Brownian motion competes with self-propelled 
motion. For example, a molecular motor mechanism be- 
comes superior if at room temperature, and in a medium 
with viscosity close to that of water, a bacterium needs 
more time to diffuse a body length than it does to swim 
the same distance. A passive diffusive mechanism oper- 
ating alone simply becomes inefficient. 

A solution common to most cell machinery is to have 
molecular motors operating on a track that constrains 
the motion to essentially one dimension along a periodic 
sequence of wells and barriers. The energy barriers signif- 
icantly suppress the diffusion, while thermal noise plays 
a constructive role by providing a mechanism, thermal 
activation (Hanggi et al. 19901, by which motors can 



escape over the barriers. The energy necessary for di- 
rected motion is supplied by asymmetrically raising and 
lowering the barriers and wells, either via an external 
time-dependent modulation (e.g., due to the coupling 
with other motors) or by energy input from a nonequi- 
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librium source such as a chemical reaction, like the ATP 
hydrolysis. Thus, in agreement with the reasoning un- 
derlying the analysis 2 of the Smoluchowski-Feynman 



stylized ratchet engine (Feynman et al. 1963 Smolu 



chowski 1912 1, under appropriate nonequilibrium condi- 



tions, structural anisotropy can sustain directed motion. 
Such a device clearly does not violate the second law of 
thermodynamics because the very presence of nonequi- 
librium renders inoperative those limiting (thermal equi- 
librium) restrictions. 

In the case of a bacterium, as for any ordinary heat en- 
gine, the relevant state variables, namely, its position and 
the phase of the flagellum stroke, always cycle through 
one and the same periodic time sequence; the two vari- 
ables are tightly coupled and almost synchronized. In 
clear contrast to this familiar scenario, the state vari- 
ables of molecular motors are often loosely coupled due 
to the prominent action of fluctuations, a salient feature 
nicely captured by Hanggi who originally coined the term 



Brow nian motors in the feature (Bartussek and Hanggi 
1995| ). 3 

Important hallmarks of any genuine Brownian motor 



are ( Astumian and Hanggi 2002| Hanggi et al. 2005 1: (i) 
The presence of some amount of (not necessarily thermal) 
noise. The intricate interplay among nonlinearity, noise- 
activated escape dynamics and non-equilibrium driving 
implies that, generally, not even the direction of trans- 
port is a priori predictable; (ii) Some sort of symmetry- 
breaking supplemented by temporal periodicity (typi- 
cally via an unbiased, non-equilibrium forcing), if a cycli- 
cally operating device is involved. Therefore, not every 
small ratchet device falls under the category of Brownian 
motors. This holds true especially if the governing trans- 
port principle is deterministic, like in mechanical ratchet 
devices of macro- or mesoscopic size. 

The following prescriptions should be observed when 
designing an efficient Brownian motor: (a) Spatial and 
temporal periodicity critically affect rectification; (b) All 
acting forces and gradients must vanish after averaging 
over space, time, and statistical ensembles; (c) Random 
forces (of thermal, non-thermal, or even deterministic 
origin) assume a prominent role; (d) Detailed balance 
symmetry, ruling thermal equilibrium dynamics, must be 
broken by operating the device away from thermal equi- 
librium; (e) A symmetry-breaking mechanism must ap- 
ply. There exist several possibilities to induce symmetry- 
breaking. First, the spatial inversion symmetry of the 



2 Note in this context also the insightful examination of Feynman's 
analysis by |Parrondo and Espanol| l |1996| . 

3 The notion ot "molecular motor" is reserved within this review 
for motors specifying biological, intracelluar transport. Like- 
wise, the notion of "Brownian ratchet" or "thermal ratchet" is 
reserved for the operating principle of protein translocation pro- 
cesses. The latter term seemingly has been introduced by |Simon| 
et al.\ ( |1992} to describe isothermal trapping of Brownian parti- 
cles to drive protein translocation, see also in Wang and Ostcr 
12002k 



periodic system itself may be broken intrinsically; that 
is, already in the absence of non-equilibrium perturba- 
tions. This is the most common situation and typically 
involves a type of periodic, asymmetric ratchet potential. 
A second option consists in the use of an unbiased driv- 
ing force (deterministic or stochastic, alike) possessing 
non- vanishing, higher order odd time-correlations. Yet 
a third possibility arises via collective effects in coupled, 
perfectly symmetric non-equilibrium systems, namely in 
the form of spontaneous symmetry breaking. Note that 
in the latter two cases we speak of a Brownian motor dy- 
namics even though a ratchet potential is not necessarily 
involved. 



The reasoning of using unbiased thermal fluctuations 
and/or unbiased nonequilibrium perturbations to drive 
directed motion of particles and alike has seen several 
rediscoveries since the visionary works by Marian von 
Smoluchowski ( 1912 1 and Richard P. |Feyn man et al. 



(1963). From a historic perspective, the theme of di- 



rected transport and Brownian motors was put into the 
limelight of the statistical and biological physics research 
with the appearance of several ground-breaking works, 
both theoretical and experimental, which app eared dur- 
ing the (1992-1994) period. Most not a bly, Ajdari and 
|Prost|([T992|;|Astumian and Bier]Ql994D;|Bartussek et al. 



|Prostl 
P9941 
|nasco 
( [1994 



T994] ; |Chauwin et a?.|(|1994[);|Doering et al | ( |1994[);|Mag 



( |1993D;|Millonas and Dykmanj ( |1994f ; ]Prbst et al 



(1994); Rousselet et al. (1994) helped ignite a tumul- 



tuous activity in this topical area which kept growing 
until the present days. The readers may deepen their 
historical insight by consulting earlier introductory re 
ports such as those published by |Astumian| dT99~7]) ; |Astu- 



mian and Hanggi (2002); Hang gi and Bartussek| ( 1996 



|Hanggi et al. ( 2005[); Jiilicher et al. (19971; Reimann 
( |2002[ ); |Reimann and Hanggi| ( |2002| ~"^ 



This review focuses on the recent advances in the sci- 
ence of non-biological, artificial Brownian motors. 



In contrast to those reviews and feature articles men- 
tioned above, which beautifully cover the rich variety of 
possible Brownian motor scenarios and working princi- 
ples, our focus in this review is on non-biological, artifi- 
cial, mostly, solid state based Brownian motors. In this 
spirit the authors have attempted to present a compre- 
hensive overview of the present status of this field, includ- 
ing newest theory developments, most compelling experi- 
mental demonstrations, and first successful technological 
applications. Some closely related topics, such as the 
engineering of synthetic molecular motors and nanoma- 
chincs based on chemical species, are only briefly dis- 
cussed herein because comprehensive, up-to-date reviews 
have been published recently by research groups very ac- 
tive in that area QBalzani et al] |2006| [Kay et al] |2007| 
Kottas et al. 1 120051). 
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II. SINGLE-PARTICLE TRANSPORT 

Signal rectification schemes in the absence of noise 
have been known since long ago, especially in the elec- 
trical engineering literature. However, rectification in 
a nanodevices cannot ignore fluctuations and Brownian 
motion, in particular. New experiments on both bio- 
logical and artificial devices showed how noise rectifica- 
tion can actually be utilized to effectively control particle 
transport on the small scales. By now, noise rectifica- 
tion has become one of the most promising techniques 
for powering micro- and nanodevices. 

In order to set the stage, in the next subsection we first 
consider the case of systems where rectification cannot 
occur. In the following subsections we then single out all 
ingredients that do make rectification possible. 

Let us consider a Brownian particle with mass m, co- 
ordinate x(t), and friction coefficient 7 in one dimension, 
subjected to an external static force F and thermal noise 
The corresponding stochastic dynamics is described 
by the inertial Langevin equation 



= -V'(x) -mjx + F + ^t), 



(1) 



where V(x) is a periodic potential with period L, namely 
V(x + L) = V(x), and ' indicates the differentiation with 
respect to x and x the differentiation with respect to 
time t. Thermal fluctuations are modeled by a stationary 
Gaussian noise of vanishing mean, = 0, satisfying 

the fluctuation-dissipation relation 



mm) = 2D p s(t), 



(2) 



where the momentum-diffusion strength reads D p — 
m/ykT, with k denoting the Boltzmann constant, and T 
is the temperature of an equilibrium heat bath. 

In extremely small systems, particle dynamics and fluc- 
tuations occurring in biological and liquid environments 
are often described well by the overdamped limit of Eq. 
Q , - for an illustrative discussion of this fact see the ac- 
count given by |PurceH ( |1977[ )-, in terms of a the massless 
Langevin equation, which is driven by a position-diffusion 
D x = KT/{rnrt) = D, i.e., 



x = -V(x) + F + £(t) , 
with a corresponding noise correlation 

mm) = Ms® ■ 



(3) 



(4) 



In this overdamped regime, the inertia term mx can 
be dropped altogether with respect to the friction term 
—rwyx (Smoluchowski approximation). Here has 
been scaled to unity for convenience, i.e., D = kT. 

In the absence of an external bias, i.e. F — 0, the equi- 
librium stochastic mechanism in Eq. (fTl) cannot sustain 
a non-zero stationary current, i.e., (x(t)) = 0, no mat- 
ter what V{x). This can be readily proven upon solving 
the corresponding Fokker-Planck equation with periodic 




FIG. 1 (Color online) (a) Tilted periodic potential of Eq. 
([5}, V(x) - xF, with F = 0.1; (b) Locked-to -running transi- 
tions. The thresholds F\ and F3 of the hysteretic loop (dashed 
curves) and the zero temperature step at F% (solid curve) 
are marked explicitly. Parameter values: m = Vb = 1, i.e. 
F 3 = 1, and 7 = 0.03. 



A. Symmetric substrates 

Let us consider first the case when the periodic sub- 
strate with the potential V{x) is symmetric under reflec- 
tion, that is V(x — xq) — V(—x + xq) for certain xq with 
xq < xq < L. The most studied example is the symmet- 
ric washboard potential ( Riskenl 



1984) 



V(x) = -V sm(2irx/L), 



(5) 



displayed in presence of a static tilt-force F in Fig. [TJl. 
The particle mobility fi(F) = (x) / F is symmetric for 
F — > — F, namely, n(F) = fi(—F). For this reason in 
this section we restrict ourselves to F > 0. 



1. dc drive 

The Brownian motion, Eq . (|lj), in the washboard po- 
tential, Eq. ([5J, is detailed in Risken (19841 textbook. To 
make contact with Risken's notation one must rewrite 
Eq. (nj) in terms of the rescaled quantities x lz 



F 



and (77)^, so that 



x, 



-7 x + uiq cos x + F + £(t). 



(6) 



boundary conditions ( jMel'nikov 1991 Risken 19841 



The angular frequency wo = (2n/L) \JVq /m characterizes 
the oscillating motion of the particle at the bottom of 
the potential wells. To reconcile Eq. ([6| with Eq. 
for the potential in Eq. ([5]), it suffices to scale m = 1 
and L = 2tt, as assumed throughout this section, unless 
stated otherwise. 

The dynamics of Eq. ^ is characterized by ran- 
dom switches occurring either between two locked states, 
dwelling a well minimum, or switches between a locked 
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state with zero average velocity and a (downhill) run- 
ning state with a finite, average asymptotic velocity 
(x) = F/j. In terms of mobility, locked and running 
states correspond to jfi = and 7/x = 1, respectively. 

In the overdamped regime, Eq. the particle is 

trapped most of the time at a local minimum of the tilted 
substrate as long as F < F3 = ujq ; for F > F3 there exist 
no such minima and the particle runs in the F direction 
with average speed approaching F/j. This behavior is 



described quantitatively by the mobility formula (Risken 
[19841 ): 



(t(L,F))F> 



(7) 



where (t(L,F)), the mean first-passage time of the par- 
ticle across a substrate unit cell in the F direction, can 



be computed explicitly for any choice of V(x) (Hanggi 



et al. 1990 1 



In the undcrdamped regime 7 « wo the locked-to- 
running transition depends crucially on the presence of 
noise, no matter how weak. In the noiseless case, £(t) = 
0, the average speed of a Brownian particle with coordi- 
nate x(t) depends on the initial conditions according to 
the hysteretic cycle illustrated in Fig. [TJd: The transition 
from the locked to the running state occurs on raising F 
above F3 (depinning threshold) , while the opposite tran- 
sition takes place on lowering F below F\ — (4/7r)7a;o (re- 
pinning threshold). Of course, for sufficiently large values 
of 7, say, in the damped regime with 7 > 70 = (7r/4)u;o, 
the distinction between Fi and F$ becomes meaningless; 
the locked-to-running transition occurs at F — F 3l no 
matter what the initial conditions. 

At zero temperature, T = 0+, we have a totally dif- 
ferent scenario, as the stationary dynamics of x{t) is 
controlled by one threshold F2 — 3.3576 ... 7W0, only 
(Risken 19841: For F < F2 the Brownian particle re- 



mains trapped in one potential well; for F > F2 it falls 
down the tilted washboard potential with speed close to 
F/j. At the threshold F2 the quantity 7/i jumps from 
zero to (very close to) one, stepwise (Fig. [lj>). Note that, 
at variance with F3, F2 indicates a dynamical transition 
occurring in the presence of a relatively small tilt. Nev- 
ertheless, at low damping, switches between locked and 
running states correspond to long forward particle jumps, 
which can span over many substrate unit cells L; the dis- 
tributions of the relevant jump lengths exhibit persistent 



non-exponential tails ( 


Borromeo and Marchesoni 


2000 


Costantini and Marchesoni| 
Pollak et al.\ |1993| |Shushin 


1999 


|Ferrando et al. 


1995 


2002 


1. 



For finite but low temperatures, kT <C Vq, the transi- 
tion from the locked to the running state is continuous 
but still confined within a narrow neighborhood of F?; 
the relevant locked-to-running transition threshold Fth 
is numerically identifiable with high accuracy as a con- 
vex function of 7; F2 and F3 are respectively the 7 — > 
and 7 — > 00 asymptotes of the numerical curve -Fth (7) 



2. ac drive 

Suppose now that the external force F(t) acting on 
the unit-mass Brownian particle is periodic in time. The 
simplest case possible is represented by a noiseless par- 
ticle, £(t) = 0, moving on a sinusoidal potential, Eq. (J5j, 
under the action of a harmonic force 



F(t) = Aicos(fiit + 0i). 



(8) 



The wide class of devices thus modeled may be assim- 
ilated to a damped-driven pendulum operated at zero 
noise level, a chaotic system investigated at depth in 
the 1980s (iBaker and Gollubl 119901). The dynamics of 



a massive particle in a sinusoidal potential was repro- 



duced in terms of a climbing-sine map ( 


Geisel and Nier- 


wetberg 


1982 


Grossmann and Fujisaka 


19821: Run- 



ning orbits, leading to a spontaneous symmetry break- 
ing, can be either periodic or diffusive, depending on 
the value of the map control parameter (viz. the am- 
plitude of the sine term). The phase-space portrait of 
the actual damped-driven pendulum was computed by 



Huberman et al. 



(Risken 1984) 



who revealed the existence of 
delocalized strange attractors with an intricate structure 
on all scales , later recognized to be fractal objects ( |Baker 
and Gollub 19901. This means that, despite the global 
reflection symmetry of the dynamics in Eqs. (JT|) and pi, 
for sufficiently small 7 the particle drifts either to the 
right or the to left, depending on the initial conditions, 
but with equal probability in phase-space. 

Coupling the particle to a heat bath, no matter how 
low the temperature, changes this scenario completely 
The action of the noise source amounts to scram- 
bling the initial conditions, which therefore must be av- 
eraged on. As a consequence, trajectories to the right 
and to the left compensate one another and the symme- 
try is restored: No Brownian drift is expected in the zero 
temperature limit. 

A system symmetry can be broken by adding a dc com- 
ponent Aq to the external force, 

F{t)=A + A 1 cx)s(n 1 t + (/>!). (9) 

The most evident effect of such a periodic drive is the 
appearance of hysteresis loops ( Borromeo et al. 1999[ ) in 
the parametric curves of the mobility n(t) versus F(t). 
For Oi <C an d Ao -C F t h, the mobility hysteresis 
loop is centered on the static mobility curve h(Aq) and 
traversed counterclockwise; with decreasing fli, its major 
axis approaches the tangent to the curve n(A ). Hystere- 
sis loops have been observed even for f2i much smaller 
than the relevant Kramers rate, the smallest relaxation 
rate in the unperturbed stationary process of Eq. ^ 
with F = 0. The area encircled by the hysteretic loops 
is maximum for Aq ~ F t h, that is close to the transition 
jump. Even more interestingly, it exhibits a resonant 
dependence on both the forcing frequency and the tem- 
perature, thus pointing to a dynamical stochastic res- 
onance (Gammaitoni et al. 19981 mechanism between 



7 



locked and running states (Borromeo and Marchesoni 
2000b. 



Finally, Machura et al. ( 2007a I showed that under spe 



cial conditions, involving small Brownian motion at tem- 
perature T and appropriate relaxation constants 7 < 
^1 J$ the damped process in Eq. (j6j) occasionally 
exhibits the phenomenon of absolute negative mobility: 
The ac cycle averaged drift velocity (x) may be oriented 
against the dc bias A , as the result of a delicate inter- 
play of chaotic and stochastic dynamics. This observa- 



tion was corrobora t ed the oretically by Speer et al. (|2007 1 
and Kostur et al. ( 2008 1 and experimentally by INagel 



et al. (20081. Yet, another artificial Brownian motor sys- 
tem where this phenomenon can likely be validated ex- 
perimentally are cold atoms dwelling in periodic optical 
lattices. 



3. Diffusion peak 

Diffusive transport of particles or, more generally, 
small objects is a ubiquitous feature of physical and 
chemical reaction systems ( Burada et al. 2009 ) . Directed 



Brownian motor transport is typically controlled both 
by the fluctuation statistics of the jittering objects and 
the phase space available to their dynamics. For exam- 
ple, as the particle in Eq. (JsJ) drifts with average speed 
(x) in the direction of the external force F, the random 
switches between locked and running state cause a spa- 
tial dispersion of the particle around its average position 
The corresponding normal diffusion constant, 

Mi) 2 } - {At)) 2 



D(F) = lim ■ 

t — >-OC 



2t 



(10) 



was computed numerically as a function of F at constant 
temperature by Costantini and Marchesoni ( |1999 1. A 
peak in the curves of D versus F is detectable in the 
vicinity of the transition threshold F t h, irrespective of 
the value of 7. In particular, at low damping, where 
Fth ~ Fa, and low temperature, kT <C u)q, the peak of D 
at F — F2 is very pronounced; on increasing the damp- 
ing, the diffusion peak eventually shrinks down to just a 
bump corresponding to the threshold Fth — F3. In any 
case, in an appropriate neighborhood of F t h the diffu- 
sion constant can grow larger than Einstein's diffusion 
constant for the free Brownian motion in one dimension, 
Do = fcT/7. On increasing T the diffusion peak eventu- 
ally disappears, no matter what 7. 

A refined analytical formula for the diffusion peak was 
obtained by Reimann and collaborators who regarded the 
locked-to-running transition in the overdamped regime as 
a renewal process, that is ( Reimann et al.\ 2001 2002): 



D(F) 



(11) 



L 2 (t*(L,F))-(t(L,F))* 
2 (t{L,F))* 

where the n-th moments of the first passage time, 
{t n (L, F)), can be computed explicitly. For F — ► 0, Eq. 
( 11 1 reproduces the zero-bias identity D/D = jfi (Festa 
and d'Agliano||1978[). 



4. Single-file diffusion 

When a gas of particles is confined to a linear path, an 
individual particle cannot diffuse past its neighbors. This 
constrained ID geometry is often called a "single-file" , 



Harris 


1974 


Jepsen 


1965 



particles moving on a tilted sinusoidal substrate in Eq. 
(l6| of length L. If the particle interaction is hard-core 
(zero radius) , the file constituents can be labeled accord- 
ing to a given ordered sequence and the long-time dif- 
fusion of an individual particle gets strongly suppressed. 



In early studies ( jHarris 1974 
and Percusl [IMniLevittl 11973 



Jepsen 1965 


Lebowitz 


Marchesoni and Taloni| 



the thermodynamic limit (L, N — > 00 with constant den- 
sity p — N/L) was calculated analytically. Those results 
were generalized to the diffusion of a single-file of driven 
Brownian particles on a periodic substrate by |Taloni and| 



Marchesoni ( 2006 1 who derived the subdiffusive law 



lim [(x(tf) - Mi)) 2 ] = - 

t— >oo n 



D(F)t 



(12) 



with D(F) given in Eq. (111. This result applies to any 



choice of the substrate potential V{x) (Burada et al. 



2009 1 and to the transport of composite objects ( Hein- 



salu et al. 20081, as well. Excess diffusion peaks have 
been obtained experimentally in the context of particle 
transport in quasi-lD systems (Sec. 



IIII 



B. Rectification of asymmetric processes 

Stochastic transport across a device surely can be in- 
duced by applying a macroscopic gradient, like a force 
or a temperature difference. However, under many prac- 
tical circumstances this is no viable option: (a) Current 
induced by macroscopic gradients are rarely selective; (b) 
A target particle that carries no charge or dipole, can 
hardly be manipulated by means of an external field of 
force; (c) External controls, including powering, critically 
overburden the design of a small device. Ideally, the op- 
timal solution would be a self-propelled device that op- 
erates by rectifying environmental signals. In the quest 
for rectification mechanisms of easy implementation, we 
will start from the symmetric dynamics of Eq. ([lj and 
add the minimal ingredients needed to induce and control 
particle transport. 

If V{x) is symmetric under reflection, the only way to 
induce a drift of the Brownian particle consists in driv- 



deterministic or random ( 


Astumian and Hanggi 


2002 


Chialvo et al. \ |1997| |Hiinggi et al. 1996 Luczka et al. 


1995 


Reimann and Hanggi 


20021. Here, "symmetric" 



means that all force moments are invariant under sign 
reversal F — > — F. Note that the condition of a vanish- 
ing dc component, lim^oo j J Q F(s)ds = 0, would not be 
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sufficient. For instance, a bi-harmonic signal with com- 
mensurate frequencies and arbitrary phase constants, al- 
though zero-mean valued, is in general nonsymmetric. 

On the contrary, particles in an asymmetric potential 
can drift on average in one direction even when all per- 
turbing forces or gradients are symmetric. However, as 
pointed out in Sec. |II.D[ to achieve directed transport in 
such a class of devices, the external perturbation is re- 
quired to be at least time-correlated, like in the presence 
of a non-Markovian noise source (correlation ratchets) or 
a time periodic drive (rocked and pulsated ratchets) . 

The interplay of time and space asymmetry in the 
operation of a Brownian motor has been established 
on firm er mat hematical groun ds by |Yevtushenko et ~aL\ 
(2001) and by Reimann (2001). Let us slightly general- 
ize the overdamped dynamics in Eq. ([3| to incorporate 
the case of time dependent substrates, that is 



x = -V'[x,F 2 (t)]+F 1 (t) + Z{t). 



(13) 



The potential V\ x, F 2 (t )] is termed sup ersymmetric 
(Jung and Hanggi 1991| Marchesoni et al. 1988) if, for 
an appropriate choice of the x and t origins, 



V[x,F 2 (t)] = V[x + L/2 ) F 2 {-t)]. 



(14) 



Analogously, the additive drive Fi(t) is supersymmetric 
if for an appropriate t origin, 



-F 1 (t) = F 1 (-t). 



(15) 



Should Fi(t), with i = 1,2, be stationary noises, clearly 
no restriction can be set on the t origin; the equalities in 
Eq. " 



|14| and in Eq. ( 15 1 must then hold in a statistical 
sense, meaning that the two terms of each equality must 
be statistically indistinguishable. 

Let us consider now the time reversed process z(t) = 
x(—t) + L/2. By definition, (i) = —(&), whereas, 
on simultaneously imposing the supersymmetry condi- 
tions p"4| and (15 1, the Langevin equation (13 1 yields 
(i) = (x); hence7(i) — 0- As a consequence, a nonzero 
rectification current requires that either the substrate or 
the additive drive (or both) are non-super •symmetric. 

We remark that the above theorem has been proven 
only for zero-mass particles, that is, when the Smolu- 
chowski regime mx = applies. In the presence of in- 
ertia, instead, rectification may occur, under very spe- 
cial conditions, also in fully supersymmetric devices, as 
shown lately by Machura et al.\ ( 2007a I for a rocked co- 
sine potential. 



C. Nonlinear mechanisms 

In this Section we review transport on symmetric sub- 
strates driven by asymmetric forces. The rectification 
mechanisms outlined below can be traced back to the 
nonlinear nature of the substrate; for this reason, at vari- 
ance with Brownian motors, they work also, if not more 
effectively, in the absence of noise. 



We remind the reader that these mechanisms have 
been introduced and demonstrated experimentally in the 
most diverse fields of physics and engineering. Direct ap- 
plications to various categories of artificial nanodevices 
will be discussed in the subsequent sections. 



1. Harmonic mixing 

A charged particle spatially confined by a nonlinear 
force is capable of mixing two alternating input electric 
fields of angular frequencies f2i and r2 2 , its response con- 
taining all possible higher harmonics of f2i and f^- For 
commensurate input frequencies, i.e., / O2 = n/m with 
n and m coprime integer numbers, the output contains 



a dc component, too (Schneider and Seeger 



19661; har- 



monic mixing thus induces a rectification effect of the 
(n + m)-th order in the dynamical parameters of the sys- 
tem (Goychuk and Hanggi 1998 Marchesoni 19861. 



Let us consider the overdamped stochastic dynamics 
of Eq. ([3]) in the potential of Eq. ([5]), driven by the 
bi-harmonic force 

F(t) = A x cos(«ii + 0i) + A 2 cos(fi 2 i + <fe) (16) 

Let the two harmonic components of F(t) be commensu- 
rate with one another, meaning that fii, fl 2 are integer- 
valued multiples of a fundamental frequency fi , i.e., 
Oi = nflo and Cl 2 — to£!o- For small amplitudes and 
low frequencies of the drive in Eq. (16), a simple ex- 
pansion of the mobility function, Eq. (|7j), in powers of 
F(t), yields, after time averaging, the following approx- 
imate expression for the non-vanishing dc component of 
the particle velocity: 



2/1 



n f A x 



1! r> \ 



A2 

2 



cos(A mi „), (17) 



where JL is the positive (n + m — l)-th derivative of fJ.(F) 
at F = Q, A m! „ = nifi 2 — m<j>i, and m + n is an odd 
number. Harmonic mixing currents j — (x)/L clearly re- 
sult from a spontaneous symmetry breaking mechanism: 
indeed, averaging j over (f)\ or (j) 2 would eliminate the 
effect completely. This is true under any regime of tem- 
perature and forcing, as proven by means of standard 
perturbation techniques ( |Breymayer and Wonneberger 
" 1986| 



1981 Marchesoni 



In particular, rectification of 
two small commensurate driving signals, Ax, A 2 -C ^q, is 
a noise assisted process and, therefore, is strongly sup- 
press ed for kT <C ujq, when Jl decays exponentially to 
zero (Risken, 1984). Moreover, we anticipate that ac- 



counting for finite inertia effects requires introducing in 
Eq. ( 17 1 an additional damping and frequency dependent 
phase-lag, as discussed in Sec. |IV.A| 

In all calculations of (x), including Eq. (17 1, the re- 
flection symmetry of V(x) plays no significant role; recti- 
fication via harmonic mixing is caused solely by the non- 
linearity of the substrate. However, it must be noticed, 
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2. Gating 

A periodically-driven Brownian motion can also be rec- 
tified by modulating the amplitude of the substrate po- 
tential V(x). Let us specialize the overdamped dynamics 
in Eq. ( p~3[ ) as follows: 

(19) 



x = -V'{x)[l + F 2 {t)] + F 1 {t)+i{t) 



0.5 



D 



1.5 



2.5 



FIG. 2 Gating mechanism. Equation ( 19 1 has been simulated 
numerically for fii = f?2 = 0.01, <j>i = 4>2 = 0, A\ = A2 = A, 
and V(x) = Wq(1 — cos a;) with wo = 1; the net velocity (x) has 
been plotted versus D = kT for different amplitudes A. For a 
comparison we show (x) versus D for fii = 0.01, n 2 /ni = V2 
and A = 0.5 (crosses). Inset: (x) versus <f>2 for (j>i = 0, 
A = 0.5, and D = 1. Note the resonant behavior of (i) for 
subthreshold drives, A < 1, where a low noise level enhances 
rectification; increasing D for A > 1 degrades the rectification 
effect. 



that a symmetric nonlinear device cannot mix rectangu- 
lar waveforms, like 



F(t) = A 1 sgn[cos(n 1 t+0 1 )]+A2Sgn[cos(n 2 t+0 2 )], (18) 



with Ai,A 2 > and sgn[. . . ] denoting the sign of [. . . ] 



To avoid interference with possible harmonic mixing ef- 
fects we follow the prescription of Sec. |II.C.1| namely we 
take Fi(t) = AjSgn[cos(fM + fa)], i = 1,2, with Ai > 
and < A 2 < 1 . Mixing of the additive F x (t) and the 
multiplicative signal F 2 (t) provides a control mechanism 
of potential interest in device design. 

In the adiabatic limit, the Brownian particle moves 
cyclicly back and forth subjected to opposite dc forces 
with amplitude A\\ the substrate potential V(x,t) = 
V(x)[l + F 2 (t)] switches, in turn, periodically between 
the two symmetric configurations V±(x) = V(x)(l±A 2 ). 
The relevant mobilities fi±(Ai) can be easily related to 
the static mobility function n(A) for the tilted potential 
V(x) studied in Fig , [I] namely ( Savel'ev et al. 2004a 



Savel'ev et cd] |2004[ ) 



li ± {A 1 ) = (l±A 2 )n 



l±Ao 



(20) 



withT->T/(l±A 2 ). 

In particular, for any pair of commensurate frequencies 
ill and il 2 such that f^/^i = (2m — l)/(2n — 1) (with 
m, n positive integers), the net particle velocity mediated 
over an integer number of cycles of both Fi(t) and F 2 (t) 
can be cast in the form 



(-l) m+n Ai 
(2m- l)(2n- 1) 



A/j,(Ai,A 2 ) p(A 

n.m ) 1 



USayel'ev et |2004a| |Saverev et al\\M)^\ . As shown wnere A/j,(A 1 ,A 2 ) = fX-(Ai) - ^+{A X ) and p(A r 



in Sec. II.D.4| asymmetric devices do not exhibit this pe 



(21) 

n) = 



culiarity. Moreover, we underscore that for incommensu- 
rate frequencies, £l±/Q 2 not a rational number, harmonic 
mixing takes place only in the presence of spatial asym- 
metry, as reported in Sec. |II.D.4| 

For practical applications, any commensuration con- 
dition, like fli/fl 2 — n/m for harmonic mixing (but see 
also Sees. II. C. 2 and II.D.4[ ) is affected by an uncertainty 
that is inverse proportional to the observation time, i.e., 
the code running time for numerical simulations or the 
data recording time for real experiments. As such time 
is necessarily finite, only a limited number of commensu- 
ration spikes can be actually detected. 

We conclude by reminding that the notion of harmonic 
mixing has been introduced long ago, e.g., to interpret 
the output of certain charge wave density experiments 
(Schneider and Seeger 19661 and to design laser ring 



A n ^ rn \/'K — 0.5 is a modulation factor with A T , 



gyroscopes (Chow et al. 1985) and annular Josephson (Borromco and Marchesoni 



junctions (|Ustinov et al. 20041; applications in the con 



(2n—\)4>2 — (2m— 1)4>\, mod(27r). Note that for different 
choices of the ratio O2/O1, no induced drift is predicted 
(see Fig. [2}. 

As a consequence, a relatively small modulation of the 
sinusoidal potential amplitude at low temperatures re- 
sults in a net transport current as illustrated in Fig. [2] 
Let us consider the simplest possible case, fii = and 
4>i = < / ) 2 : As the ac drive points to the right, the ampli- 
tude of V(x, t) is set at its maximum value ojq(1 + A 2 ); at 
low temperatures the Brownian particle can hardly over- 
come the substrate barriers within a half ac-drive period 
n/fli. In the subsequent half period i*i(t) switches to the 
left, while the potential amplitude drops to its minimum 
value U!q(1 — A 2 ): The particle has a better chance to 
escape a potential well to the left than to the right, thus 
inducing a negative current with maximum intensity for 
Oi of the order of the Kramers rate (resonant activation 

2004[ )). Of course, ampli- 



text of nano-particle transport is more recent (see Sees. 
II.C.3|and|rvT). 



tude and sign of the net current may be controlled via 
the modulation parameters A 2 and <p 2 , too (see inset of 
Fig.|}. 
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3. Noise induced transport 

Induced transport in the symmetric dynamics can be 
achieved also by employing two correlated noisy signals. 

(i) Noise mixing. |Borromeo and Marchesoni| ( |2004[ ) 
showed that the gating process in Eq. (|19[) can be driven 
also by two stationary, zero-mean valued Gaussian noises, 
Fi(t) = rji(t) with i = 1,2. The two random drives may 
be cross- and auto-correlated with 



A, 



\t\ 



(J7i(t)»ft(0)) = ^/Q~Q] — exp I - — I , i,j = 1,2. 

Tij \ TV,' / 

(22) 

Without loss of generality, one sets An = A 2 2 = 1 and 
A12 = A21 = A, and to avoid technical complications, 
= t. Of course t — > corresponds to taking the 
white noise limit of T]i(t), Eq. @. The parameter A 
characterizes the cross-correlation of the two signals; in 
particular, A = 0: independent signals; A = 1: identical 
signals, r]2 (f) = 771(f); A = — 1: signals with reversed 
sign, 772(f) = — /71 (t) . Two such signals may have been 
generated by a unique noise source and then partially 
de-correlated through different transmission or coupling 
mechanisms. 

In the white noise limit, t — ► 0, the Fokker-Planck 
equation associated with the process of Eqs. (19 1 and 



( 22 1 admits of a stationary solution in closed form. Non- 
vanishing values of (x) for A 7^ are the signature of 
a stochastic symmetry breaking due to nonlinear noise 
mixing; the sign of A, similarly to the relative phase A m)Tl 
in the gating set-up of Sec. |II.C.2| determines the direc- 
tion of the particle drift. The interpretation of this result 
is straightforward. Let us consider, for instance, the case 
of the symmetric potential in Eq. (J5J), rocked and pul- 
sated by the same signal, i.e., 771(f) = 772(f): For A = 1, 
when pushed to the left (rji < 0), the Brownian particle 
encounters lower substrate barriers than when pushed 
to the right (rji > 0) , hence the negative average drift 
current detected by means of numerical simula tion (|Bor-| 
romeo ef al.\ 2006| |Borromeo and MarchesonTj |2005b| . 

The magnitude of the i nduced current is controlled 
by the correlation time r (Borromeo and Marchesoni 
2005b I . While one can ea sily estimate (x) for A ^ 
(white noise (Risken 1984)) or r 



and r 

(strongly correlated noise (Hanggi et al 



1989)), 



00 
the 

intermediate r values are accessible solely through 
numerical simulation. In Fig. [3][a) we report (i) versus 
t for different noise intensities: the two sets of curves at 
Q and Q/t fixed, illustrate well the noise mixing effect 
for t — > and r — > 00, respectively. 

(ii) Noise recycling. If the noises r/i(t) are generated by 
the same source and then coupled to the diffusing particle 
through different paths, it may well happen that they are 
simply delayed in time. Under stationary conditions, we 
can assume that 772(f) = 771 (f — Td), with 77.; (f) given in 
Eq. p2]) and, for simplicity, A = 0. In the notation of 
Refs. (Borrom eo ef al.\ 2006 Borromeo and Marchesoni 



Q/T = 0.25 
Q/T = 0.5 
■ « Q = 0.5 

0--0 Q=I.O 




FIG. 3 (Color online) Equation |l9| , with Fi(t) = rj,(t), 
i — 1,2, and V(x) — — sin a;, has been simulated numerically; 
the n oises rji(t), Eq. ( 22 l, have same str ength, Qi — Q2 ~ Q- 
After Borromeo and Marchesoni (2005b): (a) Nonlinear noise 
mixing] Characteristics curve {x)-t for different intensities 
Q of the noises and A = 1; D = 0.1. (b) Noise recycling: 
Characteristics curve v-Ta, where Td is the relative time delay 
of r/i (see text). Data are for different r and Q = 1; D = 0, 
see in Borromeo et al. (20061. 



2007b), 771(f) represents the primary noise source and 



772(f) a recycled noise to be used as a control signal. 

By the same argument as for noise mixing, we expect 
that the Brownian dynamics gets rectified to the left with 
negative velocity (x(Td)}. Note that: (£(— r^)) = (x(rd)) 
and (x(rd)} — > — (x(Td)), upon changing the relative signs 
of ?7j. The dependence of the characteristic curve (i^)) 
on the time constants Td and r displayed in Fig. [3^b), 
is important in view of practical applications. Indeed, 
in many circumstances, it would be extremely difficult 
to recycle a control signal 772(f) so that <§; r; stated 
otherwise, measuring (i(0)) requires a certain degree of 
experimental sophistication. On the contrary, if we agree 
to work on the resonant tail of its response curve | (x(rd)) \ , 
a noise controlled rectification device can be operated 
with less effort; its net output current may be not the 
highest for Td > t, but is still appreciable and, more 
importantly, stable against the accidental floating of the 
time constant Td- 

In this sense both schemes discussed in this section are 
a simple-minded attempt at implementing the operation 
of a Maxwell's daemon: The ideal device we set up is in- 
tended to gauge the primary random signal 771 at the sam- 
pling time f and, depending on the sign of each reading, 
to lower or raise the gate barriers accordingly at a later 
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time t + Td, i.e., open or close the trap door. The rectify- 
ing power of such a daemon is far from optimal; lacking 
the dexterity of Maxwell's "gate-keeper" (ILeff and Rex 



like an automaton. 



D. Brownian motors 



2003 Maruyama et al. 20081, it only works "in average 



As detailed in Sec. |II.B[ a necessary condition for the 
rectification of symmetric signals, random or periodic in 
time, alike, is the spatial asymmetry of the substrate. 
Rectification devices involving asymmetric substrates are 
termed ratchets. In most such devices, however, noise 
(no matter what its source, i.e. stochastic, or chaotic, 
or thermal) plays a non-negligible, or even the dominant 
role. Under such conditions one speaks of Brownian mo- 



tors i 


Astumian 


1997 Astumian and Hanggi| 2002 


Bar- 


tusse 


t and Hanggi 


1995 


Hanggi and Bartussck 


1996 


Hanggi et al.\ 2005 


2002 


|Kay et al. 2007| Linke 


2002 


Reimann et al. 


199C 


> |. The label Brownian motor should 



II.D.2 



not be abused to refer to all small ratchet-like devices. 
For instance, the rocked ratchets of Sec. |II.D.1| work 
quite differently in the presence or in the absence of noise, 
whereas the pulsated and correlation ratchets (sometimes 
also referred to as "stochastic" ratchets) of Sees 
and |II.PT3| work only in the presence of noise. 

The hallmarks of genuine Brownian motors arc listed in 
Sec. |U In this Section we discuss in detail noise rectifica- 
tion and directed transport on asymmetric periodic sub- 
strates and potentials. We caution the reader that: (1) 
Strict periodicity is not a requirement for the operation 
of a Brownian motor. The ratchet system may contain 



small amounts of disorder ( 


Harms and Lipowsky 


1997 


Kafri and Nelsonj |2005| |Martinez and Chacon[ 


2008 


Popescu et al. 


2000 


) or even be non-periodic ( 


March- 



csom , 



1997) 



(2) Spatial asymmetry can also result as 
a collective effect, for instance, in the extended systems 
consisting of interacting, symmetric dynamical compo- 
nents, introduced in Sec. [V] 

The archetypal model of ratchet substrates in ID is the 



double-sine potential proposed by Bartussek et al. ( 1994 1 



V(x) 



/2ttx\ 


+ ^sin 

















(23) 



or, in the rescaled units of Eq. J 6]), V(x) = — w^sina; + 
\ sin(2x)]. In the sketch of Fig. B the barriers are skewed 
to the right with ratchet length > /_ . A Brownian par- 
ticle with Langevin equation ([I]) moving on such a sub- 
strate is characterized by an asymmetric mobility func- 
tion, fi(F) ^ fi(—F); as the particle mobility depends 
now not only on the amplitude, but also on the orienta- 
tion of the drive F, symmetric substrate perturbations 
are expected to induce a net current in either direction. 

Broadly speaking, ratchet devices fall into three cate- 
gories depending on how the applied perturbation cou- 
ples to the substrate asymmetry. 




FIG. 4 (Color online) Ratchet mechanism, (a) Sketch of the 
potential in Eq. (23| with L = 1 and Vo = L/2tt. The 
three consecutive extremal points, X- — —0.19, xq — 0.19 
and x+ = 0.81, define a potential well; the barrier height 
is AV = V(x±) — V(xo) — 0.35 and its asymmetry is 
quantified by the difference — I- with l± = \x± — xo\. 
The curvature of the potential at the bottom of the well is 
ujI = (3^/2)^0 - 10.1. (b) Rocket ratchet. Probabil- 
ity density p(x, t) for a Brownian particle initially centered 
around x = (red) and then driven for t — 10 by a dc force 
Aq — —0.5 (green) and 0.5 (grey). The backward displace- 
ment is strongly suppressed, hence the positive natural orien- 
tation of this ratchet, (c) Pulsated ratchet. The ratchet po- 
tential is switched "on" and "off" periodically with t° n — 1, 
t oB = 3 and T 2 = t on + t oB (see text). The particle, ini- 
tially set at x = xo, relaxes first in the starting well during 
t° n (curve 1, red, ti — t° n ), then diffuses symmetrically only 
driven by noise for t° (curve 2, blue, ti = T%), and finally 
gets retrapped in the neighboring wells as the next cycle be- 
gins (curve 3, green, t% = T2 + t° n ). As the left side peak of 
p{x,tz) is more pronounced than the right peak, the natural 
orientation of this ratchet is negative. In the simulations of 
(b) and (c) the ratchet potential V(x) is the same as in (a) 
(also drawn to guide the eye) and the intensity of the noise 
in Eq. (El is D — 0.1 AV. Courtesy of Marcello Borromeo. 



1. Rocked ratchets 

Let us consider first, for simplicity, the Langevin equa- 
tion ([3]) with the sinusoidal drive of Eq. ^ . When ap- 
plied to the reflection-symmetric sine potential of Eq. 
([5]), F(t) breaks instantaneously the symmetry of the 
potential by tilting it to the right for F(t) > 0, and 
to the left for F(t) < 0. However, due to the spatio- 
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temporal symmetry of the process, V(x) — V(—x) and 
F(t) = F(—t), over one drive cycle, T\ — 2ir/Cli, the 
drift to the right and to the left compensate one another: 
the net current is null. When placed in the double-sine 
potential of Eq. ( 23 1 , instead, an overdamped particle is 
apparently more movable to the right than to the left. In- 
deed, depinning occurs at F — F 3 r = 3/4 and F = —F^ 
for F(t) oriented respectively to the right and to the left; 
moreover from the asymmetry condition l + > Z_ folows 
immediately the inequality F^r < F^l. The depinning 
thresholds of the potential of Eq. ( 23 1 in the rescaled 
units of Eq. (|6) are F 3R = 3/4 and F 3L = 3/2. 

If the forcing frequency Sli is taken much smaller than 
all the intrawell relaxation constants, ui^Ti ^> 1 and 
T\ (t(L,0)), than the net transport current can be 
computed in the adiabatic approximation by averaging 
the instantaneous velocity x(t) = /j,\F(i)]F(t), se e Eq. 



(|7[), over one forcing cy cle, that is (Bartussek et al. 1994 



Borromeo et al. 



2002| ) 



J 



[X) 

L 



1 



sgn[F(s)] 
(t[L,F(s)}) 



ds. 



(24) 



The adiabatic current in Eq. p4[ ) is positive definite, 
as the natural ratchet direction is defined by the choice 
Fsr < Fsl- However, its range of validity is restricted 
to extremely low temperatures T and drive frequencies 
ill. On raising fli, the rectification current develops a 
more complicated dependence on the noise intensity D 
and the drive amplitude A±, as proven by the numerical 
simulations reported in panels (a)-(c) of Fig. [5] This is 
due to the fact that for forcing periods T\ shorter than 
the mean first-passage time (t(L, 0)), the driven parti- 
cle oscillations span over fewer substrate cells. A few 
properties of the rocked ratchet current in Fig. [5] are re- 
markable: 

(a) At finite temperature, the maximum rectification ef- 
fect occurs right in the adiabatic regime, Oi — > 0, and 
with natural orientation, panels (a) and (b); 

(b) For Oi > ujq, the drive and the interwell oscilla- 
tions may combine to reverse the current orientation. 
Such current inversions are restricted to select D and 
Ai ranges, but never in the absence of noise, panel (a) 
and (b); 

(c) In the noiseless regime, T = 0, the incommensuration 
between forced oscillation amplitude and substrate peri- 
odicity causes the quantized locking structure of panel 

(c) . Such a structure disappears either in the limit 
17 1 — > 0, as the forced oscillation amplitude grows in- 
finitely large, or in the presence of noise, when, due to 
randomness, the finer steps for Ai > F$i, merge into the 
broad oscillations of panel (b) ; 

(d) As the temperature vanishes, the adiabatic current 
in Eq. (24 1 gets suppressed for Ai < F^r and enhanced 



elsewhere[see Fig. |8ja)]. For F 3 r < A\ < F 3 l the par- 
ticle can only move to the right, so that j grows almost 
linearly with A\, while for A\ > F%l the particle can 
drift in both directions, thus making j decrease. 

Rectifiers, like a rocked ratchet, can work against an 
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FIG. 5 Rocked ratchet. Rectification current j versus the 
noise intensity D, see Eq. Q, panel (a), and versus ac mod- 
ulation strength A = Ai, panels (b) and (c) for the ratchet 



potential in Eq. (23 1 with L = 1. The relevant simulation 
parameters are: (a) Ai = 0.5 and fii = 0.01,1,2.5,4 and 
7 (top to bottom); (b) D = 0.1 and Qi = 0.01,1,4 and 7 
(top to bottom). In both panels (a) and (b) the solid line for 
Sli = 0.01 coincides with the adiabatic limit in Eq. (24 1; (c) 
deterministic regime with D = 0, and fii = 0.01 (dashed line) 



and 0.25 (solid line). After Bartussek et al. (19941 



external dc drive. On adding a dc component Aq to the 
sinusoidal force, namely, on applying the drive ^ in- 
stead of we can determine the value of Aq, termed 
stopping force, such that for a certain choice of the other 
perturbation parameters, Ai, f2i and D, the net current 
vanishes. 

Rocked ratchets with massive particles exhibit strong 
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inertial effects also capable of reversing their current. A 
noiseless inertial rocked ratchet is naturally subject to 
developing chaotic dynamics; indeed, its current turns 
out to be extremely sensitive to the initial conditions 



(Borromeo et al. 


2002 


Jung et al. 


1996 


Machura et al. 


2007a 


Mateos 2000 


20031. However, an appropriate 



evel suffices to stabilize the rectification properties 
of an inertial ratchet and makes it a useful device with po- 



tential applications in science and technology (Machura 
et aL[|2004b| |Marchesoni et aL[|2006a| ). 

In the notation of Sec. |II.A.1[ slightly generalized 
to account for the asymmetry of the potential in Eq. 
(23 1, we define two additional pairs of thresholds, besides 



the depinning thresholds F 3R ^: the rcpinning thresholds 
Fiji l an d the underdamped transition thresholds i^iJ.i; 
with Fm < FiL for i = 1,2,3 [the suffix i?, L denoting 
the orientation of the drive]. Under the adiabatic condi- 
tion fii <C 7, luo, the following approximations have been 

~ ( |2002| ); |Marchesoni et al. 



(2006a) 



obtained in Borromeo et al. 



(a) in the zero-temperature limit, T — * 0: 

j{Ax) = j R (Ai), F 2R <Ai< F 2L 

= jR{Ai)-j L (A 1 ) A 1 >F 2L (25) 

and j(Ax) = for Ax < F 2 r. The right/left current 



components are j R ,L(Ax) =2JA\- F^ RL /(TT-fL). 
(b) in the noiseless regime, T = 0: 

j(Ax) = jn(Ax), F 3R <Ax< F 3L 

= 3r{Ax) - Jl{Ax), Ax>F 3L (26) 

and j(Ax) = for Ax < F 3R . Here, j R , L (Ax) = 



A\ - F? 



1R,L 



A x ~ F%rl\/(k1 l )- Note that j(Ax) 



in Eq. ( 26 1 is discontinuous at the depinning thresh- 
olds F 3Ri l as expected for the response curve of any un- 
derdamped rocked ratchet operated in the deterministic 



regime (see Sec. II. A. 1) 



2. Pulsated ratchets 

Let us consider now the overdamped process of Eq. 
with Fx(t) = and F 2 (t) = A 2 sgn[cos(fl 2 t + fa)]. 



As discussed in Sec. |II.C.2[ modulating the amplitude of 
a symmetric potential V(x) in the presence of an uncorre- 
cted, time-symmetric perturbation £ (t) does not suffices 
to induce rectification. This state of affair changes when 
F 2 {t) couples, instead, to an asymmetric substrate, like 
our reference ratchet potential ( 23 1 . For the sake of sim- 
plicity, we set A 2 = 1, so that the effective substrate 
potential [1 + F 2 (t)]V(x) appears to switch on and off at 
every half period T 2 /2 — 7r/f2 2 - The interplay between 
time modulation and spatial asymmetry generates a non- 
zero drift current, which, at low (but not too low!) fre- 
quencies, is oriented in the positive direction. Devices 
operated under similar conditions are termed pulsated 
ratchets. 



Rectification by a pulsated ratchet can be explained 
qualitatively by looking at the cartoon of Fig. [6] The 
probability density of a particle initially placed at the 
bottom of a potential well (top), is mostly confined to 
that well, as long as T 2 /2 is not exceedingly longer than 
the escape time (t(L,0)}, which is the case at low tem- 
peratures. When, during the following half cycle, the 
substrate is switched off (middle), the particle, still sub- 
ject to noise, diffuses freely with Einstein's constant D 
and, for T 2 3> L 2 /D, its probability density approaches 
a Gaussian spanning over many a unit cell. By switching 
on again, the substrate cuts such a probability density 
into smaller peaks of different size, one for each well sur- 
rounding the initial one. For potential barriers skewed 
like in figure, l + > Z_, wells on the right are expected to 
be more populated than wells on the left, hence a positive 
net current. This is also the case of the potential in Eq. 



( 23 1 , which has positive natural orientation when pul- 
sated, and negative orientation when rocked. The mech- 
anism illustrated in Fig. [6] can work even in the presence 
of dc force, Fi(t) — A , pointing in the opposite direc- 
tion, Aq < 0. Such an upward directed motion is clearly 
powered by the thermal fluctuations £(£). 

The directed particle current is clearly bound to vanish 
in the adiabatic limit, O2 — *■ 0, when thermal equilibrium 
is approached. A similar conclusion holds true for very 
fast substrate modulations, fl 2 ^ ^0 ( m rescaled units). 



In perturbation analysis (Savel'ev et al. 2004c I, one finds 



the noteworthy result that the current decays to zero in 
both asymptotic regimes remarkably fast, namely like 
in the slow modulation limit, and fl 2 , in the fast modula- 
tion limit respectively. Moreover, for intermediate mod- 
ulation frequencies, the current in a pulsated ratchet is 
generally not oriented in its natural direction; current in- 
versions are possible when the forcing frequency matches 
some intrinsic relaxation rate of the process. 
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FIG. 6 (Color online) Pulsated ratchet. Rectification current 
j = (x) versus (dimensionless units), for the ratchet po- 



tential of Eq. ( 23 1 with L — 2-k and Vo = 1. Other simulation 
parameters are: A2 = 0.8 and D — 0.5 The low frequency 
f2|, and high frequency, f2^" 2 asymptotic limits are denoted 
by straight lines. Courtesy of Marcello Borromeo. 



Finally, pulsated ratchets can be operated under very 
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general amplitude modulations F^i), as well. This 
ratchet effect proved robust with respect to (a) modifica- 
tions of the potential shape (Reimann 20021; (b) changes 
in the switching sequence. A generic discrete modula- 
tion signal Fz(t) is characterized by different residence 
times, namely, t on in the "on" state, +A2, and f ff i n the 
"off" state, — A 2 , respectively. For periodic modulations, 
ton + t s = T2 defines the so-called duty cycle of the de- 
vice (Ajdari and Prost 1992 Bug and Berne 19871; for 



random modulations, instead, t on and t a g must be inter- 
preted as the aver age residence times i n the respective 
"off" state (Astumian and Bier 1994 Faucheux 



"on" 



et al. 1995 1. In both cases the modulation is asymmetric 
for t on ^ t B', (c) replacement of pulsated with flashing 
substrates. The substrate is made switch, either period- 
ically or randomly in time, among two or more discrete 
configurations, which do not result from the amplitude 



modulation of a unique substrate profile (Borromeo and 
Marchesoni] p98] [CEenl \TWf\ |Gorman et al.\ |1996| |Lee 



et al. 2004 1 . The rectification process is controlled by the 
spatial asymmetry of the single substrate configurations 
or by the temporal asymmetry of the switching sequence, 
or by a combination of both; (d) modulation of the tem- 
perature (Bao 1999 Reimann et al. 1996). Modulating 



the intensity of the ambient noise £(t) or its coupling 
to the device, corresponds to introducing the time de- 
pendent temperature T(t) = [1 + F2(t)]T, for an appro- 
priate control signal F^i) with amplitude smaller than 
one. Pulsating the potential amplitude Vq or the tem- 
perature T yield very similar rectification effects, being 
the net current of a pulsated ratchet mostly controlled by 
the barrier-to-noise ratio; (e) the addition of constant or 
spatially modulated damping terms (Luchsinger (2000); 
|Suzuki and Munakata| p003 1; see also Reimann (20021). 
Inertial effects make the rectification mechanism sensi- 
tive to the particle mass, so that selective trans port and 
segregation of mixed species becomes possible (Derenyi 



|and Ajdari 1996 Derenyi and Vicsek 19951 



3. Correlation ratchets 

To better understand the role of asymmetry in the rec- 
tification of a nonequilibrium process, we consider now 
an overdamped Brownian particle with Langevin equa- 
tion |3]) diffusing in the ratchet potential of Eq. (23) 
subject to a zero-mean, colored Gaussian noise F = 7j(t), 
with 



Q 



■ exp 



(27) 



The substrate asymmetry is capable per se of rectifying 
the correlated fluctuation r)(t), even in the absence of 
external modulations or thermal noise. The early predic- 
tion of this transport effect, based on simple perturbation 



arguments (Bartussek et al. 


1994 


Doering et al. 1994 


Luczka et al.\ |1995| |Magnascol |1993| Millonas and Dyk-| 


man 1994 


1, kindled a widespread interest in the ther- 



path-integral analysis (Bartussek et al. 1996) led later 



to the following low-noise estimates of the correlation 
ratchet current: 

(a) in the weak color limit, uj^t <C 1, 



(28) 



where j±(r) = tk exp[— t 2 Qc±/(D + Q) 2 ] and tk is the 
equilibrium Kramers rate for the particle to exit a poten- 
tial well through one side, r K = uj q uj± exp[— AV/(D + 
Q)]/(2tt). Here, xq and x± are the extremal points of 
V{x) defined as in Fig. g <J}\ = V"(x ), w± = \V"(x ± )\, 
and c ± =Q[V"(x)]VTx)dx; 
(b) in the strong color limit, r — » 00, 



J(r) 



Jo 



1 — exp 



QL 
2D 2 t 



IJ 



, (29) 



where jo is a positive definite constant. Equation (28 1 
can be regarded as the difference of two exit currents 
3± ( T ) > respectively through x + (forward) and X- (back- 
ward), in the presence of colored noise (Millonas and 



Dykman 1994) 



J 



0Cr> 




FIG. 7 Correlation ratchet. Rectification current j versus r 
and Q for D = 0.1. The parameters of the ratchet poten- 
tial, Eq. (23 1, are as in Fig. [4] More data can be found in 



Bartussek et al. ( 1996 1 



The orientation of the current of a correlation ratchet 
is extremely sensitive to the substrate geometry. In the 
strong color limit, the condition £_) > between the 
two ratchet lengths guarantees that the current is posi- 
tive, no matter what the choice of V(x); all correlation 
ratchets work in the natural direction of the correspond- 
ing rocked ratchet. 



modynamics of molecular motors. A more sophisticated 



In the weak color limit, the current in Eq. ( 28 ) has the 
same sign as the difference c_ — c + = J Q [V" (x)]V (x)dx, 
which, in turn, depends on the detailed profile of V(x), 
regardless of sgn[Z + — /_]. Numerical simulation of cor- 
relation ratchets with potential given in Eq. ( |23[ ), con- 
firmed that, consistently with the inequality c_ — c + > 0, 
j(r) is positive definite, as illustrated in Fig. [7] as 
j(r — * 0) = j(r — ► 00) = 0, optimal rectification 
is achieved for an intermediate noise correlation time, 
ujqT ~ 1. Modifying the potential profile, e.g., by adding 
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driven by two rectangular wave forms, 
(d) gating mechanism, Eq. (19 1. The 



FIG. 8 (Color online) Net current in a piecewise linear potential with amplitu de AV 
Eq. (181: (a) one-frequency rocked ratchet; (b), (c) harmonic mixing of Sec. II.C.l 
ratchet potential is: V{x) = xAV/L 2 for -L 2 < x < and V{x) = xAV/Li for < x < Li, with Li = 0.9 and L = 1\~AV = 1 
in (a)-(c) and = 2 in (d). Panel (a): Response curve j(A) = (x)/L to a low-frequency rectangular force with amplitude A 
at zero temperature D = (dashed curve), and low temperature D/ AV = 0.05 (solid curve). Panel (b): Numerical simulations 
for a doubly rocked ratchet with (pi = <f>2 — 7r and fii 

(green line and green crosses). The off-set Javg, Eq. (301, is indicated by the green line; the spikes of Eq 
integer- valued odd harmonics are marked with green crosses (x); Panel (c): Adiabatic approximation for <pi — 
(main panel) and <j>\ = 3n/2, <j>2 = tt/2 (inset). In both cases A\ = 3, A<i = 2, D = 0.6; Panel (d): Numerical simulations 
versus adiabatic approximation, Eqs. (21 1 and (31 1, for a rocked-pulsated ratchet with A\ = 4, A 2 = 0.5 and fii = 0.01; noise 
level: D — 0.4. Main panel: <pi — 4>2 = tt (adiabatic approximation); inset: simulation (open circles) versus the fully adiabatic 
approximation (x) for cj>i = n and (f>2 = 0. After Borromeo and Marchesoni (2005b I. 



0.01 (open circles) versus adiabatic approximation, Eqs. (211 and 1 30 1 



(21 1 at some selected 
2 = 3tt/2 



relevant drift currents get modified as follows: 



an appropriate higher order Fourier component without 
changing ?_)_—/_, suffices to reverse the sign of c_ — c+ (i) Current off-sets. At variance with Sees. II.C.l and 

and, thus, introduce at least one current inversion (Bar- |II.C.2[ two ac drives applied to an asymmetric device 

w also 



tussek et al. 1996 1 



Correlation ratchets have the potential for technolog- 
ical applications to "noise harvesting" , rather than to 
nano-particle transport. The low current output and its 
extreme sensitivity to the substrate geometry and the 



particle mass (Lindner et al. 1999 Marchesoni 1998) 



makes the design and operation of correlation ratchets 
as mass rectifiers questionable. However, the asymme- 
try induced rectification of nonequilibrium fluctuations, 
no matter how efficient and hard to control, can be ex- 
ploited by a small device to extract from its environment 
and store the power it needs to operate. 



4. Further asymmetry effects 



are expected to induce a net current j avg = 
for incommensurate frequencies f2i, ^2- This is a mere 
ratchet effect that in most experiments is handled as a 
simple current off-set. This conclusion is apparent, for 
instance, in the case of low-frequency rectangular input 
waves, where the device output current can be expressed 
as a linear combination of two known ratchet currents 
dSavel'ev et~al\ |2004a| |Savel'ev et a/.[|2004| ), namely: 
(a) Harmonic mixing current (Sec. II.C.l I. When the 
drive is a double rectangular wave, Eq. (18 1, j a vg can 



be regarded as the incoherent superposition of currents 
from two rocked ratchets driven by rectangular waves 
with amplitudes A\ + Ai and \A\ — A<i\^ respectively, i.e. 



1 



3a 



{jaQA! - A 2 \) + j (Ax + A 2 )]; 



(30) 



The rectification mechanisms introduced in Sec. III. CI 
apply to any periodic substrate, independently of its spa- 
tial symmetry. In the presence of spatial asymmetry, the 



(b) Gating current (Sec. II.C.2I. The rocked-pulsated 
device of Eq. (19 1 can be regarded as the incoherent 



superposition of two ratchets with potentials V±(x) 



1G 



(1 ± A2)V(x), respectively, both rocked with amplitude 
A\\ hence 



1 



(31) 



In Eqs. (30) and (311 jo,±(A) are the net currents of Eq. 
( 24 I for a rocked ratchet driven by a rectangular wave 
with amplitude A and vanishing frequency, the suffixes 

0, ± referring to the regular, high/low amplitude poten- 
tial configurations, V(x) and V±(x). 

(ii) Asymmetry induced mixing. As anticipated in Sec. 
II.C.1[ a double rectangular wave is not capable of rec- 
tifying a Brownian particle in a symmetric potential, 
that is, (x) = also under harmonic mixing conditions, 
^2/^1 = m/n with m,n coprime integers and m + n 
odd. This is no more the case if V(x) is asymmetric. 
The non-zero odd moments of the process x(t), deter- 
mined by the substrate asymmetry, generate additional 
harmonic couplings of the "odd" harmonics of the drive, 

1. e., for fia/^i = (2m — l)/(2n — 1) with 2m — 1 and 
2n — 1 coprime, as shown in Fig. [8] The total rectifica- 
tion current, including both the incommensurate term of 
Eq. (30 1 and this new harmonic mixing spikes, was cal- 
culated analytically in |Savel'ev et al.\ ( |200 z 



et al. (20041: 



Savel'ev 



J = Ja 



(-l) m+ >(A n , m ) 



(2m 



l)(2n- 
Eq. 



1) 



Jodd, 



(32) 



(211 and j odd 



where p(A„ m ) is defined after 
l\j Q (\A x -A 2 \)-j Q {A x + A 2 )]. 

The competition between current off-set, item (i), and 
asymmetry induced spikes, item (ii), may cause surpris- 
ing current reversals for special values of the driving fre- 
quencies, as observed in several experimental setups (cf. 
Sees. Ill IV and VII.Bl. Moreover, replacing the rectan- 
gular waveform, Eq. (18), with a more conventional lin- 



ear superposition of two sinusoids of the same frequency, 
Eq. (16 1, leads to an even more complicated interplay 
of nonlinearity and asymmetry induced harmonic mixing 
( |Savel'ev~ £ al. , 2004a I. A special limit of bi-harmonically 
driven rocked ratchet is discussed in Sec. III. E. 21 

A brilliant demonstration of the combination of har- 
monic mixing and asymmetry effects in the context of 
particle transport at the nanoscales, has been recently re- 
ported by Kalman et al. ( 2007 1 , who managed to drive di- 



lute ions through conical nanopores by applying a bihar- 
monic rectangular voltage. Harmonic mixing and gating 
currents could be separated, as predicted by the theory, 
and the resulting commensuration effects are displayed 
in Fig. [9] 



E. Efficiency and control issues 



The analysis of Sec. |II.D| suggests that the output of a 
real ratchet device is hard to control experimentally, let 
alone to predict. As a matter of fact, if we use a ratchet 
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FIG. 9 (Color online) Net current through a single conical 
nanopore for the sum of two voltage signals applied across 
the pore, 0.5sgn[sin(fiit)]+0.5sgn[sin(fi2i+7r)], with different 
^2/^1. The measurements were taken in a 0.1M KC1 solution 
at pH 8.0. For further details see |Kalman et al\\2QQ7 ). 



to rectify an assigned signal, the only tunable variables 
are the temperature T, the damping constant 7 (or the 
mass m), and the substrate profile V(x). Under certain 
experimental circumstances these quantities may prove 
not directly accessible to this purpose or inconvenient to 
change. This is why control of particle transport in small 
devices sometimes requires introducing auxiliary signals, 



auxiliary particle species (Sec. V.C I, or even the use of a 



feedback control scheme to optimize the Brownian motor 



current ( Craig et al. 


2008a|b||Feito and Cao 


2007 


2008 


Son et al. 


2008). 



1. Optimization 

The most common definition of efficiency for a loaded 
Brownian motor is 



770 = (i)A /(P)i 



(33) 



where (x)Aq is the average mechanical work done per 
unit of time against a working load Aq and (P)m is the 
average net power pumped into the system by the ex- 
ternal drives, no matter how applied. The quantity r) 
has been interpreted in terms of macroscopic thermo- 
dynamics by Sekimoto ( |1998 1 (for further developments 
on issues of energetics and efficiency see also |Parrondo| 
and de Cisneros (2002)), and references therein. In their 
scheme, a ratchet operates like a Carnot cycle, where the 
lower temperature is determined by the thermal noise, 
Eqs. ([2] and [4]), and the higher temperature is related 
to the magnitude of the external modulation. Operating 
far from thermal equilibrium, the efficiency is typically 
smaller than the upper limit of a thermal reversibly op- 
erating Carnot cycle. 

The ensuing question, whether a ratchet can be oper- 
ated at the maximal Carnot efficiency, spurred an intense 
debate (see for a review, Sec. 6.9 in Reimann (2002)). 
The issue has been clarified by an analysis performed 
within the validity regime of the framework of linear irre- 
versible thermodynamics by den Broeck (20071: The key 
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towards obtaining maximal Carnot efficiency for a Brow- 
nian motor is zero overall entropy production. This can 
be achieved by use of architectural constraints for which 
the (linear) Onsager matrix has a determinant equal to 
zero, implying vanishing (linear) irreversible heat fluxes. 
Maximizing efficiency subject to "maximum power" also 
leads to the same condition of a vanishing determinant 



coefficient also has a maximum (Sec. II.A.3 1 . That lead 



of the Onsager matrix (i.e perfect coupling) (|den Broeck 



2005), yielding the Curzon-Ahlborn limit (Curzon and 



Ahlborn 19751, which in turn at small temperature dif- 
ference just yields half the Carnot efficiency. For the 
archetype Smoluchowski-Fcynman Brownian motor de- 
vice (Feynman et al. 1963 Smoluchowski 19121 the ef- 



ficiency at maximum power has recently been evaluated 
for different coupling schemes by Tu ( |2008 1: Interest- 
ingly enough , the typical upper bound set by |Curzon 
and Ahlborn ( |1975 1 can then even be surpassed. 

In spite of the combined efforts of theorists and exper- 
imenters, as of today the question remains unanswered. 
At present, ratchet devices operating under controllable 
experimental conditions hardly achieve an efficiency i]q 
larger than a few percent. 



On the other hand, definition (33 1 of efficiency is not 



always adequate to determine the optimal performance 
of a Brownian motor. Firstly, r/Q assumes that work is 
being done against a load Aq, which is not always the 
case. For example, one clearly finds a vanishing effi- 
ciency whenever no load Aq) is present. Secondly, the 
only transport quantifier used in Eq. ( |33[ ) is the drift 
velocity (v), whereas the fluctuations of xjt) — v(t), i.e., 
the variance a 2 — (v 2 ) — (v) 2 , are also of practical im- 
portance. Note that in this Section we use v(t) to denote 
the particle velocity in Eq. ^ as a proper stochastic 
process. If a v > (v), and even more so if a v 3> (v), the 
Brownian motor can possibly move for some time against 
its drift direction (v). 

A load-independent rectification efficiency, r] r , has 



been introduced by Suzuki and Munakata (2003) and 



then generalized by Machura et al. (2004b I. They com 



puted T) r as the ratio of the dissipated power associated 
with the directed motion of the motor against both the 
friction and the load, and the input power from the time- 
periodic forcing. The result assumes the explicit form 



Vr 



(vKA + (v)) 
\( V 2) + ( V )A -D \ 



(34) 



where Dq is the free diffusion coefficient of Sec. |II.A.3| 
This definition holds evidently also for Aq — 0, while nu- 
merical evidence indicates that (v 2 ) > D , consistently 



with the inequality r\ r < 1. Note that the definition (34) 
assumes that x(t) is a damped process with finite 7 (no 
matter how large). In this way r\ r accounts explicitly also 
for the power dissipated as velocity fluctuations during 
the rectificatio n process; in particular, it increases upon 
decreasing a v . Machura et al. (2004b I noticed that, for 



a rocked ratchet, a v exhibits pronounced peaks in corre- 
spondence with the rectification thresholds (large 
damping) and F2R.L (low damping), where the diffusion 



to the general conclusion that transport by a Brownian 
motor can be optimized by operating away from activa- 
tion thresholds, in regimes of large net currents, where 
the velocity fluctuations are intrinsically small. 




FIG. 10 (Color online) Mobility versus ^2/^2 for a vibrated 
ratchet with A\ = 0.5, fii = 0.01, 4>\ — 4>2 = 0, and different 
values of the noise intensity D (see legend, main panel). All 
simulation data have been obtained for SI2 = 10, but the black 
crosses where we set D = 0.12 and f^2 = 20. Bottom inset: 
simulation data for fi as in the main panel with an additional 
curve at D = 0.6. Top inset: fj, versus D for A2 — 0, A\ = 
0.5, and Qi — 0.01; circles: simulation data; solid curve: 
adiabatic formula (11.44 ) in ( |Risken| [l984[ ). After |Borromeo| 
and Marchesoni (2005a I. 



2. Vibrated ratchets 



Besides the remarkable exceptions presented in Sec. 
IV[ transport control in a rocked ratchet cannot be ob- 
tained experimentally by deforming the substrate poten- 
tial "on demand" . As shown by |Borromeo and March- 



csom 



(2005a I, this goal can more easily be achieved by 
means of an external tunable signal. Let V(x) be as- 
signed a fixed profile, say, the standard double-sine in Eq. 
(23 1. The ratchet current can be still varied by inject- 



ing an additional control signal with frequency f^, that 
is, by replacing the harmonic drive in Eq. (|8| with the 
bi-harmonic drive in Eq. (16); in the adiabatic regime 



f2 2 3> ^1 the system response is very different than in 
Sec. |II.D.4| A rocked ratchet operated under such con- 
ditions is termed vibrational ratchet. 

Following the perturba tion approach of Baltanas et al. 
( |2003[ ); peckmanj ( p000| ; |Landa and McClintock| ( |2000[ ), 
the variable x(t) in Eq. (|3j can be separated as 
x(t) — ► x(t) + ip(t): in the remainder of this Sec- 
tion x(t) will represent a slowly time-modulated stochas- 
tic process and ip(t) the particle free spatial oscillation 
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ijj(t) — ^osin(n 2 i + </ ) 2), with = ^2/^2- On averaging 
out ip(i) over time, the Langevin equation for the slow 
reduced spatial variable x(t) can be written as 



where 
V(x) 



-Vo 



V(x)+A l coa(Sl 1 t + ^ 1 ) + ^(t), (35) 

(36) 



Jo(ipo) siax + -Jo(2ipo) sin 2x 



and Jo (a;) is the Bessel function of zero-th order. The 
symmetry of the effective vibrated potential ( 35 1 is re- 



stored if one of its Fourier components vanishes, namely 
for A 2 /n 2 = |ji,ii,|i2, |i3,i2, \h, \fa,hi---> where 
jx = 2.405, j 2 = 5.520, j 3 = 8.654, j 4 = 11.79, 
j5 = 14.93 . . . , are the ordered zeros of the function 
Jo{x); correspondingly, the ratchet is expected to vanish, 
as confirmed by the simulation data displayed in Fig. [10] 
Not all zeros of this sequence mark an inversion of the 
ratchet current. For instance, for A2/Q2 < \jx the cur- 
rent is oriented in the natural direction of the effective 
potential, Eq. 



(|36|); for |ji < A 2 /fl 2 < h, the coeffi- 



cient of sin 1x changes sign and so does the ratchet natu- 
ral orientation, or polarity; on further increasing A 2 /fl 2 
larger than ji , the sign of both Fourier coefficients in Eq. 
(36 1 become reversed: this is equivalent to turning V(x) 



upside-down (beside slightly re- modulating its profile), 
so that the polarity of V(x) stays negative. Following 
this line of reasoning one predicts double zeros (i.e. no 

current inversions) at tp = ji, J2, J3, J4, 

This control technique can be easily applied to the 
transport of massive Brownian particles both in symmet- 
ric ( Borromeo and Marchesoni 2007a I - see inset of Fig. 



10 - and asymmetric devices (Borromeo and Marchesoni 



2006 ) . Likewise, the use of a delay in a feedback contro 
signal ([Craig et q?.[ |2008a|b| |Feito and Cao[|2007 



2008 



Son et al. 2008p can, via synchronization mechanisms 



efficiently improve the performance of Brownian motor 
currents; a scheme that as well can be implemented read- 
ily in experimental flashing ratchets ( Craig et al. 2008a I . 



III. TRANSPORT IN NANOPORES 

Membranes in biology encase cells and their organelles 
and allow the compartmentalization of cellular processes, 
thereby operating far from thermal equilibrium, a condi- 
tion that is essential to life. Membranes separate two 
phases by creating an active or passive barrier to the 
transport of matter between them. As a first classifi- 
cation, membranes can be divided into biological and 
artificial membranes. The latter term is applied to all 
membranes made by man with natural, possibly mod- 
ified, materials and with synthetic materials (synthetic 
membranes). Synthetic membranes can be further di- 
vided into organic (made with polymers) and inorganic 
membranes (made with alumina, metals, etc.). 

Transport across a membrane occurs through channels, 
or pores (Berezhkovskii et al. 2003 Kolomeisky 20071. 



In several cases, the membrane transport properties must 
be regarded as a collective effect, where the function of 
an individual channel is influenced by the presence of 
(possibly diverse) neighboring channels. This is the case, 
for instance, of ion pumps ( |Im and Roux 2002 L auger | 
1991 1 in cellular membranes (Sec. III.A) or for coupled 
ion channels which experience a common transmembrane 
voltage. Single molecule techniques, however, allow for 
the study and characterization of rectification properties 
of individual entities, such as single ion channels. Nowa- 
days, synthetic membranes with assigned pore density 
and patterns are made commercially available. Moreover, 
by means o f increasing sophistica ted growth methods, 
irradiation (Fleischer et al.\ 19751, and nanofabrication 



techniques (|Dekker[ |200 7 Hcaly et al. 
|2003P 



2001 Storm etal\ 



2007 Li et al. 



the cross section of membrane 
pores can be modulated along their axis. As a result, 
transport in artificial ion channels fabricated from asym- 
metric single pores of the most diverse geometries has be 
come accessible to experiments ( Healy et al. 2007 Siwy 



and Fulihski 2004), see also in Sec. III.B 



In this Section we focus on devices where the size of 
the transported particles is comparable to the pore cross 
section. Transport in these devices can be analyzed in 
terms of the one-particle ID mechanisms illustrated in 
Sec. [TTJ Larger pores, commonly used to channel fluids 
(liquid or gaseous) or colloidal particles suspended in a 
fluid, arc considered in Sec. 



VI in the context of microflu- 



idic devices. 



A. Ion pumps 

As remarked already before, biology teaches us use- 
ful lessons that can guide us in the design of artificial 
Brownian motors. Biological membranes knowingly form 
lipophilic barriers and have embedded a diverse range of 
units which facilitate the selective movement of various 
ionic and polar segments, or the pumping of protons and 
electrons across channel like membrane openings. These 
biological nanodevices make typically use of electrochem- 
ical gradients that enable them to pump a species against 
its concentration gradient at the expense of yet another 
gradient. Such devices thereby convert nondirectional 
chemical energy, from the resource of the hydrolysis of 
adenosintriphosphate (ATP), into directed transport of 
charged species against electrochemical gradients. 

Although many details of the underlying mechanism 
are far from being understood, these machines apparently 
make use of mechanisms that characterize the physics of 



Brownian motors ( Astumian 2007 1 . A characteristic fea- 



ture is that changes in the binding affinity at selective 
sites in the transmembrane region must be coupled to 
conformational changes which, in turn, control motion 
into the desired direction. Thus, in contrast to rocked 
or pulsated ratchets, where the particle-potential inter- 
action acts globally along the whole periodic substrate 
landscape, the transport mechanism at work here can be 



19 



better described as 


an "information" ratchet (|Astumian 


and Derenyi| 


1998 


Parrondo and de Cisneros 2002 


). In- 



deed, the effective potential bottlenecks to transduction 
of Brownian motion get modified locally according to the 
actual location of the transported unit; as a result, in- 
formation gets transferred from the unit to the poten- 
tial landscape. For example, this scheme can be use d to 
model the pumping of Ca ++ ions in Ca ++ -ATPase (Xu 



et al. 2002|. Related schemes have been invoked also in 
the theoretical and experimental demonstration of pump- 
ing of Na + - and K + -ions via pulsed electric field fluctu- 
ations in Na- and K-ATPase ( Freund and Schimansky^| 



Geier[ |1999| |Tsong[ |2002| |Xie et al.\ |1994[ ), or for the 



operation of a catalytic wheel with help of a ratchet- 



like, electro-conformational coupling model (Rozcnbaum 



et al. 2004 1 



Yet another mechanism can be utilized to pump elec- 
trons in biomolecules. It involves nonadiabatic elec- 
tron tunneling in combination with asymmetric, but non- 
biased, non equilibrium fluctuations, as proposed bylGoy- 
chuk (2006). The nonequilibrium fluctuations originate 



from either random binding of negatively charged ATP or 
externally applied asymmetric, but non-biased, electric 
stochastic fields. Likewise, unbiased noncquilibrum two- 
state fluctuations (telegraphic noise) can induce direc- 
tional motion across an asymmetric biological nanopore, 
as numerically investigated for an aquaglyceroporin chan- 
nel, where water and glycerol get transported by means 



of a rocked ratchet mechanism (Kosztin and Schultcn 



2004) 



Finally, we stress that in such realistic complex biolog- 
ical channels the physical implications of an externally 
applied control are often difficult to predict. This is 
due to the multi-facetted consequences any control action 
may have in terms of chemical variations, conformational 
changes, polarization effects, and the alike. Moreover, 
at variance with single-molecule type experiments, single 
channel recordings are not easily accessible when dealing 
with biological molecules. Nanopores of lesser complexity 
are thus synthetic nanopores which can be fabricated by 
use of bottom-down nanoscience techniques - the theme 
we review next. 



B. Artificial nanopores 



With the recent advances of track-etching (Fleischer 



et al. 



et al. 



2001 



1975| and silicon technologies (Dckkcr 
" 2003| 



Storm et al. 



2007 Li 



charge transport in a 
single nanopore became experimentally accessible. This 
is a substantial leap forward with respect to ion pump 
and zeolite transport experiments, where experimental 
data are taken over a relatively high channel density. 
Fabricated nanopores in polymer films and silicon ma- 
terials are being investigated in view of their potential 
applications as biomimctic systems, that is, for modeling 
biological channels, and as biosensors. 



Siwy and coworkers (Constantin and Siwy 2007 Siwy 
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FIG. 11 (Color online) Ion current recorded experimentally 
by |Constantin and Siwy| ( |2007| > at 0.1M KC1, pH= 5.5 
through a single conical nanopore with diameters were 5 nm 
and 1000 nm, respectively. The pore rectification power is 
n = 217 at 5V. Inset: geometry of a conical nanopore with 
schematic representation of surface charge distribution creat- 
ing a bipolar nanofluidic diode. Figure provided by Zuzanna 
Siwy. 



eTaTl[2005||Siwy and Fulihskij |2002l | Vlassiouk and Siwy| 



20071 have recently prepared a nanofluidic diode, which 
had been predicted to rectify ion current in a similar 
way as a bipolar semiconductor diode rectifies electron 
current ( |Daiguji et al] |2005[ ). This diode is based on a 
single conically shaped nanopore track-etched in a poly- 
mer film with openings of several nm and 1 /im, respec- 
tively (see sketch in Fig. 111. The surface charge of the 



pore is patterned so that two regions of the pore with 
positive and negative surface charges create a sharp bar- 
rier called the transition zone. This nanofluidic diode is 
bipolar in character since both positively and negatively 
charged ions contribute to the measured current. Ma- 



jumdar and coworkers (Karnik et al. 2007) fabricated a 



similar nanofluidic diode with a sharp barrier between a 
positively charged and a neutral side of the pore. The 
presence of only one type of surface charge causes the 
latter device to be unipolar. 

Ion rectification was achieved by applying a longitudi- 
nal ac voltage (Fig. 11 1; the system thus operate as a one- 
cell ID rocked ratchet (Sec. II.D.ll, where the spatial 



asymmetry is determined by the interaction of a single 
ion with the inhomogeneous charge distribution on the 
pore walls. The rectification power of the pore is defined 
as the ratio of the ionic currents recorded for positive and 
negative driving voltages, i.e., r](V) = \I(V)\/\I(— V)|. 
Of course, due to their asymmetric geometry, conical 
pores can rectify diffusing ions also for uniform wall- 
charge distributions (Siwy et al. 20051; however, the 



corresponding rj factor would be at least one order of 
magnitude smaller than reported here. 

As a serious limitation of the present design, it is not 
possible to control the rectification power of a given con- 
ical nanopore without introducing changes to its built-in 
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electro-chemical potential. An alternate approach has 
been proposed by |Kalman et al. (2007), where two su- 
perposed rectangular voltage signals of zero mean were 
used to control the net ion current through a nanopore of 
preassigned geometry. By changing the amplitude, fre- 
quency and relative phase of these signals these authors 
used the gating effect of Sec. |II.C.2| to gain control over 
the orientation and the magnitude of ion flow through 
the pore. Their experimental data were found in ex- 
cellent agreement with the theoretical prediction of Eq. 
( 32 1 . The magnitude of ion current variation achieved 



by asymmetric signal mixing was comparable with the 
incommensurate current off-set Eq. ( 30 1 , which means 



that the nanopore diode could be operated disregarding 
the details of its intrinsic rectification power. 

We conclude this Section on artificial nanopores an- 
ticipating that asymmetric micropores etched in sili- 
con membranes also work as microfluidic ratchet pumps 



for suspended colloidal particles (Kettner et al. 2000 



Matthias and Miillcr 20031. For instance, entropic 



effects on the rectification efficiency of the conical 
nanopores of Siwy and coworkers have been analyzed 
by Kosinska et al. (20081. However, in-pore diffusion 



in a liquid suspension requires a fully 3D analysis of the 
pumping mechanism, which sets the basis for the fabrica- 



tion of more complicated ionic devices ( van der Heyden 



et al. 


2004 


Stein et al. 


2004) 



is reviewed in Sec. |VlJ 



2004). This category of devices 
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FIG. 12 (a) (Color online) Electric detection of individual 
DNA molecules with a nanopore. A constant voltage bias in- 
duces a steady-state ionic current through a single nanopore 
(left panel, from simulation); adding DNA to the negatively 
biased compartment causes transient reductions of the ionic 
current (right panel, from experiment). This reduced conduc- 
tance is associated with the translocation of DNA through the 
pore, which partially blocks the ionic current. After |Aksimen-| 
|tiev et aT](|2004[). (b ) Schematic diagram of the entropic trap 
inlHan et al.\ <|1999h. 



C. Chain translocation 

The opposite limit of composite objects passing 
through a much narrower opening is often called "translo- 
cation " . Using the experimental set-up sketched in Fig. 
lla), 



Kasianowicz et al. ( 1996 1 measured, for the first 



time, the blockage currents of single-stranded RNA and 
DNA clcctrophoretically driven through the a trans- 
membrane pore. Recent developments of this technique 
demonstrated single nucleotide resolution for DNA hair- 



pins (|Ashkenasy et al. \ |2005| |Gerland et al] |2004| |Ver 



coutere et al. 2003 1 , thus raising the prospect of creat- 



ing nanopore sensors capable of reading the nucleotide 
sequence directly from a DNA or RNA strand, see for 
recent comprehensive reviews flHealyj |2007| |Movileanu 



2008; Zwolak and Di Ventra 20081. Nowadays, translo- 



cation mechanisms arc investigated in both protein chan- 



nels (mainly the bacterial a- hemolysin pore (Kasianowicz 



et al.\ 1996[ )) and synthetic pores. Both these approaches 
have advantages and disadvantages ( I )ekker[ 2007 Healy 
20071. Protein channels can be engineered with almost 
angstrom precision, but the lipid membrane in which 
they may be incorporated is very mechanically unsta- 
ble. Synthetic pores, on the other hand, offer robustness 
to the system, which allows to better characterize the 
physical aspects of translocation phenomena. 

As a main difference with ion transport, the transloca- 
tion of a long polymer molecule in a ID device involves 
entropic effects, which become important when the open- 
ing cross section grows comparable with the radius of gy- 
ration, Rq, of the polymer. These effects were predicted 



in (Arvanitidou and Hoagland 19911 to account for the 



conformation changes a chain undergoes to move past a 
conduit constrict i on and, finally, o bserved by |Han and] 
Craigheadlp000|;|Han et aZ.|(fT999 1 in an artificial chan- 



nel. As a model pore-constriction system, these authors 
fabricated a channel consisting in a periodic sequence of 
regions of two different depths, as shown in Fig. 12 'b). 
The thick regions were 1mm deep, i.e., comparable with 
the i?o of the double stranded DNA molecules they used 
in their experiment, whereas the depth of the thin region, 
90nm, was much smaller than Rq. 

The thick regions act like "entropic traps" , as the DNA 
molecules are entropically prevented from entering the 
thin regions. For the same reason, a chain caught in be- 
tween two traps, tends to fall back into the trap that 
contains the most of it. In the presence of an external 
electric field, escape from a trap is initiated by the intro- 
duction of a small portion of DNA into an adjacent thin 
region, just enough to overcome the escaping activation 
barrier. This initiation process is local in nature, and 
the energy barrier does not depend on the total length of 
the trapped DNA molecule. Once a DNA molecule is in 
the transition state (once a proper length of "beachhead" 
is formed), it readily escapes the entropic trap, regard- 
less of the length of the remaining molecule in the trap. 
Quite counter intuitively, Han and coworkers found that 
escape of DNA in longer entropic traps occurs faster than 
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in shorter ones. 

A theoretical interpretation of these results was given 
by |Park and Sungl (119991), who treated the dynamics of 



a flexible polymer surmounting a ID potential barrier as 
a multidimensional Kramers activation process (Hanggi 



et al. 1990). To determine the activation free energy, 
Park and Sung computed the free energy of the polymer 
at the transition state. For a small-curvature barrier, 
the polymer keeps its random coil conformation during 
the whole translocation process, giving rise to essentially 
the same dynamics as that of a Brownian particle. For 
a large-curvature barrier, on the other hand, a confor- 
mational transition (coilstretch transition) occurs at the 
onset of the barrier crossing, which significantly lowers 
the activation free energy and so enhances the barrier 
crossing rate. As the chain length varies, the rate shows 
a minimum at a certain chain length due to the compe- 
tition between the potential barrier and the free energy 
decrease by chain stretching. Synthetic nanopores can 
thus be used as Coulter counter devices to selectively 
detect single DNA molecules, resolving their length and 
diameter. 

Spatial asymmetry can rectify translocation through 
synthetic and biological pores, alike. This is the case, 
for instance, of a hydrophobic polymer translocating 
across a curved bilayer membrane. Extensive simula- 



tion (Baumgartner and Skolnick 19951 showed that the 



polymer crosses spontaneously and almost irreversibly 
from the side of lower curvature to the side of higher 
curvature, so as to maximize its conformational entropy 



("entropy" ratchets, Sec. V.A.I). Moreover, at variance 



with artificial channels, in protein translocation through 
a biomembrane, the chemical structure of the pores can 
come into play by helping rectify the thermal fluctua- 
tions of the stretched molecule: When specific prede- 
termined segments of the protein cross the membrane, 
chemicals acting as chaperons bind on the segments to 



prevent their backward diffusion (Hartl 1996 Jiilichcr 



et al. 


1997 


Nigg 


1997) 



The ensuing chemical asym- 
metry then competes with the entropic asymmetry to 
determine the translocation current. 



D. Toward a next generation of mass rectifiers 

1. Zeolites 

Zeolites are three dimensional, nanoporous, crystalline 
solids (either natural or synthetic) with well-defined 
structures that contain aluminum, silicon, and oxygen 
in their regular framework. The silicon and aluminum 
atoms are tetrahedrally connected to each other through 
shared oxygen atoms; this defines a regular framework 
of voids and channels of discrete size, which is accessible 
through nanopores of well-defined molecular dimensions. 
The negative electric charge of the zeolite framework is 
compensated by (inorganic or even organic) cations or 
by protons (in the acidic form of the zeolites). The ions 



are not a part of the zeolite framework, and they stand 
in the channels. The combination of many properties - 
such as the uniform cross section of their pores, the ion 
exchange properties, the ability to develop internal acid- 
ity, high thermal stability, high internal surface area - 
makes zeolites unique among inorganic oxides and also 
leads to unique activity and selectivity. As a result, ze- 
olites can separate molecules based on size, shape, po- 



larity, and degree of unsaturation a mong others (Karger 
2008b| |Karger and Ruthven[|1992[). 



Nuclear magnetic resonance (NMR) provides direct ac- 
cess to the density and the mobility of molecules in ze- 
olites. By implementing a time-dependent NMR tech- 
nique, termed pulsed field gradient NMR, Karger and 
coworkers (see for a review |Karger| ( |2008a[ )) investigated 
the problem of particle diffusion in a narrow zeolite 
nanopore. For instance, they observed that, owing to 
the constrained pore geometry (0.73nm across), the mean 
square displacement of CF4 molecules (0.47nm in diam- 
eter) diffusing in zeolite AIPO4-5, increases linearly with 
the square root of the observation time rather than with 
the observation time itself. That was an early experimen- 
tal demonstration of the single-file diffusion mechanism 
introduced in Sec. |II.A.4| In the meantime, however, 
real zeolite crystals have been repeatedly found to deviate 
notably from their ideal textbook structure ( Schemmert 



et al. 19991, so that, their pore cross-sections ought to 
be regarded as longitudinally corrugated. To what extent 
this may affect particle diffusion is still matter of ongo- 



ing research ( Karger et al] 2005; Taloni and Marchesoni 
20061. 



Most zeolite structure types exhibit 3D pore networks. 
The most studied example are a synthetic zeolites of type 
MFI. Their pore network is formed by mutually inter- 
secting straight (in crystallographic y-direction) and si- 
nusoidal (in crystallographic re-direction) channels. Al- 
though there is no corresponding third channel sys- 
tem, molecular propagation has been observed also in z- 
direction. Experimental data are consistent with a simple 



law o f correlated diffusion anisotropy (Fenzke and Karger 
1993), a 2 z /D z — a x /D x + a^/D y , where at and Di, with 
1 = x,y,z, denote respectively the lattice and the diffu- 
sion constants in the i-direction. This means that the 
molecular "memory" is shorter than the mean traveling 
time between two adjacent network intersections. 

With the advent of synthetic zeolites and zeolitic mem- 
branes, researchers are fascinated by the option that the 
existence of different channel types within one and the 
same material may be used for an enhancement of the 
performance of catalytic chemical reactions. As the diffu- 
sion streams of the reactant and product molecules tend 
to interfere with each other, rerouting them through dif- 
ferent channels may notably speed up a catalytic reactor; 
hence, the idea of reactivity enhancement by "molecular 



bauer et al. 2000) 



traffic control" (Derouane and Gabelica 1980 Neuge 
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2. Nanotubes 

Another interesting category of artificial nanopore are 
carbon nanotubes (Dresselhaus et al. 19961. Single- 



walled carbon nanotubes are cylindrical molecules of 
about 1 nm in diameter and 1-100 /im in length. They 
consist of carbon atoms only, and can essentially be 
thought of as a layer of graphite rolled-up into a cylin- 
der. Multiple layers of graphite rolled in on themselves 
are said to form a multi-wall carbon nanotubc. The elec- 
tronic properties of nanotubes depend strongly on the 
tube diameter as well as on the helicity of the hexag- 
onal carbon lattice along the tube (chirality). For ex- 
ample, a slight change in the pitch of the helicity can 
transform the tube from a metal into a large-gap semi- 
conductor, hence they potential use as quantum wires in 



nanocircuits (Collins and Avouris 2000 Dckker 19991 



An even wider range of geometries and applications be- 
came available recently with the synthesis of various 
oxide nanotubes (Remskar 20041. With their hollow 



cores and large aspect ratios, nanotubes are excellent 
conduits for nanoscale amounts of material. Depending 
on the filling material, experimenters have thus realized 
nanoscale magnets, hydrogen accumulators, thermome- 
ters and switches. 

Nanotubes also provide an artificial substrate for con- 



trollable, reversible atomic scale mass transport. Regan 
et al. ( 2004 I attached indium nanocrystals to a multi-wall 



carbon nanotube and placed it between electrodes set up 
in the sample chamber of an electron microscope. Apply- 
ing a voltage to the tubes, they observed that the metallic 
particles at one end of the tube gradually disappeared, 
while those at the other end grew. While the details 
of the underlying driving mechanisms (thermo- versus 
electro-migration) remain unclear, they concluded that 
the voltage dictates the directionality of that nanoscale 
mass conveyor. Experimenters even succeeded to synthe- 
size carbon nanotubes encapsulating metallic atoms and 
characterize the electro-mechanical properties of such 
nanochannels (see e.g., Gao and Bando (2002)). In the 



next future, carbon nanotubes will also be combined to 
form molecular "gears" , whose feasibility has been proven 

A conceptual 



so far only by simulation (Drcxler 



1992) 



example is provided by a double- walled carbon nanotube 
consisting of two coaxial single-walled nanotubes with 
different chirality, immersed in an isothermal bath. In 
the presence of a varying axial electrical voltage, this 
system would exhibit a unidirectional ratchet-like rota- 
tion as a function of the chirality difference between it 
constituents ( |Marchesoni 1996| . 

Nanotubes have also been used to realize prototypes 



of thermal diodes for phonon transport (Chang et al. 
2006). With this concept in mind, one can further devise 



a phonon Brownian motor aimed to ratchet a net heat 
flux from "cold to hot" , as numerically demonstrated in a 
recent work by Li et al. ( 2008 1 . This class of devices has 



the potential to allow an efficient control of heat fluxes 
at the nanoscales. 



In spite of the recent advances in nanotechnology, (in- 
ner or outer) transport along nanotubes is still controlled 
by external gradients. Nevertheless, nanotube based 
ratchets, though not immediately available, are like to 
be one of the next frontiers in artificial Brownian motor 
research. Moreover, synthetic nanotubes are potential 



building blocks also for nanofluidic devices (Holt et al. 



|2006|, as discussed in Sec. VI 



IV. COLD ATOMS IN OPTICAL LATTICES 



Optical lattices are periodic potential for atoms ere- 



ated by the interference of two or more laser fields ( 


Gryn- 


berg and Robilliard 


2001 


Jessen and Deutsch 


19961. In 



near resonant optical lattices the laser fields produce si- 
multaneously a periodic potential for the atoms and a 
cooling mechanism. The optical potential for an atom in 
an optical lattice is given by the light shift, or ac Stark 
shift, of the atomic energy level that acts as the ground 
state of an assigned optical transition. The simplest case 
is represented by the J g — 1/2 ~ > J e = 3/2 atomic tran- 
sition in a ID configuration. 

Let us consider two counter-propagating laser fields, 
detuned below the atomic resonance, with orthogonal 
linear polarizations and same intensity and wavelength 
A. Their interference results into a spatial gradient with 
polarization ellipticity of period A/2. This in turn pro- 
duces a periodic potential for the atom. For instance, the 
atomic hyperfine ground states \g, ±) = \ J g = 1/2, M = 
±1/2) experience periodic potentials in phase opposition 
along the direction x of light propagation, namely 

V±{x) = y(-2±cosfc2:), (37) 

where k — 2n/X and the depth of the potential wells 
Vq scales as Il/A, with II the total laser intensity and 
A the detuning from atomic resonance. As the laser 
fields are near to resonance with the atomic transition 
J g = 1/2 — > J e = 3/2, the interaction with the light fields 
also leads to stochastic transitions between the Zeeman 
sublevels \g, ±) of the ground state. The rate of this tran- 
sitions can be quantified by the photon-atom scattering 
rate V which scales as Il/A 2 . It is therefore possible to 
vary independently the optical lattice depth Vq and V by 
changing simultaneously 1^ and A. 

The stochastic transitions between ground states also 
lead to damping and fluctuations. The damping mecha- 
nism, named Sisyphus cooling, originates from the com- 
bined action of light shifts and optical pumping which 
transfers, through cycles of absorption/spontaneous 
emission involving the excited state J e , atoms from one 
ground state sublevel to the other one. Moreover, the 
stochastic transitions between the two potentials V±(x) 
also generate fluctuations of the instantaneous force ex- 
perienced by the atom. In conclusion, the equilibrium 
between cooling and heating mechanisms determines a 
stationary diffusive dynamics of the atoms, a process con- 
fined to the symmetric ground-state optical lattice V±(x). 
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Random amplitude pulsations of one symmetric po- results in Gommers et al. (20051; Schiavoni et al. (20031 



tential, Eq. (37), do not suffice to produce rectification. 
The only missing clement to reproduce the rocked set-ups 
of Sees. 11. C and II. D is the additive ac force F(t). In 
order to generate a time-dependent homogeneous force, 
one of the lattice beams is phase modulated, and we will 
indicate by <p(t) its time-dependent phase. In the lab- 
oratory reference frame the phase modulation of one of 
the lattice beams results into the generation of a moving 
optical lattice V±(x — 4>(t)/2k), reminiscent of t he peri- 
odic sieves in ( Borromeo and MarchesonH |2007a I . In the 
accelerated reference frame x — > x — 4>(t)/2k, the optical 
potential would be stationary; however, an atom of mass 
m experiences also an inertial force in the z-direction pro- 
portional to the acceleration of the moving frame, namely 
F{t) = (m/2k)(j>(t). This is the homogeneous ac drive 
needed to rock a cold atom ratchet. 

The above described harmonic mixing/rocking ratchet 
set-up was recently used by Renzoni and coworkers 



( Gommers et al. 


2005 


2007 


2006 2008 Schiavoni et al. 


2003 


) to investigate experimentally the relationship be- 



tween symmetry and transport in ID atom traps. We 
discuss here separately the experimental results for two 
different cases: (1) a bi-harmonic driving including two 
harmonics at frequencies £Iq and 2f2o, and (2) a multi- 
frequency driving obtained by combining signals at three 
different frequencies. 



A. Bi-harmonic driving 



These authors (Gommers et al. 2005 Schiavoni et al. 



2003) generated a bi-harmonic drive, Eq. (161, with 
fi 2 = 2f2i and <fii = 0, and searched for a harmonic 
mixing current proportional to sin(02 — 4>o). At variance 



with Eq. 17 a phase lag 4>o was introduced to account 
for the finite atom dissipation. Indeed, by generalizing 
our argument of Sec. |II.C.1| one can easily prove that 
increasing the damping makes <po drop from down to 
— 7r/2. In the experiment, the atom dissipation could be 
tuned continuously, without changing the optical poten- 
tial constant Vo, by varying simultaneously II and A at 
a constant /l/A ratio; thus, harmonic mixing was exper- 
imentally accessible in both limiting regimes of zero and 
infinite damping. 

The experimental results of Fig. [13] for cesium atoms 
cooled in the /iK range, clearly demonstrate the mecha- 
nism of harmonic mixing at work. In fact, the net atom 
velocity plotted in panel (a) is well fitted by v/v r — 
A/\ sin(</>2 — </>o), where the dependence of 4>o on V is 
as in panel (b), and is a characteristic function of the 
system (Gommers et al. 20051. For the smallest scat- 



tering rate examined in the experiment, no current was 
generated at (f>2 — 1ft, with I integer, as expected from 
Eq. (17 1. On the other hand, the magnitude of the phase 
shift </>o increases with increasing the scattering rate, thus 
causing current generation also for </> 2 = It:. Besides con- 
firming the notion of harmonic mixing, the experimental 



also suggest that symmetry breaking can be controlled by 
dissipation. Closely r elated results have been reported by 
|Ustinov et al. ( |2004| ) in their attempt to control fluxon 
ratcheting in Joscphson junctions by means of harmonic 
mixing. 



B. Multi-frequency driving 

Recent experiments with multi-frequency driving 
(Gommers et al. 2007 20061 aimed to investigate the 



transition from periodic to quasi-periodic driving, and to 
examine how the analysis of Sec. |II.C| is modified in this 
transition. The multi-frequency driving was obtained by 
adding a sinusoidal component to the ac drive employed 
in the previous Section, i.e., 

F(t) = A x cos(^it)+A 2 cos(2f2i+ <j) 2 ) + A 3 cos(tt 3 t + (j) 3 ) 

. (38) 

For O3 /fii irrational the driving is quasi-periodic. 
Clearly, in a real experiment O3 /fii is always a ratio- 
nal number, which can be written as p/q, with p, q two 
coprime positive integers. However, as the duration of 
the experiment is finite, by choosing p and q sufficiently 
large it is possible to obtain a driving which is effectively 
quasi-periodic on the time scale of the experiment. 

Let us consider first the case of periodic driving, with 
il 3 /ili rational. Adding a third harmonic with phase 
constant <j) 3 =/= in Eq. ( 38 1 breaks the time symmetry of 



F(t) so that directed transport is allowed also for cf>o — 
and <j>2 = Itt- In other words, for tp 3 ^ the third 
driving component leads to an additional phase shift of 
the current as a function of 4>2- As a results, the current 
will retain its sinusoidal phase dependence sin(02 — </>o)j 
as above, with the difference that now <^ accounts for 
the phase shifts produced both by dissipation and by 
the added driving component. Following the discussion 
sketched in Sec. |II.D.4[ one easily concludes that in the 
quasi-periodic regime the third harmonic at frequency 
fl 3 is not relevant to characterize the time symmetry of 
the forcing signal, which is entirely determined by the 
bi-harmonic terms at frequency f2i and 2^2 ■ 

In later experiments ( Gommers et al. 2007 [2006 ) , the 
transition to quasi-periodicity was investigated by study- 
ing the atomic current as a function of </> 2 for differ- 
ent p/q. By increasing p and q the driving was made 
more and more quasi-periodic on the finite duration of 
the experiment, with the quantity pq an easily accessible 
measure of the degree of quasi-periodicity. The data for 
the average atomic current were fitted with the function 
f/^max = sin(02 — 4>o) and the resulting value for 4>o were 
plotted as a function of pq. For small values of pq, i.e., 
for periodic driving, the added component at frequency 
SI3 lead to a shift which strongly depends on the actual 
value of pq. For larger values of pq, i.e., approaching 
quasi-periodicity, the phase shift ip tends to a constant 
value independent of <p 3 , which coincides with the purely 
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FIG. 13 Cold atoms in an optical lattice, (a) Average atomic velocity as a function of the phase <j>2 for different values of 
r s = [uj v I (2-7r)] 2 / A (a quantity proportional to the scattering rate F'). The data are labeled by the lattice detuning A, as the 
vibrational frequency at the bottom of the wells was kept constant, uj v /2ti = 170kHz; the forcing frequency is £li/2n — 100kHz. 
The lines are the best fit of the data with the function v/v r — A a sin(<j) — cj>o), where v r = hk/m is the so-called atom recoil 
velocity, (b) Experimental results for the phase shift <j>a as a function of T s . After Gommers et al. (20051. 



dissipative phase-shift measured in the bi-harmonic driv- 
ing case of Sec. |IV.A| The experimental results confirm 
that in the quasi-periodic limit, the only relevant sym- 
metries are those determined by the periodic bi-harmonic 
driving and by dissipation. 



C. More cold atom devices 

Cold atoms in optical traps proved to be a playground 
for rectification experiments. 



(i) Gating effect. For instance Gommers et al. (20081 
modified the experimental set-up described above, to 
demonstrate experimentally a gating ratchet with cold 
rubidium atoms in a driven near-resonant optical lattice. 
As suggested in Sec. |II.C.2[ a single-harmonic periodic 
modulation of the optical potential depth with frequency 
^2 was applied, together with a single-harmonic rocking 
force with frequency VL\. The modulation of the optical 
potential depth Vq was obtained by modulating the in- 
tensity II of the laser beams. This also resulted in an un- 
avoidable modulation of the optical pumping rate, which 
affected only the fitting phase shift cf>Q. Directed motion 
was observed for rational values of f^/fii, a result due 
to the breaking of the symmetries of the system. 

(ii) Pulsated ratchets. Although a bona fide rocked 
ratchet for cold atoms in an optical lattice has not been 
realized (see e.g. |Ritt et al. (2006)), yet, a clever varia- 



tion of ID optical lattice described above allowed an early 



demonstration of randomly pulsated ratchet ( Mennerat- 
|Robilliard et al.\ [T999 ) . These authors set the polariza- 
tions of the laser beams at an angle 9 ^ tt/2 and applied a 
weak Zeeman magnetic field orthogonal to the optical lat- 
tice. The effect of the magnetic field consisted in remov- 
ing the degeneracy of the ground states by adding a A/2 
wavelength components to both potentials in Eq. ( 37 1 , 
which thus acquired different asymmetric profiles. As a 



consequence, the random optical transitions between the 
modified potentials V±(x) turned out to propel trapped 
cold rubidium atoms in the vertical direction, with the 
sign that depended on 6 and on the orientation of the 
magnetic field. 

In Sec. |II.D.2| we reported that rectification can occur 
on symmetric ID substrates that shift instantaneously 
back and forth in space with a fixed amplitude, namely, 
on substrate that are subjected to a time-discrete phase 
modulation (periodic or random, alike). In this case, 



breaking the supersymmetry condition (14 1 requires ap- 
propriate asymmetri c space-dependent switching rates 
(Gorman et al. |1996| . Based on this rectification scheme, 
Sjolund et al. (|2006[) realized an simple flashing ratchet 



for cold atoms. It consisted again of a /zK cold gas of 
cesium atoms switching between two symmetric ground 
state optical lattices, coupled via optical pumping. In 
the presence of induced friction, although small, and for 
appropriate laser detunings and intensities of the laser 
beams, the degeneracy of the ground states was removed 
by making one ground state, say V+, long lived and the 
other one, V-, short lived. If we further consider that 
the two optical lattices were shifted one relative to the 
other, V-(x) = V + (x — x ), then we conclude that the 
switching rates between potentials were state and posi- 
tion dependent with A; >+(x) ^> >-(x). In this setup, 

the atoms execute stationary time-asymmetric sequences 
of random jumps between V + and V_. During the time 
spent in the short-lived lattice V— , they experience a po- 
tential with an incline that depends on xq. Thus, their 
diffusion is strongly enhanced in one specific direction, 
and correspondingly reduced in the opposite direction. 
In the experiment, the spatial shift xq and the transition 
rates between the two optical lattices could be adjusted 
at will. A directed motion with constant velocity was ob- 
served in the absence of additional forcing terms, i.e., for 
F(t) — 0, except for specific system parameters, where 
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symmetry was restored. Moreover, these authors showed 



that their device, when operated in 3D (Ellmann et al. 



2003), can generate directed motion in any direction. 



V. COLLECTIVE TRANSPORT 

The studies on rectification mechanisms reviewed in 
Sec. [IT] have been conducted mostly for a single particle; 
however, often systems contain many identical particles, 
and the collective interactions between them may sig- 
nificantly influence the transport of particle aggregates 
(Aghababaie et al. 19991. For example, interactions 
among coupled Brownian motors can give rise to novel co- 
operative phenomena like anomalous hysteresis and zero- 
bias absolute negative resistance ( jReimann et al. 19991. 
The role of interactions among individual motors plays a 
particular important role for the occurrence of unidirec- 
tional transport in biological systems, e.g. see Sec. 9 in 
the review by Reimann (20021. 



In this section we shall restrict ourselves to the prop- 
erties of ID and 2D systems of point-like interacting par- 
ticles, although an extension to 3D systems is straight- 
forward. Let the pair interaction potential W(r) be a 
function of the pair relative coordinates r = (x,y), char- 
acterized by an appropriate parameter g, quantifying 
the strength of the repelling (g > 0) or attracting force 
(<? < 0), and by a constant A, defining the pair interac- 
tion length. An appropriate choice of W(r), with g > 0, 
can be used to model hardco re particles with diameter A 
(iSavel'ev et al.\ 120031 |2004b|). We term the pair interac- 



tion long-range, if A is larger than the average interparti- 
cle distance, and short-range in the opposite regime. In 
the absence of perturbations due to the substrate or the 
external drives, long-range repelling particles rearrange 
themselves to form a triangular |Abrikosov (19571 lattice. 



When placed on a disordered substrate or forced through 
a constrained geometry, such an ideal lattice breaks up 
into lattice fragments punctuated by point-like defects 
and dislocations and separated by faulty boundaries. In 
the presence of a drive, such boundaries act like ID easy- 
flow paths, or rivers, for the movable particles. This 
process, often referred to as "plastic" flow, is analyzed 



numerically in Reichhardt et al. ( 1998 1 



Rectification of interacting particles occurs as a com- 
bined effect of the configuration of the device and the 
geometry of its microscopic constituents. In this con- 
text the dimensionality of the underlying dynamics is also 
important: (1) Transport in some 2D geometries can of- 
ten be reduced to the ID mechanisms illustrated in Sec. 
[ll| (reducible 2D geometries); (2) Under certain circum- 
stances, however, a bona-fide 2D rectification may occur 
at a non-zero angle with the driving force (irreducible 
2D geometries). Finally, we remark that particle inter- 
action can be exploited to rectify particles of one species 
by acting on particles of another species, alone, either 
by an appropriate drive or a specially tailored substrate 
geometry. 



(a) 

V_^^^ 2d 




f\ L ^ x ^\ 




$..2;i 


o ° ocP °» °°° o 

° ° ° °„ °° o ° <?° ° ° 

8 %o o ° °°<*> o< 
o °°° °ooV* °o °°oV oj* 
= ° ' oooo°°° V«°«o % 
oo o o £5 o o 

° Wo -t . „ . 



FIG. 14 ID reducible asymmetric geometries: (a) Chan- 
nel with asymmetric profile y(x) — yo — yi[sin(2nx/L) + 
\ sin(47rx/L)]; the parameters y\ and j/o control the width 
2d and the bottleneck 2a; (b) Disordered pattern of circular 
defects generated with periodic distribution y(x) along the 
x-axis 



A. Asymmetric ID geometries 

We consider now two examples of 2D reducible geome- 
tries where the transport of massless repelling particles, 
subjected to an external ac drive, occurs as a collective 
effect, namely, under conditions incompatible with the 
rectification of a single particle. In particular, we show 
that collective effects make transport of interacting par- 
ticles even in the absence of an ad hoc ratchet substrate. 



1. Boundary effects 

Let us consider the asymmetric channel sketched in 
Fig. " 



14 'a) filled with n repelling particles per unit cell. 
The corresponding particle density is p — n/ai, where a; 
is the area of the channel unit cell. The walls are rigid and 
the particle-wall collisions are taken elastic and relatively 
short-range (interaction length not much larger than A). 
The profiles ±y(x) of the upper (+) and lower (— ) walls 
arc modeled by an appropriate double-sine function; 2d 
and 2a denote, respectively, the width of the channel and 
of its bottlenecks. Due to the repulsive interactions, the 
particles are pressed against the walls which corresponds 
to an effective asymmetric spatial modulation. As a con- 
sequence, when driven by an ac force, the particles are 
more likely to flow to the left than to the right. The ensu- 
ing rectification mechanism is reminiscent of a ID rocked 
ratchet with an inverted potential of Eq. ( 23 1 . Numerical 



simulations, besides supporting this prediction, clearly 
indicate an optimal or resonant, temperature regime in 
which the particle drift is maximized. This observation 
can be explained by noticing that at low temperatures, 
the ac drive causes a moving particle to migrate to the 
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center of the channel where it no longer interacts with 
the boundaries; while at high temperatures the driving 
force becomes irrelevant and thus the particle is no longer 
pushed through the channel bottlenecks periodically in 
time. The first simulation evidence of this phenomenon 
was produced by Wambaugh et al. (19991. Their results 



obtained for the special case of magnetic vortex channel- 
ing are discussed in Sec. |VII.A| 

At high enough particle densities, p, the net current in 
the channel is expected to be suppressed, as the repelling 
particles end up clogging the bottlenecks and thus ham- 
pering collective longitudinal oscillations. In the opposite 
limit of, say, n = 1 , the arguments above seem to rule out 
rectification, as a single oscillating particle would be con- 
fined in the inner channel of radius a. However, this con- 
clusion holds only in the absence of thermal fluctuations, 
T = 0. At finite T, reducing the diffusion in a 2D chan- 
nel to a I D proc ess implies defining an effective diffusion 
constant (|lI.A.3|, D(x)/D = [I + t/O) 2 ]" 1 / 3 , equivalent 



to a periodic asymmetric modulation of the temperature 
( Reguera et al. 2006 ) (a Seebeck ratchet in the notation 
of Reimann '( 2002 )). A particle crossing a channel bottle- 
neck perceives a lower effective temperature on the left, 
where the wall is steeper, than on the right, so that it 
gets sucked forwards; the opposite happens in correspon- 
dence to the largest channel cross-sections. As long as 
/ dx/D(x) 0, the oscillating motion of a single parti- 
cle can indeed be rectified, with the sign that dep ends 
on the details of the wall profile (|Ai and Liul 20061. No 



matter how weak, such a mechanism, termed "entropic" 



ratchet (Slater et al.\ 1997 1, supports the conclusion of 
Wambaugh et al. ( 1999^ that rectification in a channel 
occurs only at finite T, as a certain amount of noise is 
needed for the particle to explore the asymmetric geomet- 
ric of the device. Moreover, the argument above hints at 
the occurrence of an optimal channel density p, as de- 



tected in real superconducting devices (Sec. VILA I 



2. Asymmetric patterns of symmetric traps 



Olson et al. (20011 proposed a new type of 2D ratchet 
system which utilizes gradients of pointlike disorder, 
rather than a uniformly varying substrate potential. Let 
us consider a 2D sample containing a periodically grad- 
uated density of point defects, as sketched in Fig. [l4^b) . 
Each defect is depicted as a circular micro-hole, which 
acts as a symmetric short-range particle trap of finite 
depth. In real experiments, such defects can actually 
be created by either controlled irradiation techniques 



or direct- write electron-beam lithography (Kwok et al. 
2002 ) . The defect density pi (x) was chosen to be uniform 



along the vertical axis, and to follow a typical double- 
sine asymmetric profile of Eq. ( 23 I along the horizontal 



axis. Let us now inject into the sample repelling masslcss 
particles with average density p and interaction length A. 
The defect radius controls particle pinning by defects and 
was taken much smaller than A. For a sufficiently high 



particle-to-defect density ratios, the particles fill most of 
the pinning sites and create an effective repulsive poten- 
tial. If we further assume long-range particle pair interac- 
tions W(r), such a 2D potential surely gets insensitive to 
the details of the defect distribution; it only retains the 
periodic horizontal modulation of pi(x), thus resulting 
in a mean-field ratchet potential V(x), like in Eq. (23 1. 



The effective amplitude Vq is a function of at least three 
length scales: the interaction constant, the average par- 
ticle distance and the average defect spacing. A certain 
fraction of the particles does not become pinned at indi- 
vidual defects but, subjected to an applied ac drive, can 
move in the interstitial regions between pinning sites. Al- 
though the moving interstitials do not directly interact 
with the short-ranged defects, they feel the long-range 
interaction of the particles trapped at the pinning sites 
and described by the mean-field potential V(x). As a 
consequence, a horizontally applied ac drive can induce 
a longitudinal particle transport, as proven in |Qlson et al.\ 
(2001) by means of numerical simulation. 



Two conditions are instrumental to the onset of a rec- 
tification current: (1) a finite temperature, T > 0, and 
(2) a defect filling fraction close to or larger than one, 
both conditions being required for interstitial particles 
to get and stay unpinned. As in Sec. |V.A.1[ condition 
(1) implies the existence of an optimal rectification tem- 
perature. The dependence on the forcing frequency and 
amplitude are as discussed in Sec. |II.D.1| These ideas 



have been implemented to control transport of colloids 
and charge carriers in experimental setups where point 
defect gradients could be engineered at will (cf. Sees. 



VI.A 


VII.B 


and 


VIII.C 



B. 2D lattices of asymmetric traps 

Examples of substrates sustaining transverse rectifica- 
tion are asymmetric potential barriers, height q > 0, 
or wells, depth q < 0, cither isolated or arranged into 
ID chains and 2D lattices. Similar lateral displacement 
devices, also known as bumper arrays, have been pro- 
posed to separate particles by exploiting their mass and 
size d ispersion (Heller and Bruus 2008| Savel'ev et al. 



2005 1 . Consider for instance the pyramidal potential bar- 
riers/wells with isosceles triangular cross-section. This 
geometry, sketched in Fig. \TE\ is a generalization of the 



experimental setup by |Villegas et al. (J2003J. 

The substrates discussed here combine two types of 
asymmetry: the triangular shape of their building blocks 
and the asymmetry associated with the pyramidal struc- 
ture of each block. The latter asymmetry affects the 
motion of the particles only if the drive is strong enough 
to push them across the barriers/wells. |Savel'ev et al. 



(2005) have numerically simulated the dynamics of a 
gas of repelling massive particles driven across ID or 2D 
lattices of barriers/ wells for different parameters of the 
substrate lattice (Fig. [15] bottom panel), of the drive 
F(t) and of the particle pair potential W(r). Although 
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FIG. 15 (Color online) Top panel: asymmetric 2D arrays of 
potential energy barrier/ wells, top view. The parameters I, lo, 
a x , a y , A x , A y define the geometry of the array. The dashed 
and dash-dotted arrows represent trajectories perturbed by 
thermal noise or particle-particle repulsion. Bottom panel: 
(a) Transverse net velocity (v y ) versus Ai for massless inter- 
acting particles ac driven through a gapless triangular chain: 
Oy = 2, a x = 6 I = 1, A x = 4, A„ = 0, and T — 0. The 
repulsive potential W(r) is a wedge function with g = 1 
and half-width A = 0.05. Green solid circles are for weak 
short-range interacting particles, while blue open circles are 
for strong long-range interacting particles. Particle distribu- 
tions for strong long-range (b,c) and weak short-range (d,e) 
interactions for the Ai values indicated by arrows in (a). See 
for details |Savel'ev et~al\ ( |2005[ ) . 



Savel'ev et al. ( 2005 ) have shown that inertial effects of- 
ten enhance transverse rectification, we restrict our pre- 
sentation to the case of massless particles, certainly the 
most relevant for technological applications. We consider 
next two distinct operating regimes: 

(i) Deterministic setups. An ac force applied along the 
y-axis, i.e., parallel to the symmetry axis of the pyramid 
cross-section, was known to induce a longitudinal particle 



drift in the drive direction (see for more details in |Zhu 
et al. (20041 and Sec. VII). Indeed, driving a (distorted) 



lattice of repelling particles along the crystallographic 
axis of the substrate, oriented parallel to the symmetry 
axis of the substrate blocks, makes the system reducible 
to a mean-field ID dynamics in that direction, along the 
line of Sec. |V.A| Due to the pyramidal shape of its build- 
ing block, the reduced ID substrate is spatially asym- 
metric, which explains the reported longitudinal particle 
flows. 

(ii) Diffusive setups. Under appropriate conditions, 
instead, dc or ac forces applied perpendicularly to the 
symmetry axis may induce a transverse particle drift in 
the y-direction. For all geometries considered, the net 
velocity of a gas of non-interacting overdamped particles 
driven perpendicularly to the symmetry axis vanishes for 
T — > 0. Indeed, sooner or later each particle gets cap- 
tured in a horizontal lane between two triangle rows and 
then keeps oscillating back and forth in it forever. At fi- 
nite temperature, instead, fluctuations tend to push the 
particle out of its lane, thus inducing the net transverse 
currents reported in the earlier literature (Bier et al. 



2000||Duke and Austin[ |1998| |Huang et q?.||2004p . More 
remarkably, Savel'ev et al. ( 2005[ ) reported that the in- 
teraction among particles not only contributes to trans- 
verse rectification but actually plays a dominant role if 
the interaction length A, or the particle density are large 
enough. In particular, particle-particle interaction was 
shown to control transverse rectification for both weak 
short-range and strong long-range inter-particle forces 
(sec frame (a) of Fig. 
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Therefore, transverse rec- 
tification of interacting colloidal particles (short range) 
and magnetic vortices (long rang e) ar e expected to differ 
appreciably (see Sees. VIA and VII). 



To illustrate the key mechanism of transverse rectifi- 
cation, we consider only rectification by pyramidal bar- 
riers, q > 0, subjected to the square ac force F(t) = 
AiSgn[cos(fM )]. For a fully d e tailed analysis the reader 



is referred to Savel'ev et al. (2005). Let us consider 



first the geometry in frames (b)-(e) of Fig. 15 (bottom 



panel) , where the pyramids are stacked up in a close row 
(A y = 0). The particle interactions have a strong impact 
on the equilibrium particle distribution; a transition from 
an ordered lattice-like to a disordered liquid-like phase is 
displayed in frames (b)-(e). Nevertheless, the transverse 
currents (v y ) show a qualitatively similar A\ dependence 
for both weak short-range and strong long-range interac- 
tions; only the decaying tail is longer for the latter ones. 
The region of the linear growth of (v y ) for the case of 
weak short-range interaction corresponds to the regime 
when a constant fraction of particles (less than 1/2 be- 
cause of the geometry of the system) is rectified into the 
y-direction. This regime applies for increasing A\ until 
when particles start crossing over the triangles [frames 
(d) and (e)]. These results can be easily generalized to 
describe transverse rectification in any 2D lattice of tri- 
angular shaped barriers/ wells. 

Savel'ev et al. ( |2005 1 also simulated the case of a 2D 
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array of pyramids, top panel of Fig. [15] where a gap 
between triangles along the y-axis, A y > 0, turned out 



to make the net current sensitive to the particle-particle 
interaction length A. In the limit of low T, if A is smaller 
than a certain threshold value A c ~ A y , then the net 
current vanishes; if A exceeds A c then the current rapidly 
increases. This effect has been advocated to separate 
particles according to their interaction length. More pre- 
cisely, particles having an interaction length smaller than 
A c would pass through the array, or sieve of barriers sep- 
arated by Ay. In contrast, particles with a longer in- 
teraction length A > A c would be sifted sidewise. On 
connecting several such sieves with different gap A y , one 
can construct a device capable of separating the different 
fractions of a particle mixture. 



The overall conclusions of Savel'ev et al. (20051 do 
not change on replacing barriers with wells of the same 
shape, nor do so for lattices of asymmetric pins of differ- 



ent aspect-ratio and geometries 


((Bier et al 


2000 


Chep- 


elianskii and Shepelyansky||2005 Duke and Austin 


1998 


Ertas 


1998; 


Huang et al. 


2004 


Zhu et al. 


2003a 


20041. 



et al. (19981, who investigated the magneto-transport 
properties of a square lattice of triangular antidots of 
the type shown in Fig. [15] Antidot lattices are 2D elec- 
tron gases with appropriately placed sub-micron voids. 
These authors fabricated their lattices by electron beam 
lithography on shallow high electron mobility transis- 
tor structures, grown on semi-insulating gallium-arsenide 
substrates. In particular, they reported evidence that un- 
der far-infrared irradiation, electron sloshing in-between 
the antidot rows lead to a net electric current along the 
symmetry axis of the antidots. Significant transverse 
effects were observed in the presence of an orthogonal 
magnetic field. Moreover, in most applications to the 
electrophoresis of macromoleculcs, the particles moving 
through these sieve are suspended in the fluid where they 
diffuse, thus involving additional microfluidic effects (Sec. 

vib. 



C. Binary mixtures 

We address now the problem as how to induce and 
control the net transport of passive particles (target or A 
particles) , namely, particles that are little sensitive to the 



applied drives and/or substrates. Savel'ev et al. (2003 



2004b ) proposed to employ auxiliary B particles that (1) 
interact with the target species and (2) are easy to drive 
by means of external forces. By driving the auxiliary 
particles one can regulate the motion of otherwise passive 
particles through experimentally accessible means. 

Savel'ev and coworkers considered a mixture of two 
species of pointlike overdamped Brownian particles A 
and B at temperature T, diffusing on the ID periodic 
substrates described by the potentials Va(x) and Vb(x), 
respectively. Particles of type A interact pair-wise with 
one another as well as with the B particles via the po- 



tentials Waa and Wab (= Wba), while Wbb describes 
the interaction of the B pairs. The pair interaction is 
quantified by the tunable strengths g^, with i,j — A,B, 
whereas the interaction constants Ay play no significant 
role as long as they are conveniently small. 

To illustrate the rectification mechanisms in a binary 
mixture, we review in some detail the case of a pulsated 
device with oscillating temperature, i.e. a temperature 
ratchet studied originally by Reimann et al. (19961, see 
also Sec. 



II.D.2I 



The model can be made even simpler 
by assuming that (i) One subsystem (say, the auxiliary B 
particles) is subject to an asymmetric ratchet potential 
Vb(x), while the other one (the target A particles) is 
not, Va{x) = 0. This may happen, for instance, in a 
mixture of neutral and charged particles with the ratchet 
potential being produced by an electrical field; (ii) The 
interaction among particles of the same type is repulsive, 
i.e., g A A > 0, g B B > 0. 

Under these operating conditions, the B particles con- 
dense naturally at the minima of the asymmetric sub- 
strate potential [see Figs. 16 ^a) and (b), bottom pan- 



els]. In addition, if the B particles repel the A particles, 
Qab > 0, then the latter ones will accumulate in the re- 
gions where the density of the B particles is minimum, 
that is near the maxima of the substrate potential Vb{x). 
Vice versa, for attractive AB interactions, gAB < 0, the 
B particles concentrate around the minima of Vb(x). 
Therefore, the target A particles feel an effective poten- 
tial VX (x), which has opposite spatial asymmetry with 
respect to Vg S (x) for gAB > 0, and has the same asym- 
metry for gAB < [see Figs. 16 a) and (b), top panels]. 
Note that for low occupation numbers tia and tib, the 
bare potential Vb (x) is only marginally affected by parti- 
cle interaction, i.e., Vg S (x) ~ Vb(x). During the lower T 
half cycle, all particles are confined tighter around the 
minima of the relevant effective potential, while, dur- 
ing the higher T half cycle, the particles of both species 
diffuse more easily out of the V% B (x), V B )S (x) potential 
wells. As the asymmetry of the ratchet potentials (x) 
and Vg S (x) for repelling A and B particles is opposite, so 
is the orientation of their currents. On the contrary, two 
attracting A and B species drift in the same direction. 
This implies that the transport of both particle species 
can be effectively and separately controlled by regulating 
their number ua and ns per unit cell, without the need 
of tuning their substrate. 

As gAA, gBB > 0, the potential wells of V^{x) 
(Vg ff (x)) tend to flatten out when increasing the cor- 
responding density 115 (via)- This results in the decay 
of the associated ratchet current (ug) ((vb))- In con- 
trast, the ratchet asymmetry and current of one species 
may be enhanced by increasing the density of the other 
mixture component. For instance, in the case of AB at- 
tractive forces, the A particles tend to concentrate in the 
regions of higher B densities, that is around the minima 
of Vb, thus attracting even more B particles and making 
the Vq S (x) potential wells deeper. In short, one can en- 
hance the transport of the target particles by adding a 
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a) Repelling A-B species 



b) Attracting A-B species 
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FIG. 16 (Color online) Spatial dependence of the effective 
potentials Vj^ (x), Vg S (x) at different densities of the A and B 
particles. In both panels, particles of the same type repel one 
another; the interaction between particles of different species 
is repulsive in (a) and attractive in (b). There is no substrate 
for the A particles, Va{x) = 0, whereas the ratchet potential 
Vb{x) is piecewise linear: Vb(0 < x < 1+) = x/l+, Vb(1+ < 
x < I) = (1 - x)/l-, with I- = 0.2 and 1+ = 0.8. The other 
coupli ng parameters are gA A = ffss = \9ab\ = 1 and T = 1. 
[After |Savel'ev et al. | ( [20051 )]. 



certain amount of auxiliary particles. All these proper- 
ties have been reproduced analytically in the framework 



of the nonlinear Fokker-Planck formalism (Savcl'cv et al. 



2003 2004c) 



The influence of the interspecies interaction on the 
transport of a binary mixture has been investigated fur- 
ther by Savel'ev et al. (2004b I, where dc and ac external 



forces were applied at constant temperature, to the par- 
ticles either of one species, only, or of both, simultane- 
ously. Their main results can be summarized as follows: 
(1) with increasing the intensity of an applied dc driv- 
ing force, there is a dynamical phase transition from a 
"clustered" motion of A and B particles to a regime of 
weakly coupled motion (Marchesoni et al] 2006b); (2) 



by applying a time-asymmetric zero-mean drive to the B 
species only, one can obtain a net current for both the 
A and B species, even in the absence of a substrate; (3) 
when two symmetric ac signals act independently on the 
A and B particles, and only the particles of one species 
feel an asymmetric substrate, then, the two species can 
be delivered in the same or opposite direction by tuning 
the relative signal phase (for both attractive and repul- 
sive AB interactions) . Superconducting devices based on 
these two-species transport mechanisms have been real- 
ized experimentally in recent years (Sec. VII.Cl. 



VI. MICROFLUIDICS 



Microfluidics (Squires and Quake 2005) is playing 



an ever growing role in controlling transport of parti- 
cles, or even whole extended phases, on the micro- and 
nanoscales. The ability to manipulate the dynamics of 
liquids is crucial in various applications such as for a 



lab-on-a-chip technology. The Brownian motor concept 
has recently been invoked in different contexts to face 
this challenge. In this regard we remark that the iner- 
tial forces of small suspended particles are typically quite 
small ( |Purcell[|1977[ ). Pointlike particles can then be con- 
sidered to be advectively transported by the fluid velocity 
field at the particle's actual position. Then, for an in- 
compressible liquid the particle dynamics is volume con- 
serving and consequently no dynamical attractors emerge 
(Kostur et al. 20061. However, for extended particles the 



local velocity of a surface point need not coincide with 
the fluid velocity which would act at this point in the 
absence of the particle. Notably, for extended objects 
with internal degrees of freedom the volume of the state 
space is no longer conserved by the dynamics, thus giv- 
ing rise to attractors for stationary flow fields ( |Kostur| 



et al. 2006). Of course, Brownian diffusion provides an 
additional transport mechanism that must also be taken 
into account. 

We remark here that in spite of the intrinsic hydro- 
dynamical effects, microfluidic devices do not fall into 
the category of collective ratchets as denned in Sec. [Vj 
because here transported objects are not required to in- 
teract with one another. The suspension fluid still plays 
a central role: (a) it powers particle transport, and (b) it 
determines how the particle dynamics is coupled to the 
asymmetric geometry of the substrate. Genuine collec- 
tive ratchets are discussed in Sec. IVIII 



A. Transporting colloids 

Many physical examples and technological applications 
involve particles or molecules in solution that undergo a 
directed net motion in response to the action of a ratchet. 
There, the ratchet does not induce a mean flow of the 
solvent itself. For instance, colloidal particles or macro- 
molecules, suspended in solution, move when exposed to 
a sawtooth electric potential that is successively turned 



on and off (Roussclct et al. 1994 1. Electrolytic effects can 



be avoided by shuttling microsized Brownian polystyrene 
particles by optically trapping them with repeatedly ap- 
plied on-off cycles in an optical tweezer, thus mimick- 



ing a flashing Brownian motor i 


Faucheux et al. 


1995 


Faucheux and Libchaber 


1995 


Marquet et al. 


2002). 



flux reversal of diffusing colloidal spheres in an optical 



three-state thermal Brownian motor (Lee et al. 2004 1. 
This modus operandi of a Brownian motor can there- 
fore be put to work to pump or separate charged species 
such as fragments of DNA. A micromachined silicon-chip 
device that transports rhodamine-labeled fragments of 
DNA in water has been demonstrated with a flashing on- 
off Brownian motor scheme by Bader et al. (19991. Yet 



other devices are based on i deas and experimental real- 
izations of entropic ratchets (Chou et al. 1999 Duke and 
| Austin| |1998| |Ertasl |1998| |van Oudenaarden and Boxer 



|1999| |Slater et al] |1997| |Tessier and Slater} |2002[ ) which 
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FIG. 17 Top Panel: Scanning-electron-microscope picture of 
a single pore with ratchet longitudinal profile. The length 
of one period is 8.4/xm. The maximum pore diameter is 
4.8/im and the minimum pore diameter is 2.5/im. Bottom 
panel: Scanning-electron-microscope picture of a silicon wafer 
pierced by practically identical pores about 1.5mm apart and 
1mm in diameter. This picture illustrates the enormous po- 
tential for particle separation of a parallel 3D ratchet architec- 
ture. Figure provided by Frank Miiller, Max-Planck-Institute 
of Microstructure Physics, Halle, Germany. 



make use of asymmetry within geometric sieve devices to 
transport and separate polyelectrolytes. 

A similar concept is employed to selectively filter meso- 
scopic particles through a micro-fabricated macro-porous 
silicon membrane (Fig. 17 1, containing a parallel array of 



etched asymmetrical bottleneck-lik e pores (|Kettne r et al. 
|2000[ IMatthias and Miiller ~ 



2003 Miiller et al. 20001 



The working principle and the predicted particle flow 
for this micronuidic Browinan motor device are shown in 
Fig. 1 1 8 1 A fluid such as water containing immersed, sus- 
pended polystyrene particles is pumped back and forth 
with no net bias through the 3D array of asymmetric 
pores of Fig. [17] Such an artificial Brownian motor 
is thus kept far away from equilibrium by the periodi- 
cally varying pressure across the membrane. Due to the 
asymmetry of the pores, the fluid develops asymmetric 
flow patterns ( Kettner et al.\ 2000 1 , thus providing the 



ratchet-like 3D force profile in which a Brownian parti- 
cle of finite size can both: (i) undergo Brownian diffu- 
sion into liquid layers of differing speed, and/or (ii) be- 
come reflected asymmetrically from the pore walls. Both 
mechanisms will then result in a driven non-equilibrium 
net flow of particles. Note that the direction of the net 
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FIG. 18 Microfhudic drift ratchet of Kettner et al. (2000). 



Top panel: A macro-porous silicon wafer is connected at both 
ends to basins. The pores with their ratchet-shaped profile 
are sketched in dark grey. The basins and the pores are filled 
with liquid; microparticles of two different species are repre- 
sented. The fluid is pumped back and forth by the piston 
on the left. Figure provided by Christiane Kettner. Bottom 
panel: theoretically predicted net particle current, v e , versus 
the particle diameter for different driving frequencies and vis- 
cosities (relative to water) vr. Note in particular the very 
sharp velocity reversal around 0.5 /im. The results are for a 
pore of infinite length consisting of periodic units with a pore 
length of L = 6/xm, a maximum pore diameter of 4^i m and a 
minimum pore diameter of 1.6/xm, see Fig. [IT] (a). For further 
details see Kettner et al. (20001. 



flow cannot be easily guessed a priori. Indeed, the di- 
rection of the Brownian motor current is determined by 
the interplay of the Navier-Stokes flow in this engineered 
geometry and hydrodynamic thermal fluctuations. 

A detailed quantitative interpretation of the experi- 
mental data is plagued by several complications such as 
the influence of hydrodynamic interactions and, possibly, 
electric response effects due to residual charge accumula- 
tion near the boundaries. A most striking feature of this 
setup, however, is the distinct dependence of current re- 
versals on particle size. The sharply peaked current-size 
characteristics curves of this directed flow, that is, the 
theoretical flow vs. size prediction in Fig. [18] and the 
experimental current-pressure characteristics in Fig. |19| 
suggest a highly selective particle separation efficiency 
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( Kettner et al. 2000 ) . This microfluidic artificial Brown- Leidenfrost drops placed on a hot surface also with a 



ian device has recently been implemented in experiments 
with suspended polystyrene spheres of well defi ned diam- 
eters (100 nm, 320 nm, 530 nm and 1000 nm) by Matthias 
and Miiller (20031: Their experimental findings are in 



ratchet-like grooved profile (Fig. 20 1. The locally asym- 



metric geometry of the support induces a directed drop 
motion. This effect has been observed for many liquids 
and in wide temperature ranges within the film boiling 



good qualitative agrcment with the theory as shown in 

Fig. [T5J 

Remarkably, this device has advantageous 3D scaling 
properties: A massively parallel architecture composed 
of about 1.7 million pores illustrated in Fig. [T7] is ca- 
pable to direct and separate suspended microparticles 
very efficiently. For this reason, such type of devices 
have clear potential for bio-medical separation applica- 
tions and therapy use. 

The separation and sorting of cellular or colloidal parti- 
cles is definitely a topic attracting wide interest in differ- 
ent areas of biology, physical chemistry and soft matter 
physics. The powerful toolbox of optical manipulation 
(Dholakia et al. 2007 Grier[[2003[ ) uses the optical forces 
exerted on colloids by focussed laser beams to move and 
control objects ranging in size from tens of nanometers to 
micrometers. If the optical forces are sufficiently strong 
to rule transport, the stochastic Brownian forces play 
only a minor role, so that transport occurs as a determin- 
istic process with an efficiency close to unity, a circum- 



(Leidenfro st) regime (Linke et al. 
dari 



200G) 



2006 Quere and Aj- 



Moreover, micro-drops confined to the gap between 
asymmetrically structured plates, move when drop shape 
or wetting properties are changed periodically, for in- 
stance by vibrating the substrate, applying an on/off 
electric field across the gap or a low-frequency zero-mean 



electric field along the gap (Buguin et al. 2002 Gorre 



et al. 1996 1. In a related experiment the motion of drops 



on a global wettability gradient resulted to be strongly 



enhanced when a periodic force was applied (Daniel et t 
20041. In that work the ratchet aspect was determined 



by the intrinsically asymmetric shape of the drops and 
the hysteresis of their contact angle on the gradient sub- 
strate. 

The theory for particle transport operated by Brown- 



ian motors is presently well developed (Astumian 1997 
Astumian and Hanggi[ |2002[ |Hanggi and Bartussek 



stance known as (optical) peristalsis (Bleil et al. 2007 
Koss and Grier| |2003 1 . Such a scheme can then be im 



1996||Hanggi et al.\ [2005] |Jiilicher eTaTj [19971 [Reimann 
and Hanggi] 2002). In clear contrast, models for ratchet- 



plemented also to experimentally realize ratchet cellular 
automata capable of performing logical operations with 
interacting objects. All together, in combination with 
this new optical technology, colloids provides ideal model 
systems to investigate statistical problems near and far 



driven transport of a continuous phase are a rather scarce 
commodity, although the concept is thought to be impor- 
tant for emer ging microfluidic applications (Squires and 
Quake[|2005l. 



away from thermal equ ilibrium ( |Babic and Bechinger[ by John and Thiele (20071 and John et al 



2005 Babic et al. 2005) 



In recent work, a model for the flow of two immisci- 
ble layered liquids confined between two walls and driven 
by a flashing "wettability ratchet" has been put forward 

~ p008l 



B. Transporting condensed phases 

Most present applications use ratchet concepts to 
transport or filter discrete objects, like colloidal particles 
or macromolecules. However, ratchets may also serve to 
induce macroscopic transport of a continuous phase us- 
ing local gradients, only. One such realization, strongly 
related to the above cases of particle transport in a "rest- 
ing" liquid phase, is the Brownian motion of magnetic 
particles in ferrofluids subjected to an oscillating mag- 
netic field ( |Engel and Reimann 2004). In contrast to 
the cases reported above, here the motion is also trans- 
mitted to the solvent by viscous coupling. 

Another class of systems does not require colloidal par- 
ticles to drive the motion of the liquid phase. In a first 
example, a secondary large scale mean flow is triggered 
in Marangoni-Benard convection over a solid substrate 
with asymmetric grooves (Stroock et al. 20031. The di- 



In 

doing so, one employs the general concept of wettability 
comprising all effective interactions between the liquid- 
liquid free interface and the solid walls. Notably, any in- 
teraction that is applicable in a time-periodic, spatially 
periodic (but locally asymmetric) fashion, is capable to 
drive a flow, i.e. will result in directed transport. There 
exist several possibilities for experimental realizations of 
spatially inhomogeneous and time switching interactions. 
A practical setup consists of thin films of dielectric liq- 
uids in a capacitor with a periodic, locally asymmetric 
voltage profile, that can be periodically switched on and 
off. A spatially homogeneous electric field would desta- 
bilize the interface between two dielectric liquids. There- 
fore, an on-off ratchet with a constant lateral force can 
result in the de wetting-spreading cycle sketched in Fig. 



21 The process underlying this motor scheme may be 



termed "electro-dewetting" . 



C. Granular flows 



rection and strength of the mean flow can be controlled 
by changing the thickness of the liquid layer and the tem- 
perature gradient across the layer. 

Yet another intriguing situation involves self-propelled 



The concept of ratchet physics had found an early ap- 



plication in the field of granular matter (Jaeger et al. 
|1996[ ). Indeed, to explain the vertical grain size seg- 
regation under vibration, the so-called Brazil-nut effect 
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FIG. 19 (Color online) Parallel Brownian motors in a macroporous silicon membrane containing ca 
pores with 17 etched modulations each (see Fig 



1.7 million asymmetric 
When the pressure oscillations of the water 



17 top panel). Left panel 

are switched on, the measured photoluminescence signal and thus the number of particles in the basin located to the right, 
see top panel of Fig. |18| increases linearly in time, corresponding to a finite transport velocity. In contrast, for symmetric 
cylindrical-shaped pores no net drift is observed. Right panel: The experimentally observed net transport behavior strongly 
depends on the amplitude of the applied pressure and qualitatively shows the theoretically predicted current inversion ( Kettner 
et al. 20001. The pressure oscillations are toggled on and off each 60s. Other experimental parameters are as follows: the 
suspended luminescent polystyrene spheres in water are 0.32/xm across, the frequency of the pressure oscillations is 40Hz and 
the root mean square (RMS) of the applied pressure during the 'on' phase is 2000Pa. A control experiment using straight 
cylindrical pores with a diameter of 2.4/im showed no directed Brownian motor transport. [Image: Max-Planck-Institute of 
Microstructure Physics, Halle, Germany.] 
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FIG. 20 (a) Video sequence (time interval 32 ms) demon- 
strating droplet motion perpendicular to a thermal gradi- 
ent. A droplet of the refrigerant tetrafluorethane (R134a; 
boiling point Tb = 26.1°C) was placed on a room- 
temperature, horizontally leveled brass surface with a ratchet- 
like grooved profile [d = 0.3mm, s = 1.5mm). See also 
http//www. aip.org/pubservs/epaps. html (b) Droplet of liq- 
uid nitrogen on a fiat surface moving with a small initial ve- 
locity towards a piece of tape (shaded). Figure provided by 
Heiner Linke. 



(Rosato et al. 19871, Bug and Berne (19871 had initially Broeck 2007 



investigated directed diffusion in sawtooth-like inclined 
grooves. In the meantime, this problem spurred sev- 
eral computer simulations aimed at exploring the won- 
derful world of granular ratchets. Another unusual phe- 
nomenon, observed both experimentally and in computer 
simulations, is horizontal size segregation. It involves the 
occurrence of horizontal flows of granular layers that are 



vibrated vertically on plates whose surface profile con- 



sists of sawtooth-like grooves (Derenyi and Vicsek, 1995 
Rapaportl 12002 1. 



The spontaneous anelastic clustering of a vibroflu- 
idized granular gas has been employed, for instance, to 
generate phenomena such as granular fountains (con- 
vection rolls) and ratchet transport in compartmental- 
ized containers, where symmetry breaking flow patterns 
emerge perpendicular to the direction of the energy in- 



put (Brey et al. 2001 


Eggers 


van der Meer et al. 




2004 
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1999 



Granu 



Farkas et al. 



2002 



ar systems exhibit 



indeed a rich scenario of intriguing collective transport 
phenomena, which have been tested and validated both 



computationally and experimentally (van der Meer et al. 



2007) 



All these systems are kept away from thermal equilib- 
rium by vertical vibrations, while the presence of random 
perturbations play a role as a source of complexity, due to 
the intrinsic chaotic nature of granular dynamics. Closer 
in spirit of the present review are, however, those de- 
vised granular Brownian motors which are truly capable 
to tame thermal Brownian moti on (|Cleuren and Van den 



2007 2008). In those se- 



Costantini et al. 
tups, an asymmetric object of mass m (Fig. 22), is free to 
slide, without rotation, in a given direction. Its motion 
with velocity v is induced by dissipative collisions with 
a dilute gas of surrounding particles of average tempera- 
ture T. On each collision, only a certain fraction (called 
restitution coefficient) of the total kinetic energy of the 
object-particle system is conserved. At variance with the 
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FIG. 21 (Color online) Panels (a) to (c) illustrate the work- 
ing principle of a fluidic ratchet based on dewetting-spreading 
cycles, (d) illustrates the spatial asymmetric periodic interac- 
tion profile Q(x) responsible for the wettability pattern and 
(e) indicates the time dependence 0(t) of the switching in 
relation to the dewetting and spreading phases (a)-(c). The 
filled circle indicates a fluid element that gets transported 
during one ratchet cycle although the evolution of the inter- 
face profile is exactly time-periodic. For further details and 
resulting directed flow pattern see John et al. (20081. 



case of elastic collisions, here isotropy, implying (v) = 0, 
and energy equipartition, implying (v 2 ) — kT/m, do not 
apply: Dissipation is responsible for breaking the time 
reversal symmetry. As a result one finds that the asym- 
metric object undergoes (approximately) a directed ran- 
dom walk with nonzero average velocity, which can be 
modeled in terms of an effective biased Langevin dynam- 
ics for the velocity variable v. Such a Langevin dynam- 
ics assumes the form of an Ornstein-Uhlenbeck process 
(Hanggi and Thomas 1982 Riskenl 1984 1 governed by an 



effective Stokesian friction and an effective external force, 



both depending on the restitution coefficient ( 


Cleuren 


and Van den Broeck 


2007 


Costantini et al. 20071. The 



directed motion of the object is accompanied by a contin- 
uous flow of heat from the gas to the object, in order to 
balance the dissipation of the collisional processes; as a 
consequence the average temperature of the asymmetric 
object is lower than the gas temperature. 

Ever since the formulation of the second law of ther- 
modynamics there has been haunting debates about its 
validity at the small scales. In this context the ratchet 



gadget invented by Smoluchowski ( 1912 1 and later beau- 
tifully popularized by Feynman et al. ( 1963 1 points at 



the heart of a controversy, which concerns the very work- 
ing principle of Brownian motors: To be consistent with 
this law, no rectification of thermal fluctuations can take 
place at thermal equilibrium for which all dynamics is in- 
trinsically governed by the principle of detailed balance. 
A steadily working Brownian motors necessarily requires 
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FIG. 22 (Color online) Taming thermal motion. Panel 
(a) shoes a schematic representation of the Smoluchowski- 
Feynman rachet device, which is capable to rectify thermal 
Brownian motion of gases held at two different temperatures 
(Feynman et al. |1963[ ). An idealized version suitable for com- 
puter studies is presented in panel (b) with the ratchet re- 
placed by a triangular unit and the pawl replaced by a rigidly 
joined blade residing in the lower gas compartment. Panel (c) 
illustrates a 3D rotational construction. Figure provided by 



den Broek and den Broeck (2008) 



a combination of asymmetry and non-equilibrium, such 
as a temperature gradient. This point has been convinc- 
ingly elucidated in recent molecular dynamic studies by 



Van den Broeck and his collaborators (Van den Broeck 


and Kawai 


2006 Van den Broeck et al. 


2004 2005 


1, who 



variously stylized a Maxwell daemon operating far away 
from equilibrium (Fig. 
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Several computer versions 
of Smoluchowski-Feynman ratchets operating in a dilute 
gas of colliding granular particles, have thus been con- 
firmed to generate directed transport in the presence of 
a finite temperature difference. With the possibility of 
today's observation and manipulation of physical, chem- 
ical and biological objects at the nano- and meso-scale, 
such devices no longer represent a theorist's thought ex- 
periment, but are rather becoming a physical reality. It 
is right on these length scales that thermal fluctuations 
cannot be ignored, as they combine with non-equilibrium 
forces to yield constructive transport. 



VII. SUPERCONDUCTING DEVICES 

Magnetic vortices in type-II superconducting devices 
provide a flexible and experimentally accessible play- 
ground for testing our stochastic transport models. In 
recent years an impressive effort has been directed toward 
the design and operation of a new generation of vortex de- 
vices with potential applications to terahertz (THz) tech- 

|1997| jZitzmann et~aT\ |2002[ ) , 



nology ( jHechtfischer et al. 



fluxon circuitry ( Shalom and Pastoriza 2005 1 , and quan- 



tum computation ( You and Nori[|2005 l, to mention but a 
few. We anticipate that magnetic vortices are inherently 
quantum objects that, under most experimental condi- 
tions, behave like (massless, point-like) classical parti- 
cles. Genuine quantum effects in the vortex transport 
are briefly discussed in Sec. |VIII| 

A magnetic field applied to a type-II superconduct- 
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ing film, with intensity H comprised between the criti- 
cal values H c i and H C 2, penetrates into the sample pro- 
ducing supercurrent vortices, termed Abrikosov vortices 
( |Abrikosov]|1957| [Blatter et q ?.||1994[). The supercurrent 



superconducting states of great relevance to the problem 



circulates around the normal (i.e. non-superconducting) 
core of the vortex; the core has a size of the order 
the superconducting coherence length £ (parameter of 
a Ginzburg-Landau theory); the supercurrents decay on 
a distance A (London penetration depth) from the core. 
In the following we assume that A ^> £. The circulat- 
ing supercurrents produce magnetic vortices, each char- 
acterized by a total flux equal to a single flux quantum 
$0 = hc/2e (hence the name fluxon for magnetic vortex) 
and a radially decaying magnetic field 



B{r) = B Q K (r/X), 



(39) 



where Kq(x) is a modified Bessel function of zero-th or- 
der and Bo = 3>o/(27rA 2 ). Far from the vortex core, 
r >■ A, the field decays exponentially, B(r)/B a = 
y // Xn/2r e~ r / A , whereas, approaching the vortex center, 
it diverges logarithmically for £ < r <C A, B(r)/B = 
ln(A/r), and then saturates at B(0) = J3 ln(A/£) inside 
the core r < £. As a consequence, the coupling between 
vortex supercurrents and magnetic fields determines a 
long range repulsive vortex-vortex interaction. 

Due to mutual repulsion, fluxons form a lattice (also 
called Abrikosov vortex lattice, usually triangular, pos- 
sibly with defects and dislocations) with average vortex 
density p approximately equal to H/^q. A local den- 
sity I of either a transport current or a magnetization 
current or both, exerts on each fluxon a Lorentz force 
F L = $ I x H/(cff). 

The distribution and microscopic properties of pinning 
centers can qualitatively influence the thermodynamic 
and vortex transport properties of the superconducting 
sample. Pinning forces created by isolated defects in the 
material oppose the motion of the fluxons, thus deter- 
mining a critical (or threshold) current, below which the 
vortex motion is suppressed. Many kinds of artificial pin- 
ning centers have been proposed and developed to in- 
crease the critical current, ranging from the dispersal of 
small non-superconducting second phases to creation of 
defects by irradiation (Altshulcr and Johansen 20041. A 



novel approach to the problem came with advances in 
lithography, which allowed for regular structuring and 
modulation of the sample properties over a large surface 
area (Harada et al. 1996). Long-range correlation in the 



position of the pinning centers resulted in the interplay 
between the length scales characterizing the pin lattice 
and the vortex lattice. Indeed, the magnitude of the com- 
mensuration effects is readily controlled by tuning the in- 
tensity of the magnetic field H : at the first matching field 
Hi = &o jai , pin and vortex lattices are exactly commen- 
surate, with one fluxon per pin lattice cell of area af, at 
the n-th matching field, H(n) — nHi, each pin lattice 
cell is occupied by n — pen fluxons; for irrational H/H\ 
pin and vortex lattices are incommensurate. Commensu- 
ration effects are responsible for a variety of dynamical 



of vortex rectification (Reichhardt et al. 19981 



As a first attempt to design a vortex ratchet |Lee et at] 
( |1999[ ) showed that an ac electric current applied to a 
superconductor patterned with the asymmetric pinning 
potential shown in Fig. |23[ can induce vortex transport 
whose direction is determined only by the asymmetry of 
the pattern. The fluxons were treated as zero-mass point- 
like particles moving on a saw-tooth potential along the 
x-axis, according to the scheme of Sec. |II.D.1| These au- 
thors demonstrated theoretically that, for an appropri- 
ate choice of the pinning potential, such a rocked ratchet 
can be used to manipulate single vortices in supercon- 
ducting samples under realistic conditions. The recti- 
fication mechanisms of an extended overdamped string, 
like a fluxon, on a ratchet potential had been anticipated 



by Marchesoni and coworkers (Costantini et al. 2002 



|Marchesoni 19961. Brownian motor ratcheting of single 



oscillating fluxons have been observed experimentally in 



asymmetrically etched Nb strips (Plourde et al. 20001, 



on planar patterns of columnar defect (Kwok et al. 



2002 1 , and in asymmetric linear arrays of underdamped 



Josephson junctions ( Lee 2003 ) , in an annular Joseph 



son junction (Ustinov et al. 2004 1, to mention just a few 



recent experiments. 
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FIG. 23 (Color online) Diagram of a superconductor in the 
presence of an external magnetic field H. A direct current 
with density / flowing along the y direction (indicated by the 
large arrow) induces a Lorentz force Fl that moves the vor- 
tex in the x direction (upper panel). The superconductor is 
patterned with a pinning potential (lower panel). The poten- 
tial is periodic and asymmetric along the x direction, and is 
invariant under translation along y\ the potential V(x) is an 
effective ratchet potential (see Sec. II. D I. 



A. Fluxon channels 

A more viable approach to achieve controllable 
stochastic transport of fluxons in superconducting de- 
vices consists in exploiting collective boundary effects as 
suggested Sec. |V.A.1| Fluxon rectification in asymmet- 
ric ch annels has been investigated by |Wa mbaug h et al. 
(1999) by means of a molecular dynamics (MD) code, 
originally developed to study magnetic systems with ran- 
dom and correlated pinning (Reichhardt et aT\ 19981. 
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They simulated a thin slice of a zero-field cooled, type- 
II superconducting slab, taken orthogonally to the mag- 
netic vortex lines generated by a tunable external mag- 
netic field of intensity H. On regarding them as fairly 
rigid field structures, magnetic vortices formed a con- 
fined 2D gas of zero-mass repelling particles with density 
p = H/Q . The simulated sample had very strong, effec- 
tively infinite pinning except the "zero pinning" central 
sawtooth-shaped channel, sketched in the inset of Fig. 
|24| In the channel fluxons moved subject to fluxon-fluxon 
repulsions, an externally applied ac driving Lorentz force, 
forces due to thermal fluctuations, and strong damping. 
The fluxons outside the channel could not move, thus 
impeding the movable fluxons in the middle to cross the 
channel walls. As a consequence, the interaction length 
of the fluxon-wall collisions was the same as of the fluxon- 
fluxon collisions. The net fluxon velocity (v) versus tem- 
perature shown in Fig. [24] is negative and exhibits an 
apparent resonant behavior, both versus T and H (i.e., 
the fluxon density p), as anticipated in Sec. V.A.1| 

These conclusions have been recently corroborated ex- 



perimentally by |Togawa et al.\ ( |2005[ ) and then by Yu 
et al. ( 2007 1 . The authors of the more recent work fab- 



ricated triangular channels from bilaycr films of amor- 
phous niobium-germanium, an extremely weak-pinning 
superconductor, and niobium nitride (NbN), with rela- 
tively strong pinning. A reactive ion etching process re- 
moved NbN from regions as narrow as 100 nm, defined 
with electron-beam lithography, to produce weak-pinning 
channels for vortices to move through easily. In contrast, 
vortices trapped in the NbN banks outside of the chan- 
nels remain strongly pinned. The vortex motion through 
such asymmetric channels exhibited interesting asymme- 
tries in both the static depinning and the dynamic flux 
flow. The vortex ratchet effect thus revealed a even richer 
dependence on magnetic field and driving force ampli- 
tu de than anticipated by th e simplified model simulated 
bylWambaugh et al.\<\1999L 



Devices like that simulated in Fig. 24 serve as excellent 
fluxon rectifiers. By coupling two or more such devices, 
one can design fluxon lenses, to regulate the fluxons con- 
centration in chosen regions of a sample, and eventually 
channel networks or circuits, for fluxons in superconduct- 
ing films. 

Fluxon channeling was observed experimentally for the 
first time in an effective ID vortex rectifier by |Villegas| 
et al. (20031. Their device has the pedagogical merit 



of showing all collective transport effects summarized in 
Sec. |V]at work. A 100/im thick Niobiun film was grown 
on an array of nanoscale triangular pinning potentials 
oriented as in Fig. 25 (top panel), with bases lined up 



along the x-axis and vertices pointing in the y-direction. 
Magneto-resistance R{H) experiments were done with a 
magnetic field H applied perpendicular to the substrate 
in a liquid Helium system. The dc magneto-resistance 
curves in the bottom panel of Fig. [25] exhibit commen- 
surability effects in which dissipation minima develop as 
a consequence of the geometrical matching between the 
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FIG. 24 (Color online) Rectified average fluxon velocity (v) 
(in units of A/r) versus T for the sawtooth channel geometry 
in the inset (width: 7A, bottleneck: A). The magnetic field H 
is directed out of the figure. A vertically applied alternating 
current square wave drives the fluxons back and forth horizon- 
tally along the channel with period r = 100 MD steps. The 
number of fluxons per unit cell is: 1 (circles), 25 (diamonds), 
50 (squares), 75 (triangl es), 100 (x), 150 (+), 250 (*). [After 
Wambaugh et al.\\W99\]. 



vortex lattice and the underlying periodic structure. At 
these matching fields, the vortex lattice motion slows 
down, and R(H) minima appear at the equally spaced 
values H(n) = nH x , with H 1 = 32 Oe. The R(H) min- 
ima are sharp at the n — 1,2,3 matching fields, but shal- 
low and not well-defined for higher n. This effect is a 
consequence of the appearance of interstitial vortices at 
increasing H beyond H 3 , three being the maximum num- 
ber of vortices contained in each triangle. 

In order to interpret the experimental results, one sep- 
arates all vortices in two groups: (i) pinned vortices, 
which move from one triangular-shaped pinning trap to 
another one and, thus, are directly affected by an effec- 
tive ID rocked ratchet substrate with positive polarity 
in the y-direction; and (ii) interstitial vortices, which are 
channeled in-between triangles and do not directly in- 



2003 Zhu 



teract with the pinning traps (Savel'ev et al. 

2003a[ ). However, as apparent in the top panel of 



et al. 



Fig. |25|and discussed in Sec. |V.C| pinned vortices, being 
strongly coupled to the substrate potential, determine a 
weaker mean-field potential that acts on the interstitials 
with opposite asymmetry (Savel'ev et aZ.[|2003|. As a 



consequence, when all fluxons are subjected to the same 
ac drive force, the fluxon-fluxon repulsion pushes the in- 
terstitials in the direction opposite to the pinned vortices 
(for an animation see http://dml.riken.go.jp/vortex). 

The effects of this mechanism are illustrated in Fig. 
[26[ At constant temperature (close to T c to avoid ran- 
dom pinning) and H multiple of H±, an ac current den- 
sity I(t) = I ac sm.(Q£) was injected along the x-axis. This 
yields a sinusoidal Lorentz force with amplitude Fl act- 
ing on the vortices along the y-axis. Despite the zero 
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FIG. 25 (Color online) Top panel: Sketch of the positions of the vortices for several matching fields H(n), i.e., for n vortices per 
unit cell. Vortices pinned on the triangles (no more than 3 per unit cell) are shown in red and interstitial vortices in blue. The 
pinning substrate is a square lattice with constant 750nm; the triangular pinning sites have been grown on a Si support, are 
40nm thick with side about 600nm long. Bottom panel: dc magneto-resistance R{H) versus H . The temperature is T = 0.98T C 
and the injected dc current density is Id c = 12.5fc.A- cm -2 . The upper inset shows a micrograph of the measuring bridge, which 
is 40/im wide. The darker central area of the inset is the 90 x 90/im 2 array of magnetic triangles. The lower inset shows the 
characteristics curves Vdc-Idc at the matching magnetic fields H(n) — nHi with Hi = 320e and n = 1 (red), 3 (cyan), 4 (blue), 
6 (green), 8 (magenta), and 10 (black). The V(Idc) curves change abruptly at magnetic fields from n = 3 to n = 4, because 
an ohmic regime appears at low currents. This is a clear signature of the presence of interstitial vortices. [After | Villegas et al 

poo!) ]. 



time average of such Lorentz force, a non-zero dc volt- 
age drop Vdc in the x-direction was measured, thus prov- 
ing that the vortices actually drift in the y-direction. 
Note that, due to the peculiar substrate symmetry, an 
ac current I{t) oscillating along the j/-axis, can rectify 
the fluxon motion only parallel to it (transverse rectifica- 
tion, Sec. VII.B ). The amplitude of the dc voltage signal 



decreases with increasing H because the effective pin- 
ning is suppressed by the inter- vortex repulsion. More- 
over, when n > 3 (corresponding to more than three 
vortices per unit cell), a reversed Vdc signal begins to 
develop with a maximum (marked by blue arrows) oc- 
curring at a lower Lorentz force Fl than the positive dc 
maxima (red arrows): The interstitials, which feel a weak 
inverted ratchet potential with respect to the pinned vor- 
tices, dominate the rectification process. This current 
reversal effect is enhanced when further increasing the 
magnetic field [panels (c) and (d)] and finally, at very 
high magnetic fields, close to the normal state, the volt- 
age reversal, although suppressed in magnitude, spans 
over the entire Fl range [panel (f)]. 



Moshch alkov a nd coworkers ( Silva et al. 2006 deVon- 
del et al\ 2007) also demonstrated longitudinal fluxon 



rectification. Silva et al. (2006) observed multiple cur- 



rent reversals in regular square lattices of asymmetric 
double-well traps periodically driven along an "easy" di- 
rection, as in the simulations by Zhu et al. ( 2003a|b l. 



More interestingly, dc Vondel et al. ( 2007 1 have detected 



inverted fluxon currents in large triangular arrays, similar 
to those reviewed here, but at magnetic fields (or fluxon 
densities) so high that a collective ratchet mechanism 
would not be plausible. They attributed their finding to 



a new intra-antidot rectification mechanism controlled by 
the magnetic flux quantization condition synchronously 
satisfied at the edge of each asymmetric antidot. 



B. Fluxon rectification in 2D arrays 

The rectification of a fluxon lattice, no matter how dis- 
torted or disordered, is an inherently 2D process. Follow- 
ing the theoretical and numerical investigations reviewed 
in Sec . V.B a number of exp e rimen t al groups (Crisan 
~ pOOM; iMenghini et al.\ (120071); IVan Look et al. 



et al. 



(2002 



to mention a few) prototyped devices aimed at 
controlling vortex ratcheting both in the direction paral- 
lel (longitudinal rectifiers) and perpendicular to an ap- 
plied ac Lorentz force (transverse rectifiers). 

Longitudinal transport in a 2D array of asymmetric 
pinning traps has been experimentally obtained, for in- 
stance, by Ide Vondel et al.\ (120051), who reported that, 



under appropriate operating conditions, fluxon rectifica- 
tion can be enhanced for Fl amplitudes comprised be- 
tween two critical pinning forces of the underlying asym- 
metric substrate. The shape of the net dc voltage drop 
Vdc as a function of the drive amplitude indicated that 
their vortex ratchet behaved in a way very different from 
standard overdamped models. Rather, as the repinning 
force necessary to stop vortex motion, is considerably 
smaller than the depinning force, their device resembles 
the behavior of the incrtial ratchets of Sec. III. D. II More 



recently, Shalom and Pastoriza (2005) reported longitu- 



dinal fluxon rectification in square arrays of Josephson 
junctions, where the coupling energies had been period- 
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FIG. 26 (Color online) (a)-(f) Net velocity (v) of vortices ver- 
sus the amplitude Fl of the ac Lorentz force for Q. = 10kHz 
and different matching magnetic fields H(n); other parameter 
values are as in Fig. |25| Red and blue arrows show the di- 
rection of the net flow, dominated respectively, by the pinned 
and the interstitial vortices (see panel (a) of Fig. |25[ ). [After 
Villegas et aT] (j2003|). 



ically modulated along one symmetry axis. Similarly, 
longitudinal fluxon rectification has been observed also 
in square arrays of magnetic antidots with size-graded 
period (Gillijns et al.\ 20071. Both are nice experimental 



implementations of the asymmetric patterns of symmet- 
ric traps introduced in Sec. |V.A.2| Finally, |Ooi et aL\ 
(2007) trapped symmetric fluxon lattices in a triangu- 
lar dot lattice, they had photo-lithographed on a thin 
Bi 2 Sr 2 CaCu20 8+ 5 (Bi2212) single-crystal film, by set- 
ting the applied magnetic field at the lowest multiples 
of Hi . By subjecting the fluxon lattice to a bi-harmonic 
Lorentz force oriented along the crystallographic axes of 
the dot lattice, they obtained a neat demonstration of 
harmonic mixing on a symmetric substrate. However, 
also the results of these experiments can be easily ex- 
plained in terms of simple ID models. 

Transverse transport, instead, requires genuinely irre- 
ducible 2D geometries. When designing and operating a 
device capable of transverse rectification, experimenters 
can vary the orientation of the pinning lattice axes, the 
symmetry axes of the individual traps (if any) and the di- 
rection of the injected current (i.e., of the Lorentz force). 
Many authors explored by numerical simulation the ge- 
ometries best suitable for the experimental implementa- 



tion of this conc ept (|01son Reichhardt and Reichhardt 



2005 Zhu et al. 20011. The most recent realizations of 



transverse fluxon rectifiers fall into two main categories: 



et al. ( 2007 1 modified the experimental set-up of Ville- 



gas et al. ( 2003 1 to investigate the rectification mecha- 
nisms presented in Sees. |V.A|and V.B They confirmed 



the numerical simulations by Savel'ev et al. (2005), who 



had predicted that the same device can exhibit cither 
longitudinal or transverse output current depending on 
orientation of the ac drive with respect to the pinning 
lattice axes. In the longitudinal ratchet configuration a 
sinusoidal driving current I(t) was applied perpendicular 
to the triangle reflection symmetry axis (x-axis) and the 
output voltage signal Vd c was recorded in the same di- 
rection. We recall that the Lorentz force induced vortex 
motion, parallel to the triangle reflection symmetry axis 
(y-axis), corresponds to a voltage drop in the direction 
of the injected current. To observe an optimal transverse 
rectification effect, the axes of the current injection and 
the voltage drop were inverted. The asymmetry of the 
traps with respect to the direction of the Lorentz force is 
responsible for the observed fluxon drift in the y direc- 
tion. 

The number n of vortices per lattice cell was controlled 
by varying the intensity of the magnetic field orthogonal 
to the device, H — nH\. On increasing n, the voltage as- 
sociated to the longitudinal ratchet current changes sign 
and gets amplified, as to be expected due to t he pres - 
ence of the n — 3 interstitials per unit cell (Sec. VILA I. 



The H dependence the transverse fluxon rectification is 
very different: (a) the transverse net current shows no 
inversions; (b) increasing the number of the interstitials 
transverse rectification is suppressed. 

(ii) Asymmetric arrays of symmetric traps. Controlled 
transport of vortices through rows of antidots was mea- 
sured by Wordenweber et al.\ fl2004|) via standard 4- 



probe Hall-type experiments (lAltshuler and Johansen 



20041. 100-150nm thin YBa 2 Cu 3 2 (YBCO) films were 



deposited on Ce0 2 buffered sapphire and then covered 
with a 50nm thick Au layer. Asymmetric lattices of sym- 
metric antidots in the shope of circular micro-holes, were 
patterned via optical lithography and ion beam etching. 
The inset of Fig. 27 shows a typical antidot lattice. A 



tunable dc current J^ c was then injected so that the corre- 
sponding Lorentz force Fx, was at an angle 7 with respect 
to the orientation of the antidot rows. 

Reference measurements on samples without antidots 
as well as temperature-dependent measurements of the 
Hall resistance Rh for different angles 7 clearly indicate 
that for T < 83K and low dc drive current, the Hall resis- 
tance is dominated by the directed vortex motion along 
the antidot rows. Under this circumstances, Rh quanti- 
fies the fluxon current in the direction of the dc current 
Idc (with negative Hall voltage), or equivalently, trans- 
verse to the drive force Fl- In view of the symmetry of 
the traps, for currents perpendicular to the antidot rows, 
no transverse fluxon current was observed; as shown in 
Fig. 27 changing the sign of 7 led to a inversion of 



(i) Symmetric arrays of asymmetric traps. Gonzalez 



the transverse current. Most notably, the detected drift 
was not restricted by a current threshold, that is, fluxons 
seemed not to jump from antidot to antidot, like individ- 
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FIG. 27 Angular dependence of Rh measured on a circular 
shaped 90nm thick sample at T = 30K, H — 143mT and 
with dc current Idc = 10mA. Inset: portion of a 20/imx 10/im 
rectangular lattice of circular antidots with radius of l/im; Idc 
is horizontally oriented from left to right corresponding to a 
Lorentz force pointing upwards. The angular dependence 
is fitted by the sinusoidal law sin 27, as sketched in the inset. 
[After |Wordenweber et al.\ \200ty] . 



ual vortices in a defective superconductor sample (Sec. 
V.A.2 ), but rather obey an Ohmic behavior. The absence 
of an activation threshold was explained in Wordcnweber 
et al. (20041 by invoking a combination of (i) finite size 
effects, as the area where vortices can be nucleated is 
severely restricted by the device geometry; (ii) screening 
effects, as trapped fluxons induce spatially non-uniform 
current distributions around neighboring antidots. Due 
to current screening, the antidot rows serve as easy-flow 
(i.e., threshold- less) paths for the driven vortices, which 
thus acquire a transverse velocity proportional to sin 27 
(see Fig. 27 inset and main panel). 

At higher temperatures and stronger drives, this effect 
becomes negligible and no transverse rectification was de- 

Moreover, on replac- 



2004) 



tected (Wordenweber et al. 
ing Idc with a symmetric ac drive I ac {t)i no net current 
is expected, neither longitudinal nor transversal. 



C. Anisotropic fluxon rectifiers 

Transport control in a binary mixture can be achieved 
in samples with no ratchet substrate. As reported in Sec. 
|V.C[ when one of the components of the binary mixture 
is driven, the moving particles drag along the other non- 
driven component, interacting with it. A time asymmet- 
ric ac drive produces rectification of both components of 
the binary mixture, which can be tuned by means of the 
ac drive parameters. The device acts like a ratchet, but 
it has no ratchet substrate. 



Savel'ev and Nori (20021 proposed to implement this 
new ratchet concept in a layered superconductor. In this 
class of materials, a magnetic field, tilted away from the 



H. 
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FIG. 28 (Color online) Vortex lens. Subjected to a superpo- 
sition of dc and time asymmetric ac in-plane magnetic fields, 
iZjl(t) = H^+H{\ c (t), the JV (horizontal cylinders) are asym- 
metrically pushed in and out of the sample. If the ac com- 
ponent of H\ 1 (t) increases slowly, the PV stacks (pancakes) 
remain trapped on the JV, and both move together towards 
the center of the lens [panel (a)]. If, on the way back, H^ c (t) 
decreases rapidly, the JV leave the PV behind them [(panel 
(b)]. Courtesy of Sergey Savel'ev. 



high-symmetry crystalline c-axis, penetrates the sample 



as two perpendicular vortex arrays (Grigorenko et al. 



2001 Matsuda et al. 2001 



"crossing" vortex lattices (Koshelev 



Ooi et al. 



1999 1, known as 
1999T One vor- 



tex sublattice consists of stacks of pancake vortices (PV) 
aligned along the c-axis, whereas the other is formed by 
Josephson vortices (JV) confined between CuC>2 layers 
(Fig. 28). A weak attractive PV-JV interaction has 



been experimentally observed (Grigorenko et al. 20011. 
The JV are usually very weakly pinned and can be eas- 
ily driven by changing either the in-plane magnetic field, 
or applying an electrical current, I z , flowing along 
the c-axis. 

Time-asymmetrically driven JV can drag along the 
P V, resulting in a net motion of the P V, as in the PV lens 
device illustrated schematically in Fig. |28j The simplest 
possible operating mode consists of slowly increasing the 
in-plane field, H\ | , from to iij" ax for a fixed value of the 
out-of-plane magnetic field, H z . The increasing in- plane 
field slowly drives the JV from the edges to the middle of 
the sample. In turn, the JV drag the PV along with them 
towards the sample center. As a result, asymmetrically 
cycling causes either pumping (focusing) or antipumping 
(defocusing) of the PV at the center of the lens. 

|Cole et~aT. ( |2006 1 actually performed a vortex lens- 
ing experiment on an as-grown single Bi2212 crystal. 
The changes in magnetic induction, arising from PV 
lcnsing/anti-lensing, were detected using one centrally 
placed element of a micro-Hall probe array (Altshuler 



and Johansen 2004 1 . In order to realize the asymmetric 
ac driven mode for the vortex lens, the following steps 
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sketched in Fig. 
was cooled in fixed H z at Hi 



29] were carried out 
II = 0; (2) 



(1) The sample 
A "conditioning" 



triangular wave was run for 4 min to equilibrate the PV 
system; (3) A time-asymmetric, zero-mean ac drive, was 
switched on for 4 min (for an animation sec http: 
dml.riken.go.jp). The magnetic induction, B Zl was then 
monitored in real time starting from step (2) at a cen- 
trally located Hall element and then related to the local 
PV density. The measured efficiency of this vortex lens, 
displayed in Fig. [29] is quite low. However, the same 
experiment, when performed on a symmetric substrate, 
yields a much higher efficiency, as recently proven by |Cole| 
et al.\ ( |2006[ ) for a film of Bi2212 patterned with a square 
array of circular antidots. 

The two main advantages of this class of vortex de- 
vices over the earlier proposals reviewed in Sees. |VII.A| 



and VII.B 



are: (1) the possibility to guide particles with 
no tailored spatial asymmetry; (2) the tuning of the vor- 
tex motion by simply changing the parameters of the ex- 
ternally applied drive. The first feature allows to avoid 
expensive and cumbersome nanofabrication processing; 
the second property becomes very important if the trans- 
port properties of a device must be frequently varied, 
something which is generally hard to achieve in standard 
ratchet devices. 
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FIG. 29 (Color online) Vortex lens operation. The curves 
show the measured percentage change of the magnetic induc- 
tion (PV density) at the samples center, when applying an 
initial conditioning signal followed by pumping/antipumping 
time-asymmetric drives. The conditioned PV density in- 
creases (decreases) during several cycles of the pumping (an- 
tipumping) drive and then saturates. As soon as the drive 
is switched off, the PV density starts to relax from the non- 
equilibrium pumping (antipumping) state towards an equilib- 
rium state. Courtesy of Sergey Savel'ev. 



VIII. QUANTUM DEVICES 

In the foregoing sections, we focused on directed trans- 
port within the realm of classical dynamical laws and 



classical statistical fluctuations. Totally new scenarios 
open up when we try to consistently incorporate quan- 
tum mechanical laws in the operation of an artificial 
Brownian motor. The quantum nature of fluctuations 
and the quantum evolution laws, as governed by quan- 
tum statistical mechanics, allow for unexpected trans- 
port mechanisms. In particular, quantum mechanics pro- 
vides the doorway to new features such as under-barrier(- 
threshold) tunneling, above barrier-reflection and the in- 
terplay of coherent (i.e. with oscillatory relaxation) and 
incoherent (i.e. with overdamped relaxation) quantum 
transport. Last but not least, quantum mechanics gener- 
ates the possibility for non-classical correlations, includ- 
ing entanglement among quantum states in presence of 
coupling. 



A. Quantum dissipative Brownian transport 

To set the stage and in order to elucidate the complex- 
ity involving directed quantum Brownian transport, we 
consider a ID quantum particle of coordinate q and mass 
m, moving in a typically time-dependent rachet-like po- 
tential landscape Vn(q, i). The particle is bi- linearly cou- 
pled with strength Cj to a set of N harmonic oscillators 
Xi, with i = 1, . . . , N. The latter set of oscillators mod- 
els the heat bath, with the oscillators being prepared in a 
canonical state with density matrix pbath, corresponding 
to an isolated, bare bath Hamiltonian and a temperature 
kT. Accordingly, the total dynamics is governed by the 
Hamiltonian: 



H 



P 

2m 



V R {q,t) 



(40) 



N 

E 



2to,- 



1 2 2 
-U>4 x. 



qCiXi 



2m,- cj 2 



The last term, which depends only on the system coor- 
dinate, represents a potential renormalization term which 
is needed to ensure that the potential Vfi(q,t) coincides 
with the bare ratchet potential also in presence of cou- 
pling to the bath degrees of freedom. This Hamiltonian 
has been studied since the early Sixties for systems which 
are weakly coupled to their environment. Only in the 



Eighties it was realized by Caldeira and Leggett ( Caldeira 



and Leggett 1983 1984 1 that this model is also applica- 
ble to strongly damped systems and may be employed to 
describe, for example, dissipative tun neling in solid state 
and chemical physics (Hanggi et al. 1990). 



One may convince oneself that the Hamiltonian in Eq. 
( 40 1 models indeed dissipation. Making use of the so- 



lution of the Heisenberg equations of motion for the ex- 
ternal degrees of freedom, one derives a reduced system 
operator-valued equation of motion, the so-called inertial 
generalized quantum Langevin equation, that is 
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mq\t) + m dsj(t — s)q(s) 
J t 

, dV R (q,t) 



dq 



r](t)-mi(t-t )q(t ). (41) 



The friction kernel is given by 

^) = 7M) = -E^-M), (42) 

»=i 1 1 

and the quantum Brownian force operator reads 

N , 

77(f) = E Cl ( Xifo) cosi^iit - fo]) (43) 
i=i ^ 

H sin(wj f-t ) 

The random force 77(f) is a stationary Gaussian opera- 
tor noise with correlation functions 



S w (f- S )^(77(f)77( S ) + 77(5)77^ 



(44) 
(45) 



E 



2 — ' 771,; 
z=l 



■ COS 



[uJi(t — s)) coth 



^2fcry 



Moreover, being an operator-valued noise, the 77 commu- 
tators do not vanish, i.e., 



[77(f), 77(a)] = -ihY] —t—sin a)) 



i=l 



(46) 



This complex structure for driven quantum Brownian 
motion in a potential landscape follows the fact that a 
consistent description of quantum dissipation necessar- 
ily requires the study of the dynamics in the full Hilbcrt 
space of the system plus environment. This is in clear 
contrast to the classical models, where the Langevin dy- 
namics is directly formulated in the system state space 
( |Hanggi and Ingold |2005[ ). Caution applies in mak- 
ing approximations to this structure, even if done semi- 
classically. The interplay of quantum noise with the com- 
mutator structure is necessary to avoid inconsistencies 
with the thermodynamic laws, such as a spurious finite 
directed current in an equilibrium quantum ratchet, i.e., 
even in absence of rocking, i.e., for Vfl(g, f) = Vn{q) 
( |Hanggi and Ingold| |2005| Machura ef q?.[ |2004a| ). Sim- 
plifications are possible only under specific circumstances 
like in the case of very strong friction, when quantum cor- 
rections can consistently be accounted by a semiclassical 
quantum Smoluchowski equation operating on the state 



space of the classical system, only ( Ankcrhold ef at. 2001 




b) \I(t) =/„ + /, sinffifj 




FIG. 30 (a) Schematic picture of an asymmetric SQUID with 
three junctions threaded by an external magnetic flux, (b) 
Equivalent electric circuit: the two junctions in series of the 
left branch behave like a single junction with u> replaced by 
ip/2. For further details see |Zapata et al. (19961. 



are often in conflict with the laws of thermodynamics; 
this is true in particular for the second law, central to 



this review ( 


Hanggi and Ingold 


2005 


Zueco and Garcia- 


Palacios 


2005 


1. 



An alternative approach to the quantum Langevin de- 
scription, which in addition allows powerful computa- 
tional methods, is based on the real time path integral 
technique. Starting from the quantum statistical repre- 
sentation of the density operator evolution of the total 
dynamics of a system coupled to its environment (s), one 
traces over the bath degrees of freedom to end up with 
a path integral representation for the reduced density 
operator, the so-called influence functional, which con- 
sistently incorporates quantum dissipation (see, for in- 

[1998 " 



stance, (Grifoni and Hanggi 



2005 Hanggi etalX 1990 Kohlcr ef al. 



|Ha nggi and Ingold 
" 2005P ). " 



B. Josephson Brownian motors 

In the following we consider the role of quantum effects 
on artificial Brownian motors made of coupled Josephson 
junctions. 



1. Classical regime 

A first case refers to the limit when quantum coherence 
and tunneling events, such as photon-assisted tunneling, 
can safely be neglected, which is often the case at suffi- 
ciently high temperatures. In this regime the description 
of Josephson quantum systems can well be approximated 
by an effective Fokker-Planck dynamics in the context 



Machura ef al. 2004a 2007b). The situation is presently of the Stewart-McCumber model (Barone and Patcrno 



less settled for weakly damped quantum Brownian mo- 
tors ( Carlo ef al. 2005 Denisov ef al.\ |2008[ ) and for sys- 
tems at high temperatures, where reduced descriptions 



19821, which includes effective parameters containing h. 
An artificial Brownian motor then can be experimentally 
realized for instance with the asymmetric superconduct- 
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FIG. 31 Theoretically evaluated dc current-voltage (Joseph- 
son relation: Vo = (%/2e) (<p(t))) characteristics for the 
SQUID ratchet of Fig. 30 Simulation parameters: (non- 
adiabatic) rocking frequency ui = 2hQ,/eRJi = 0.3, ac drive 
amplitude A = Io/Ji = 1.7, and different noise levels: D = 
ek B T/hJi = (solid), 0.01 (dashed), and 0.5 (dotted). Inset: 
magnified plot showing characteristics steps at fractional val- 
ues of ui for D = 0. Here, R and Ji denote the resistance and 
the critical current amplitude of the two identical Josephson 
junctions placed in the left branch of the device, cf. in Fig. 

EDI 



ing quantum interference device (SQUID) illustrated in 
Fig. |30]( Zapat a et al.\ 1996). In the overdamped regime, 



where capacitative effects can be ignored, such a device 
maps precisely onto the rocked ratchet of Sec. II.D.l 



Under such operating conditions, the phase of the de- 
vice is a classical variable which can be adequately de- 
scribed by the aforementioned "resistively shunted junc- 
tion" model. Thermal Brownian motion at temperature 
T is included by adding Nyquist noise. The behavior of 
the so induced rachet voltage is displayed in Fig. [3T] for a 



non-adiabatically rocked SQUID-ratchet (Zapata et al. 



199G) 



The results for this ratchet setup have been experimen- 
tally validated by use of high temperature superconduct- 
ing dc SQUID's by the Tubingen group (Sterck et al.\ 



2005 2002 Weiss et al. 20001. In the meantime, several 



variants of this scheme have been studied, both theo- 
retically and experimentally, including fluxon ratchets of 



various designs ( 


Beck et al. 2005 


Berger 


2004 Cara- 


pella and Costabile 2001 


Falo et al. 2002| Lee| 2003 


Shalom and Pastoriza 


2005 


Zapata et al. 


19981; cf. also 



Sec. |VII[ Moreover, transport of fluxons in an extended, 
annular Josephson junction has been demonstrated by 
Ustinov et al. ( |2004[ ) also as a harmonic mixing effect. 



2. Quantum regime 

To tackle true quantum effects, including quantum 
tunneling phenomena, we must resort to the theoreti- 
cal scheme outlined in Sect. IVIII.AI A first theoreti- 
cal study of quantum ratchets has been pioneered by 
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FIG. 32 The classical Brownian motor current I c i occurring 
in a dichotomously rocked ratchet potential, Eq. (231, of pe- 



riod L, is compared with the corresponding quantum Brown- 
ian motor current I qm versus dimensionless inverse temper- 
ature T. The current is measured in units of [Lujq] with 
luq = [4tv 2 Vo /(L 2 m)] 5 and temperature is measured in units 
of the (maximal) crossover temperature Tq 1 ^ occurring in the 
rocked potential landscape. Notably, in clear contrast to an 
adiabatically rocked classical Brownian motor, the quantum 
current changes sign near Tq 1 ^ /T = 2.8 which is a manifes- 
tation of true quantum tunneling. For more details we refer 



to the original work by Reimann and Hanggi ( 1998 1 



Reimann et al. (19971 (see also (Reimann and Hanggi 
1998)). For a regime with several (quasi)-bound states 
below a potential barrier one can evaluate the ratchet 
tunneling dynamics in terms of an effective action for the 
extremal bounce-solution (in combination with a fluctu- 
ation analysis around this bounce solution) to obtain the 
corresponding tunneling rates. In Fig. [32] the result of 
such a quantum calculation is compared with the clas- 
sical result for an adiabatically rocked quantum particle 
of mass to and coordinate q dwelling in a ratchet poten- 
tial Vn(q). Note that within a quantum Brownian motor 
scheme the role of mass to, i.e. the inertia, does enter the 
analysis explicitly. The most distinctive quantum signa- 
ture is a striking current reversal which emerges solely 
as a consequence of quantum tunneling under adiabatic 
rocking conditions. Moreover, in contrast to the classical 
result, the directed current no longer vanishes as T tends 
to zero and additional resonance-like features show up. 

Such a dependence of the current versus decreasing 
temperature in a quantum Brownian motor has been cor- 
roborated experimentally for vortices moving in a quasi- 
1D Josephson junction array with ratchet potential pro- 
file specially tailored so as to allow several bands below 
the barrier (Majer et al. 20031. The experimental setup 
of this quantum ratchet device is displayed in Fig. [33] 
The experimental findings are in good agreement with 
the theoretical analysis reported by Grifoni et al. (2002 1. 
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FIG. 33 (Color online) Top left: Scanning electron picture 
(with enlargement on the right) of a strongly asymmetric ar- 
ray of a long, narrow network of Josephson junctions arranged 
in a rectangular lattice. This arrangement determines the po- 
tential shape felt by vortices. Bottom left: Schematic layout. 
The Josephson junctions are represented by crosses, each net- 
work cell being delimited by four junctions. All arrays have a 
length of 303 cells and a width of 5 cells, placed between solid 
superconducting electrodes (busbars). The vortices assume a 
lower energy in cells with larger area and smaller junctions. 
Figure provided by Milena Grifoni, University of Regensburg, 
Germany. 



C. Quantum dot ratchets 

Another ideal resource to observe the interplay of (i) 
quantum Brownian motion, (ii) quantum dissipation, 
and (iii) non-equilibrium driving are semiconductor engi- 
neered quantum rectifiers. Composed of array s of asym- 
metric quantum rectifiers (Linke et al. 1998), these de- 



vices operate on geometric and dynamical length scales 
ranging between a few nanometers up to micrometers. 
This class of devices allowed the first experimental val- 
idation of the distinctive features of quantum rectifiers, 
namely, tunneling enhanced rachet current and tunncl- 
ing induced current reversals, by Linke and collaborators 




FIG. 34 (Color online) In an experimental quantum Brownian 
motor driven by an adiabatic square-wave zero-mean rock- 
ing voltage, quantum tunneling can contribute to the elec- 
tron current. Due to the underlying asymmetric potential 
structure, the two opposite components to the time-averaged 
driven current differ in magnitude, yielding a net quantum 



ratchet current (Linke et al. 1999 Reimann et al. 19971. 



The magnitudes can be individually controlled by tuning the 
temperature. This in turn causes a tunneling induced current 
reversal (occurring near 1.5K in the top graph) that can be 
exploited to direct electrons along a prion designed routes 
(Linke et al. 19991. Below the current versus T graph is a 



trated in Fig. 34 



scanning electron micrograph of the relevant quantum device. 
Figure provided by Heiner Linke, University of Oregon. 



such configurations, one can imagine to operate the de- 
vice as a quantum refrigerator to separate "cold" from 



(Linke et al. 2002, 1999). Their central results are illus- "hot" electrons in the absence of currents (Linke et al. 



Their quantum Brownian motor consisted of a 2D 
gas of electrons moving at the interface of a fabricated, 
ratchet-tailored GaAs/AlGaAs heterostructure. Its op- 
erating regime was achieved by adiabatically switching 
on and off a square wave source-drain periodic voltage of 
frequency 190Hz and amplitude of about lmV. Recall- 
ing that for an adiabatically rocked classical ratchet of 
Sec. II.D.l noise induced transport exhibits no current 
reversals, things change drastically when quantum tun- 
nelling enters into the dynamics. A true benchmark for 
the quantum behavior of an adiabatically rocked Brown- 
ian motor is then the occurrence of a tunneling induced 

1997) . 



2002 1 . At higher temperatures this and other asymmet- 
ric quantum-dot arrays, when subjected to unbiased ac 
drives, exhibit incoherent quantum ratchet currents. Ex- 
perimental such realizations using two-dimensional elec- 
tron systems with broken spatial inversion symmetry are 



reported in Refs. (Lorke et al. 


1998 


Sassine et al. 


2008 


Vidan et al. 


2004 1 . An even richer behavior of the quan- 



turn current, including for example multiple current re 
versals, emerges when this class of devices is operated in 
the non- adiabatic ac regime, as revealed by recent the- 
oretical studies (Goychuk et al. |1998a|b| |Goychuk and 



Hanggi 1998; Goychuk and Hanggi 2005} |Grifoni et al. 



2002 



Kohlcr et al. 



reversal at low temperatures (Reimann et al 



This 

characteristic feature has been first experimentally ver- 
ified with an electron quantum rocked ratchet by |Linke| 

the current reversal reported 
|34] indicates the existence of 



et al.\ ([1999). Moreover 
in the top panel of Fig 
parameter configurations where the current of a quan- 
tum Brownian motor vanishes. In correspondence to 



2005). 

As a second example of an artificial quantum dot Brow- 
nian motor, we consider an experimental double quantum 
dot device where both dots are independently coupled to 
a nonequilibrium energy source which is given by a bi- 
ased quantum point contact. Moreover, both quantum 
dots are embedded in independent electric circuits via a 
common central gate, marked by a "C" in the sketch of 
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FIG. 35 (Color online) (a) AFM micrograph of a double quan- 
tum dot device. Metal surface gates have a light color, black 
squares mark source and drain regions. The black scale bar 
marks the length of lfim. (b) Schematic diagram: A biased 
quantum point contact (QPC) provides the nonequilibrium 
fluctuation source for driving a tunneling current Tdqd in the 
asymmetric double quantum dot (DQD) (dashed two circles) 
The asymmetry is induced via gate 
Figure provided by 



( Khrapai et aZ.||2006 l. 
voltages at the plunger gates 2 and 4 
Vadim S. Khrapai, Ludwig-Maximilian University, Munich 
Germany. 



Fig. |35| On breaking its internal symmetry by tuning 
the dot gate voltages, this double quantum dot acts as a 
quantum ratchet device (Khrapai et al. 2006 1 . For weak 



interdot tunneling, detuning of the quantum dot energy 
levels results in the localization of an electron in one dot, 
so that elastic electron transfer to the other dot is en- 
ergetically forbidden. To operate as a quantum ratchet, 
the device still requires a nonequilibrium energy source. 
To this purpose, a strong tunable bias on the quantum 
point contact induces nonequilibrium energy fluctuations 
across the dividing central gate, thus promoting inelastic 
interdot tunneling inside the Coulomb blockaded dou- 
ble quantum dot. This in turn leads to a net quantum 
ratchet current flow, as plotted in Fig. [36] In the insets of 
the same figure, the interdot tunneling process is sketched 
for the right-to-left transition [m,n + 1] — > [m + l,n], 
with asymmetry energy A = E m+ i_ n — E m , n+ \, and for 
the opposite left-to-right transition with energy —A. No- 
tably, a finite ratchet current is only observed if the elec- 
tron energy states in the dots are detuned asymmetri- 
cally, i.e. when A ^ 0. In contrast, a likely inelastic 
ionization of one dot towards its adjacent lead, followed 
by recharging from the same lead, does not result in a net 
current. The nonequilibrium energy fluctuations, carried 
by the quantum point contact electrons and absorbed by 
the electrons in the double quantum dot, could either 
consist of acoustic phonons, long wavelength photons or 



plasmons (Khrapai et al. 20061 



The technology available to generate 2D electron gases 
can be generalized to ratchet not only charge but also 
spin carriers. The interplay of spatially periodic poten- 
tials, lateral confinement, spin-orbit or Zeeman-type in- 
teractions and ac driving then gives rise to directed spin 
ratchet currents, as theoretically proposed in recent stud- 
ies ( Scheid et al.\ 2007a|b[ ). Of particular experimen- 
tal relevance is the recent theoretically predicted phe- 



nomenon of a spin current which emerges via an unbi- 
ased ac charge current driving a dc spin current in a 
non-magnetic, dissipative spin quantum ratchet which 
is composed of an asymmetric periodic structure with 
Rashba spin-orbit interactions ((Matte] |2008| 



Smirnov 



et al. 20081. Remarkably, this spin current occurs al 
though no magnetic fields are present. 

Likewise, substantial rectification of a spin current can 
also be achieved by coupling impurities to spatially asym- 
metric Luttinger liquids under ac voltage rocking (Brau- 



necker et al. 20071. The transport mechanism in those 
schemes is governed by coherent tunneling processes, 
which will be addressed in more detail in the following 
subsection. 
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FIG. 36 (Color online) Experimental quantum ratchet cur- 
rent Tdqd measured through the double quantum dot as a 
function of its asymmetry energy A. The energy source is a 
nearby quantum point contact biased with Vqpc = —1.55 mV, 
see in Fig. |35| An elastic contribution to Tdqd is subtracted 



(Khrapai et al. 2006). The two insets sketch the correspond- 



ing inelastic tunneling processes which drive the ratchet cur- 
rent. Notably, this quantum ratchet current vanishes when no 
finite zero asymmetry is present. Figure provided by Stefan 
Ludwig, Ludwig-Maximilians-University, Munich, Germany. 



D. Coherent quantum ratchets 

The study of quantum Brownian motors is far from 
being complete. Actually, there is an urgent need for 
further developments. In many cases the role of incoher- 
ent tunneling and coherent transport channels can not 
be clearly separated. For quantum devices built bottom- 
up with engineered molecules or artificial molecules such 
as quantum dots, the role of coherent quantum trans- 
port typically plays a dominant role, as other interactions 
like electron-phonon processes turn out to be negligible. 
The key transport process for quantum ratcheting is then 
driven coherent quantum tunneling ( |Grifoni and Hanggi] 
1998] IKohler et aU|20fJ5b. 
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sical phase space structure ( Denisov and Flach 2001 
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FIG. 37 Time-averaged current / (measured in units of the 
lead coupling strength eF) as a function of the infrared drive 
angular frequency £1 (measured in units of the tunneling ma- 
trix element A = O.leV between sites for a laser amplitude 
A = A (Lehmann et al. 20021. The inset displays the de- 



pendence of I on the static voltage drop, V, applied along 
the molecule. Note that the current at zero voltage is finite, 
thereby indicating a ratchet effect. The driving frequency 
here corresponds to the vertical arrow in the main panel. 



1. Quantum ratchets from molecular wires 

This is the case of the quantum Brownian motors con- 
sisting of nanowires formed by asymmetric molecular 
groups and subjected to infrared light sources ( |Lehman"n 
et al. 2002 1 , or symmetric molecular wires with symme 



try breaking now provided by irradiating the wire with 
harmonic mixing signals (iLehmann et al. 2003 1 . The 



theoretically predicted quantum ratchet current through 
an asymmetric molecular wire irradiated by far-infrared 
light is displayed in Fig. [37} 

In these class of quantum ratchets transport proceeds 
coherently between two or more fermionic leads that en- 
sure the dissipative mechanism for the transported elec- 
trons to relax towards equilibrium Fermi distributions. 
Such a ratchet current can as well be used to sensitively 
probe the role of the electron correlations in the leads, as 
it is the case with a ratchet device coupled to Luttinger 
liquids rather than to Fermi liquids ( |Feldman et aL[|2005| 
Komnik and Gogolin 2003). 



2. Hamiltonian quantum ratchet for cold atoms 

Much to one's surprise, systems governed by strictly 
Hamiltonian dynamics are able to yield directed trans- 
port, as well. Directed transport may then take place 
even though no dissipation mechanisms are present. 
Trivial such examples are provided by integrable regu- 
lar dynamics subjected to unbiased ac-dri ves with fixed 
initial phase, i.e., not to be averaged over (Goychuk and 



Hanggi 2001 Yevtushenko et al., 20001. Much more in- 



triguing are time-dependent driven Hamiltonian quan- 
tum Brownian motors exhibiting a nontrivial, mixed clas- 



Flach et al] |2000| |Schanz et al.\ |2005[ |2001[ ). In these 



systems, the onset of a directed flow requires, apart from 
the necessity of breaking time reversal symmetry, also an 
additional dynamical symmetry breaking. As it turns 
out, in the semi-classical limit, the corresponding di- 
rected flow then obeys a remarkable sum rule: Directed 
currents occurring in regular regimes of the underlying 
phase space dynamics are counter-balanced by directed 
flows occurring in chaotic regimes in phase space ( |Schanz 



et al. 2005| 20011. As a result, if those individual di- 
rected currents are summed up over all disjoint regimes 
in (semi)-classical mixed phase space, no net transport 
emerges. 

The concept of Hamiltonian quantum Brownian mo- 
tors extends as well to fully quantum systems governed 



by a unitary time-evolution (Denisov et al. 


2006 


2007 


Gong et al.\ 2007 Goychuk and Hanggi 2000 


2001 


Schanz et al. 


2005, 


2001 


). In fact, the issue of pure quan- 



turn coherence in the directed transport of chaotic Hamil- 
tonian systems is presently a topic of active research. 
This in particular holds true for cold atoms loaded in 
optical lattices: If properly detuned, the intrinsic quan- 
tum dynamics is then practically dissipation-free, thus 
providing a paradigm for Hamiltonian quantum ratchet 
transport. Various experimental cold atom ratchets have 
been realized (Sec. IV). Relevant to the topic of this 
section is also the demonstration of sawtooth-like, asym- 

All 



(2007) 



metric cold atom potentials by Salger et al. 
these systems are currently being employed to effectively 
shuttle cold matter on a quantum scale. 

In the context of quantum ratchets for cold atoms one 
must distinguish the common case of "rectification of ve- 
locity", implying that the mean position of particles is 
growing linearly in time, from the case with a "rectifica- 
tion of force" , i.e. with mean momentum growing linearly 
in time. The latter cl ass is better classified as "quan tum 
ratchet accelerators" ( Gong and Brumer 2004 20061) . In 



the context of cold atom ratchets the first situation is ob- 
tained by rocking a cold atom gas with temporally asym- 
metric driving forces or temporally asymmetrically flash- 
ing, in combination with spatially asymmetric potential 
kicks (Denisov et al. 20071. The accelerator case origi- 



nates from the physics of quantum r)-kicked cold atoms 
displaying quantum chaos features such as the quantum 
suppression of classical chaotic diffusion (dynamical lo- 
calization) and the diametrically opposite phenomenon 
of quantum resonance (occurring when the kicking period 
is commensurate with the inverse recoil frequency). Un- 
der quantum resonance condition a linear (quadratic) in- 
crease of momentum (kinetic energy) takes place. These 
ratchet acceleration features, theoretica lly predicted in 
suitably modified kicked rotor models (|Kenfack et al.\ 



2008 Lundh and Wallin 2005 Poletti et al. 2007b), re 



quires tuning to exact resonance. In contrast, the accel- 
erator models obtained from either a generalization of 
a quantum kicked rotor or a generalization of a kicked 
Harper model ( |Gong and Brumer 2004 2006 Wang and 
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Gong. 2008) are generic "rectifiers of force" in the sense 
that here no need for tuning to exact resonance is nec- 
essary. Interestingly, these fully quantum ratchet accel- 
erators display unbounded linear growth of mean mo- 
mentum, while the underlying classical dynamics is fully 
chaotic, a situation where classical quantum Brownian 
motor transport necessarily vanishes according to the 



above mentioned classical sum rule ( jSchanz et al. 



2005 



2001). Thus, within this full quantum regime, which car- 
ries no clear relationship with the dynamics in the semi- 
classical regime, the quantum accelerator work b y |Gong| 
and Brumer ( |2006[ ), see also (Wang and Gong 20081, 
presents an intriguing and generic quantum mechanical 
exception to the classical sum rule. 

Early quantum resonance experiments in quantum 
ratchet have already been successfully carried out: For 
a J-kicked rotor model with time-symmetry broken by a 
2-pcriod kicking cycle (asymmetric temporal drive), di- 
rected growth of momentum has been detected by | Jones | 



et al. (20071. For a phase-dependent initial prepara- 
tion of a Bose-Einstein condensate kicked at resonance, 
a momentum acceleration has been observed by |Sad-| 
grove et al. (20071 at zero quasimomentum, while for an 



arbitrary quasimomentum directed quantum Brownian 
transport has been realized in Dana et al. (20081. The 
latter experiment also evidenced that an intrinsic exper- 
imentally non-avoidable finite width in quasimomentum 
causes a suppression of the acceleration eventually lead- 
ing to a saturation effect after short times. 

This field of Hamiltonian quantum ratchets is presently 
undergoing a racy development. For example, it is pos- 
sible to apply control schemes to relative phases for re- 
sulting single-particle quantum Brownian motor currents 
by harvesting Landau-Zener transitions ( Morales-Molina 



et al. 20081. Of particular theoretical and experimental 
interest is also the study of the role of nonlinearity on 
the size of directed quantum currents in interacting cold 
gases, as described within mean field theory by a nonlin- 
ear quantum ratchet evolution of the Gross-Pitaevskii- 
type. There, the interplay of time-dependent driving 
and nonlinear Floquet states yields new features, such 



as lifting of accidental symmetries (Poletti et al. 2007a) 



scenarios of the peristaltic type, or nanosystems which 
are driven by unbiased, but asymmetric mechanical or 
chemical causes that are tightly coupled to resulting mo- 
tion. 



A. Pumping of charge, spin and heat 

A first class of physical systems that comes to mind 
in relation to the working principles of artificial Brow- 
nian motors are nanoscalc pump devices. Pumping is 
characterized by the occurrence of a net flux of particles, 
charges, spins and alike, in response to time dependent 
external manipulations of an otherwise unbiased system. 
This mechanism is well studied and peristaltic pumps 
are being widely exploited in technological applications. 
These systems do not require a periodic arrangement of 
components nor is thermal Brownian motion an issue for 
their operation. In particular, adiabatic turnstiles and 
pumps for charge and other degrees of freedom have at- 



tracted considerable interest both experimentally (Ha 
hberger et al.\ 12001] [Kouwenhoven et al.\ |199T| |Poth 
ier et al.\ |1992| ISwitkes et al.\ |1999 1 and theoretically 



Sinitsyn and Nemenman[|2007 
lessl 119831 IVavilov et all 2001 



( Aono| 


2003 


Brandes and Vorrath] 


2002 


Brouwer 


1998 


Moskalets and Biittiker[ |2004| Shutenko et al. 




2000 



Spivak et al. 



Zhou et al. 



19951 |Thou-| 



1999) 



Realizations of artificial pumps on the nanoscale of- 
ten involve coupled quantum dots or superlattices. Most 
notably, in such peristaltic devices the number of trans- 
ferred charges, or more generally, the number of trans- 
porting units per cycle is directly linked to the cycle 
period. As a result the output current thus becomes 
proportional to the driving frequency. This observation 
leads to the conclusion that high frequency nonadiabatic 
pumping might become more effective. Indeed, nonadia- 



batic pumping of charge, spin ( Flatte . 2008 Scheid et al. 



2007a Smirnov et al. 2008) or also heat (Li et al. 20081 



does exhibit a rich phenomenology, including resonances 



(Arrachea et al. 2007 


Cota et al. 


2005 


Kohler et al. 


2005 


Moskalets and Buttikcr, 2002; 1 


3 latero and Aguado 


2004 


Rey et al. 


2007 


Stafford and Wingreen 


1996 


) and 



( |Morales-Molina and Flach[ |2008 1. 



et al. 



2008; Strass et al. 



and a resonant enhancement of di rected ratchet currents other potentially useful noise-induced features (ISanchez 



20051; hence yielding a close 



interrelation between nonadiabatic pumping in the pres- 
ence of noise and the physics of Brownian motors. 



IX. SUNDRY TOPICS 



In the following we discuss some classes of nanosystems 
and devices which feature directed transport in the spirit 
of Brownian motors. In these systems, however, rectifi- 
cation of Brownian motion does not constitute the main 
element for directed transport. We recall that an artifi- 
cial Brownian motor is mainly noise-controlled, meaning 
that such a motor operates in a hardly predictable man- 
ner. In contrast, we discuss next systems that exhibit 
directed transport predominantly as the result of strong 
coupling schemes. Typical examples are adiabatic pump 



B. Synthetic molecular motors and machines 

As already emphasized in our introduction, the field 
of Brownian motors has its roots in the study and ap- 
plications of intracelluar transport in terms of molecular 



motors ( Astumian and Hanggi 2002 Jiilicher et al. 1997 



Lipowsk y and Klumpp[ |2005| |Reimann 



and Hanggi[|2002 Wang and Oster 2002). These molec 



2002 Rcimann 



ular motors function in view of structure and motility 
by use of specialized proteins in living systems. These 
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FIG. 38 (Color online) Free energy landscape of a single 
azobenzene molecule along the reaction coordinate of the con- 
formation coordinate given by the bond angle &nnc, which 
varies from about —60° to ~ +60°. Transitions can be in- 
duced optically (Ai ~ 420 nm for the short cis-form (0.6 nm) 
and A2 ~ 365 nm for the extended trans- form (0.9 nm), re- 
spectively) from the singlet ground state So into the excited 
singlet state 5*i , from which the molecule subsequently relaxes 
fluctuation-driven into either the cis-form or the trans-form. 
Notably, the trans-form is thermally favored. The insets de- 
pict the corresponding conformations with the transition state 
located near <&nnc ~ 0°. Figure provided by Thorsten Hugel 
from TU-Munich, Germany 



biological motor enzymes are fueled by ATP hydrolysis 
and are able to efficiently perform mechanical work on 



the nanoscale inside biological cell structures ( |Howard 
20frTl|Schliwa||2002l. 



Closely related to synthetic motors are the intrigu- 
ing possibilities of devising DNA-fueled artificial motors: 
several settings render possible to biologically engineer 
nanomachines which move along a priori designed tracks 



2000) 



(Turberfield et al. 2003 Yin et al. 2004 Yorke et al. 



An offspring of this topic is the engineering of nanoma- 
chines based on interlocked molecular species. This has 
spurred a flurry of new investigations within the physical 
biology and the organic and physical chemistry commu- 
nity, aimed at building bottom-up synthetic molecular 
systems which carry out such diverse functions as molec- 
ular switches, molecular rotors and any other kind of 
molecular gears. 

A beautiful example of this class of molecular mo- 
tor is the light powered molecular machine developed by 
Hugel et al. (20021. Bistable photosensitive azobenzene 



molecules can be synthesized into long chains, each con- 
taining molecules in either their trans- or cis-form. The 
free energy landscape corresponding to different molecule 
conformational states is sketched in Fig. [38] versus the 
reaction coordinate. Polymeric chains have the advan- 
tage of scaling up the length changes corresponding to 
the two different conformations of its constituents. The 
length changes of an azobenzene polymeric chains can 
then be transformed into mechanical energy by means of 
the lever arm of a atom force microscope (AFM) which, 
in combination with the light sources, forms the core of 
the light driven molecular motor of Fig. 39 



Such a molecular motor can be operated very much 
according to an idealized thermodynamic cycle of the 
Stirling type: Individual azobenzene polymers are made 
stretch and contract by inducing optical trans- cis transi- 
tions of their constituents. This machine thus demon- 
strates opto-mechanical energy conversion in a single- 
molecule motor device. The analogy with a thermody- 
namic cycle should though not be taken too seriously. A 
thermal Stirling engine operates between baths of differ- 
ing temperature; here, the "expansion" and "compres- 
sion" of the molecular motor are performed at a fixed 
temperature and, most notably, by means of noncquilib- 
rium light sources. This laser operated molecular motor 
can as well be used as a building block to devise a bioin- 
spired molecular locomotive which can be guided forth 
and back on a preassigned track via a laser assisted pro- 
tocol fWang[|2004"|). 



In the operation of these synthetic nanoscale devices, 
especially when powered by light or chemical reactive ad- 
ditives, thermal noise does not necessarily play a major 
role, i.e. the function of these synthetic molecular mo- 
tors is ruled predominantly by deterministic forces, forces 
that depend on the mechanical and chemical properties 



of the molecules (Neuert et al. 2006 1. Yet, this area of re 



search is fascinating and we therefore encourage the inter- 
ested reader to read some recent items and tutorials ( As- 
tumian 2007 Balzani et ai, 2006; Brow ne~and Feringa| 



2006 Porto et al. 2000) and to consult the timely and 



most comprehensive reviews by Kottas et al. ( 2005 ) and 



Kay et al.\ ( |2007[ ) for more details. 



X. CONCLUDING REMARKS 

With this review we have taken the reader on a tour of 
horizon through the many intriguing and multi-facetted 
applications that Brownian motion can offer in the most 
diverse areas of nanotechnology, when combined with 
spatio-temporal symmetry breaking, nonlinearity and, 
possibly, collective interaction effects. 

The physics of classical and quantum Brownian mo- 
tion is by now well established, also in consideration of 
the breadth of the theoretical modeling and experimen- 
tal realizations produced over the last century. Many 
research activities are still spawning in interdisciplinary 
fields encompassing chemistry, biological research, infor- 
mation sciences and even are extending into social sci- 
ences and economics. The main lesson to be learned from 
Robert Brown's and Albert Einstein's work is therefore: 
Rather than fighting thermal motion, we should put it 
to work to our advantage. Brownian motors thrive from 
these ceaseless noise source to efficiently direct, separate, 
pump and shuttle degrees of freedom of differing nature 
reliably and effectively 

In writing this overview we spared no efforts in covering 
a wide range of interesting developments and potential 
achievements. In doing so we nevertheless had to make 
some selective choices of topics and applications, which 
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FIG. 39 (Color online) Operation of a light powered molecular motor: The panel on the left depicts the schematic force 
extension cycle for the opto-mechanical energy conversion cycle of a single poly-azobenzene. In analogy to a thermodynamic 
Stirling-cycle, the polymer is first stretched by the AFM-tip acting as the piston (I — ► II). Then, the application of the first 
optical excitation with A2 — 365nm shortens the polymer by inducing a (poly) trans- (poly) cis transition. This yields a first 
part of work (blue area) in bending the AFM cantilever, which is clearly not as stiff as a piston in a common Stirling motor. 
In analogy to the Stirling cycle, another amount of work is done during the relaxation of the polymer (III — > IV). Finally, a 
second optical excitation with Ai ~ 420nm is needed to reset the molecule into its starting (poly) trans-state (IV — ► I). The 
total work output of the system is the mechanical energy corresponding to the contraction (AL = NAl) of the entire polymer 
chain of TV azobenzenes against the external load. The experimental realization of a full single molecule operating cycle (Hugcl 
et al. 20021 is depicted in the panel on the right. Figure provided by Thorsten Hugel from TU-Munich, Germany 



to some extent reflect the authors' preferences and prej- 
udices. Closely related topics of ongoing research were 
not reviewed for space limits. For example, in Sec. IX. B 
we did not cover in sufficient detail the fascinating topic 
of ATP-driven molecular motors and DNA-fueled mo- 
tors, for which we refer to earlier comprehensive reviews 
and books, like in |Julicher et al. (19971, Howard (20011, 
ScMwal (12002b , and ILipowsky and Klumppl (120051. This 



surely is the case of the bottom-up design and operation 
of synthetic Brownian molecular devices, extensively cov- 
ered by the excellent review of |Kay et al. ( 2007 1 . Another 



such topic is the question of the so-called absolute nega 



in quantum systems ( 


Aguado and Platero 


1997 Grifoni 


and Hanggi, 1998; Keay et aZ.||1995; Platcro and Aguadol 


2004 


) and through mere nonequilibrium driven diffusive 



dynamics in classical systems ( 


Eichhorn et al. 


2002a b 


Kostur et al.\ 2008 Machura 


et al.\ 2007a Ros et al. 



2005) 



Our main focus was on noise-assisted directional trans- 
port, shuttling and pumping of individual or collective 
particle-, charge- or matter-degrees of freedom. The con- 
cept, however, extends as well to the transport of other 
degrees of freedom such as energy (heat) and spin modes 
(see also in Sec. (IX. A I. Both topics are experiencing a 
surge of interest with new exciting achievements being 
reported in the current literature. The concept carries 
potential for yet other applications. Examples that come 
to mind are the noise-assisted directional transport and 
transfer of informational degrees of freedom such as prob- 
ability or entropy and, within a quantum context, the 
shuttling of entanglement information. These issues im- 
mediately relate to the energetics of artificial Brownian 
motors reviewed in Sec. (II.El, including measures of 



open and closed loop control scenarios and other opti- 



mization schemes. 

We have attempted to catch the potential of Brownian 
motors in nowadays nanotechnology by putting them to 
work; we discussed how such motors can be constructed 
and characterized, and how directional, Brownian mo- 
tion driven transport can be controlled, measured, and 
optimized. Moreover, we are confident that the interdis- 
ciplinary style of this overview shall encourage the read- 
ers to bring in new approaches and motivations in this 
challenging and fast growing research area. 
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